rsvart 0.1.5

A small library for representing genomic variants and regions.
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
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111

use crate::genomic::{Contiged, Strand, Stranded};
use crate::ops::{Located, Unit};

pub trait Transposable<C> {
    fn start_on_strand(&self, strand: Strand) -> C;

    fn end_on_strand(&self, strand: Strand) -> C;
}

impl<C, T> Transposable<C> for T
where
    C: Unit,
    T: Contiged<C> + Located<C> + Stranded,
{
    fn start_on_strand(&self, strand: Strand) -> C {
        match self.strand().eq(&strand) {
            true => *self.start(),
            false => *self.contig().end() - *self.end(),
        }
    }

    fn end_on_strand(&self, strand: Strand) -> C {
        match self.strand().eq(&strand) {
            true => *self.end(),
            false => *self.contig().end() - *self.start(),
        }
    }
}

#[cfg(test)]
mod test {
    use super::*;
    use crate::genomic::{AssignedMoleculeType, Contig, SequenceRole};
    use rstest::rstest;

    #[derive(PartialEq, Eq)]
    struct TestRegion {
        strand: Strand,
        start: u8,
        end: u8,
        contig: Contig<u8>,
    }

    impl Stranded for TestRegion {
        fn strand(&self) -> Strand {
            self.strand
        }
    }

    impl Located<u8> for TestRegion {
        fn start(&self) -> &u8 {
            &self.start
        }

        fn end(&self) -> &u8 {
            &self.end
        }
    }

    impl Contiged<u8> for TestRegion {
        type Contig = Contig<u8>;

        fn contig(&self) -> &Self::Contig {
            &self.contig
        }
    }
    fn get_test_contig() -> Contig<u8> {
        Contig::new(
            "1".to_string(),
            SequenceRole::AssembledMolecule,
            "1".to_string(),
            AssignedMoleculeType::Chromosome,
            249,
            "CM000663.1".to_string(),
            "NC_000001.10".to_string(),
            "chr1".to_string(),
        )
        .unwrap()
    }
    #[rstest]
    fn test_start_on_strand() {
        let contig = get_test_contig();
        let test_region = TestRegion {
            strand: Strand::Reverse,
            start: 32,
            end: 200,
            contig,
        };
        assert_eq!(
            test_region.start_on_strand(Strand::Reverse),
            *test_region.start()
        );
        assert_eq!(test_region.start_on_strand(Strand::Forward), 49);
    }

    #[rstest]
    fn test_end_on_strand() {
        let contig = get_test_contig();
        let test_region = TestRegion {
            strand: Strand::Reverse,
            start: 32,
            end: 200,
            contig,
        };
        assert_eq!(
            test_region.end_on_strand(Strand::Reverse),
            *test_region.end()
        );
        assert_eq!(test_region.end_on_strand(Strand::Forward), 217);
    }
}