mini_myers 0.1.0

SIMD myers for short queries
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

use mini_myers::{mini_search_with_positions, TQueries};

#[allow(dead_code)]
fn generate_random_dna(l: usize) -> Vec<u8> {
    use rand::rngs::StdRng;
    use rand::{Rng, SeedableRng};
    const DNA: &[u8; 4] = b"ACGT";
    let mut rng = StdRng::seed_from_u64(42);
    let mut dna = Vec::with_capacity(l);
    for _ in 0..l {
        let idx = rng.gen_range(0..4);
        dna.push(DNA[idx]);
    }
    dna
}

fn main() {
    let queries = vec![b"ATG".to_vec(), b"TTG".to_vec()];
    let transposed = TQueries::new(&queries);
    let target = b"CCCTCGCCCCCCATGCCCCC";
    let mut results = Vec::new();
    mini_search_with_positions(&transposed, target, 1, None, &mut results);
    println!("Result: {:?}", results);
}