Function intspan::get_consensus_poa

source ·
pub fn get_consensus_poa(seqs: &[&[u8]]) -> Result<String>
Expand description
match which::which("spoa") {
    Ok(_) => {
        let seqs = vec![
        //              *
            b"TTAGCCGCTGAGAAGC".as_ref(),
            b"TTAGCCGCTGAGAAGC".as_ref(),
            b"TTAGCCGCTGA-AAGC".as_ref(),
        ];
        let cons = intspan::get_consensus_poa(&seqs).unwrap();
        assert_eq!(cons, "TTAGCCGCTGAGAAGC".to_string());

        let seqs = vec![
        //      *   **
            b"AAATTTTGG".as_ref(),
            b"AAAATTTTT".as_ref(),
        ];
        let cons = intspan::get_consensus_poa(&seqs).unwrap();
        assert_eq!(cons, "AAAATTTTGG".to_string());

        let seqs = vec![
        //           **
            b"AAAATTTTGG".as_ref(),
            b"AAAATTTTTG".as_ref(),
        ];
        let cons = intspan::get_consensus_poa(&seqs).unwrap();
        assert_eq!(cons, "AAAATTTTTG".to_string());

        let seqs = vec![
        //
            b"AAAATTTTGG".as_ref(),
        ];
        let cons = intspan::get_consensus_poa(&seqs).unwrap();
        assert_eq!(cons, "AAAATTTTGG".to_string());

    }
    Err(_) => {}
}