Function rustybam::nucfreq::nucfreq

pub fn nucfreq(bam: &mut IndexedReader, rgn: &Region) -> Vec<Nucfreq>

get count for A,C,G,T at every pos

Example

use rust_htslib::{bam, bam::Read};
let mut bam = bam::IndexedReader::from_path(".test/test_nucfreq.bam").unwrap();
let rgn = rustybam::bed::Region {
    name : "CHROMOSOME_I".to_string(),
    st :  1,
    en : 102,
    id : "None".to_string()
};
let vec = rustybam::nucfreq::nucfreq( &mut bam, &rgn);
eprintln!("{:?}", vec[0].a);
for f in vec {
  let t = vec![f.a, f.c, f.g, f.t];
  let max = t.iter().max().unwrap();
  if(*max != 0){
    assert_eq!(*max, 2);
  }
}