compact-genome 0.3.1-alpha.1

Representation of genomes
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

//! A simple representation of a genome as `Vec<u8>`.

use crate::interface::sequence::{
    EditableGenomeSequence, GenomeSequence, GenomeSequenceMut, OwnedGenomeSequence,
};

/// A genome sequence stored as vector of ASCII characters.
pub type AsciiVectorGenome = Vec<u8>;

impl<'a> GenomeSequence<'a, [u8]> for Vec<u8> {}

impl<'a> GenomeSequenceMut<'a, [u8]> for Vec<u8> {}

impl<'a> OwnedGenomeSequence<'a, [u8]> for Vec<u8> {}

impl<'a> EditableGenomeSequence<'a, [u8]> for Vec<u8> {}

#[cfg(test)]
mod tests {
    use crate::implementation::ascii_vec_sequence::AsciiVectorGenome;
    use crate::interface::sequence::{GenomeSequence, OwnedGenomeSequence};

    #[test]
    fn test_reverse_complement() {
        let genome: AsciiVectorGenome = b"ATTCGGT".iter().copied().collect();
        let reverse_complement: AsciiVectorGenome = b"ACCGAAT".iter().copied().collect();
        assert_eq!(genome.reverse_complement(), reverse_complement);
        assert_eq!(genome, reverse_complement.reverse_complement());
    }

    #[test]
    fn test_display() {
        let genome: AsciiVectorGenome = b"ATTCGGT".iter().copied().collect();
        let display_string = genome.as_string();
        let expected_string = "ATTCGGT";
        assert_eq!(display_string, expected_string);
    }
}