Function intspan::alignment_stat
source · pub fn alignment_stat(seqs: &[&[u8]]) -> (i32, i32, i32, i32, i32, f32)
Expand description
Basic stats on alignments
let seqs = vec![
// *
b"AAAATTTTGG".as_ref(),
b"aaaatttttg".as_ref(),
];
assert_eq!(intspan::alignment_stat(&seqs), (10, 10, 1, 0, 0, 0.1,));
let seqs = vec![
//* * *
b"TTAGCCGCTGAGAAGCC".as_ref(),
b"GTAGCCGCTGA-AGGCC".as_ref(),
];
assert_eq!(intspan::alignment_stat(&seqs), (17, 16, 2, 1, 0, 0.125,));
let seqs = vec![
// * ** * ** * *
b"GATTATCATCACCCCAGCCACATW".as_ref(),
b"GATTTT--TCACTCCATTCGCATA".as_ref(),
];
assert_eq!(intspan::alignment_stat(&seqs), (24, 21, 5, 2, 1, 0.238,));