kira-ls-aligner 0.1.2

Unified short- and long-read sequence aligner written in Rust 2024. It combines minimap2-style minimizers and chaining with BWA-MEM2-style exact-match anchoring and output semantics. The goal is drop-in compatibility with bwa-mem pipelines while supporting long reads efficiently.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25

#[inline]
fn comp_base(b: u8) -> u8 {
    match b {
        b'A' | b'a' => b'T',
        b'C' | b'c' => b'G',
        b'G' | b'g' => b'C',
        b'T' | b't' => b'A',
        _ => b'N',
    }
}

/// Reverse-complement DNA sequence into a new owned buffer.
pub fn reverse_complement(seq: &[u8]) -> Vec<u8> {
    let mut out = Vec::with_capacity(seq.len());
    out.extend(seq.iter().rev().map(|&b| comp_base(b)));
    out
}

/// Reverse-complement into a caller-provided scratch buffer.
pub fn reverse_complement_into<'a>(seq: &[u8], out: &'a mut Vec<u8>) -> &'a [u8] {
    out.clear();
    out.reserve(seq.len());
    out.extend(seq.iter().rev().map(|&b| comp_base(b)));
    out.as_slice()
}