zoe 0.0.31

A nightly library for viral genomics
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72

use crate::search::fuzzy_substring_match;

#[test]
fn fuzzy_matching_no_diff() {
    let haystack = b"ATAGATTTCTGCTACTCTCCTCATAAGCAGTCCGGTGTATCGAAAGTACA";
    let needle = b"TCCTCATAAG";
    assert_eq!(fuzzy_substring_match(haystack, needle, 0), Some(17));
}

#[test]
fn fuzzy_matching_one_diff() {
    let haystack = b"ATAGATTTCTGCTACTCTCCTCATAAGCAGTCCGGTGTATCGAAAGTACA";
    let needle = b"TCCTTATAAG";
    assert_eq!(fuzzy_substring_match(haystack, needle, 1), Some(17));
}

#[test]
fn fuzzy_matching_two_diff() {
    let haystack = b"ATAGATTTCTGCTACTCTCCTCATAAGCAGTCCGGTGTATCGAAAGTACA";
    let needle = b"TCATTATAAG";
    assert_eq!(fuzzy_substring_match(haystack, needle, 2), Some(17));
}

#[test]
fn fuzzy_matching_long_needle() {
    let haystack = b"ACGT";
    let needle = b"GACTAG";
    assert_eq!(fuzzy_substring_match(haystack, needle, 1), None);
}

#[test]
fn fuzzy_matching_no_haystack() {
    let haystack = b"";
    let needle = b"ACGT";
    assert_eq!(fuzzy_substring_match(haystack, needle, 1), None);
}

#[test]
fn fuzzy_matching_no_needle() {
    let haystack = b"GACTAG";
    let needle = b"";
    assert_eq!(fuzzy_substring_match(haystack, needle, 1), None);
}

#[test]
fn fuzzy_matching_many_diffs() {
    let haystack = b"NNNNNN";
    let needle = b"AC";
    assert_eq!(fuzzy_substring_match(haystack, needle, 2), Some(0));
}

#[test]
fn fuzzy_matching_u8_max() {
    let haystack = vec![b'A'; 256];
    let needle = vec![b'B'; 256];
    assert_eq!(fuzzy_substring_match(&haystack, &needle, u8::MAX), None);
}

#[test]
fn fuzzy_matching_same_len() {
    let haystack = b"TCCTCATAAG";
    let needle = b"TCCTCTTAAG";
    assert_eq!(fuzzy_substring_match(haystack, needle, 1), Some(0));
}

#[test]
fn next_aa_range_search_short_tail() {
    use crate::{prelude::Nucleotides, search::ToRangeSearch};

    let seq: Nucleotides = b"AA".as_slice().into();
    assert_eq!(seq.search_in(1..2).find_next_aa_in_frame(b'M'), None);
}