zoe 0.0.30

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

use crate::data::{
    ByteValidator,
    mappings::ByteMap,
    types::{amino_acids::AminoAcids, nucleotides::Nucleotides},
};

/// Provides methods for validating and transforming sequence data using byte
/// mappings
pub trait RetainSequence {
    /// Retains only bytes that are marked as valid in the validation mapping
    fn retain_by_validation(&mut self, validation_mapping: &ByteValidator);
    /// Retains and transforms bytes using the transformation mapping, removing
    /// any that map to 0
    fn retain_by_recoding(&mut self, transformation_mapping: &ByteMap);
}

impl RetainSequence for Vec<u8> {
    /// Allows for filtering of biological sequences using a validation mapping.
    #[inline]
    fn retain_by_validation(&mut self, validation_mapping: &ByteValidator) {
        self.retain(|b| validation_mapping[*b]);
    }

    /// Allows for filtering of biological sequences using byte re-encoding. The
    /// 0-byte is assumed to be an invalid pattern.
    #[inline]
    fn retain_by_recoding(&mut self, transformation_mapping: &ByteMap) {
        self.retain_mut(|b| {
            *b = transformation_mapping[*b];
            *b > 0
        });
    }
}

/// Provides pass through functionality for the standard library.
pub trait StdForSequences: AsMut<Vec<u8>> {
    #[inline]
    fn retain<F>(&mut self, f: F)
    where
        F: FnMut(&u8) -> bool, {
        self.as_mut().retain(f);
    }
}
impl StdForSequences for Nucleotides {}
impl StdForSequences for AminoAcids {}