Function parasailors::local_alignment_stats

pub fn local_alignment_stats(
    query_sequence: &[u8], 
    database_sequence: &[u8], 
    open_cost: i32, 
    gap_extend_cost: i32, 
    substitution_matrix: &Matrix
) -> AlignmentStats

Provides statistics for local pairwise alignment using a vectorized algorithm.

Examples

let identity_matrix = Matrix::new(MatrixType::Identity);
let query = b"AAAACCCCCCCCCCGGG";

let reference = b"AAAAAAAAAACCCCCCCCCCGGGGGGGGGGTTTTTTTTTTTNNNNNNNNN";
let stats = local_alignment_stats(query, reference, 1, 1, &identity_matrix);
assert_eq!(17, stats.score);
assert_eq!(17, stats.num_matches);
assert_eq!(17, stats.num_positive_subs);
assert_eq!(17, stats.align_length);
assert_eq!(17, stats.query_end);
assert_eq!(23, stats.ref_end);