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,));