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

use std::ops::Range;

pub type Position = u64;
pub type Length = u64;

pub trait AbstractInterval {
    /// Identifier for a genomic contig, e.g., a chromosome
    fn contig(&self) -> &str;
    /// Interval on the contig
    fn range(&self) -> Range<Position>;
}

#[derive(new, Debug, PartialEq, Eq, Clone, Hash, Serialize, Deserialize)]
pub struct Interval {
    contig: String,
    range: Range<Position>,
}

impl Interval {
    /// Mutable reference to interval on the contig
    pub fn range_mut(&mut self) -> &mut Range<Position> {
        &mut self.range
    }
}

impl AbstractInterval for Interval {
    fn contig(&self) -> &str {
        &self.contig
    }

    fn range(&self) -> Range<Position> {
        self.range.clone()
    }
}

pub trait AbstractLocus {
    /// Identifier for a genomic contig, e.g., a chromosome
    fn contig(&self) -> &str;
    /// Position on the contig
    fn pos(&self) -> Position;
}

#[derive(new, Debug, PartialEq, Eq, Clone, Hash, Serialize, Deserialize)]
pub struct Locus {
    contig: String,
    pos: Position,
}

impl Locus {
    /// Mutable reference to position.
    fn pos_mut(&mut self) -> &mut Position {
        &mut self.pos
    }
}

impl AbstractLocus for Locus {
    fn contig(&self) -> &str {
        &self.contig
    }

    fn pos(&self) -> Position {
        self.pos
    }
}