let to_channel ?(chr_map=ident) ?(increment_bp=0) t cout = 
  let f (s,i,v) = fprintf cout "%s\t%d\t%f\n" (chr_map s) (i+increment_bp) v in
  List.iter ~f t