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
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136

use std::{
    fmt::Display,
    io::{Result, Write},
    ops::{Deref, DerefMut},
    str::FromStr,
};

use crate::{
    property::{Named, Parsable, RegionCore, Scored, Serializable, Strand, Stranded},
    ChrRef,
};

use super::{Bed5, RcCowString, ToSelfContained};

#[derive(Clone, PartialEq, Eq, PartialOrd, Ord)]
pub struct Bed6<'a, T = f64> {
    inner: Bed5<'a, T>,
    pub strand: Strand,
}

impl<'a, T> Deref for Bed6<'a, T> {
    type Target = Bed5<'a, T>;

    fn deref(&self) -> &Self::Target {
        &self.inner
    }
}

impl<'a> DerefMut for Bed6<'a> {
    fn deref_mut(&mut self) -> &mut Self::Target {
        &mut self.inner
    }
}

impl<'a, T: Display> Serializable for Bed6<'a, T> {
    fn dump<W: Write>(&self, mut fp: W) -> Result<()> {
        self.inner.dump(&mut fp)?;
        write!(fp, "\t{}", self.strand)?;
        Ok(())
    }
}

impl<'a, T: Display> Serializable for Option<Bed6<'a, T>> {
    fn dump<W: Write>(&self, mut fp: W) -> Result<()> {
        if let Some(inner) = self {
            inner.dump(fp)
        } else {
            fp.write_all(b".\t.\t.\t.\t.\t.")
        }
    }
}

impl<'a, T: FromStr> Parsable<'a> for Bed6<'a, T> {
    fn parse(s: &'a str) -> Option<(Self, usize)> {
        let (inner, mut start) = Bed5::parse(s)?;
        if s[start..].starts_with('\t') {
            start += 1;
        }
        let s = &s[start..];
        let brk = memchr::memchr(b'\t', s.as_bytes()).unwrap_or(s.len());
        let strand = match &s[..brk] {
            "+" => Strand::Positive,
            "-" => Strand::Negative,
            _ => Strand::Unknown,
        };

        Some((Self { inner, strand }, start + brk))
    }
}

impl<'a, S> Bed6<'a, S> {
    pub fn new<T: RegionCore + Named<'a> + Scored<S> + Stranded>(region: &T) -> Self
    where
        S: Default,
    {
        let inner = Bed5::new(region);
        let strand = region.strand();
        Self { inner, strand }
    }

    #[inline(always)]
    pub fn set_strand(&mut self, strand: &str) {
        match strand {
            "+" => self.strand = Strand::Positive,
            "-" => self.strand = Strand::Negative,
            _ => self.strand = Strand::Unknown,
        }
    }
}

impl<'a, S> RegionCore for Bed6<'a, S> {
    #[inline(always)]
    fn start(&self) -> u32 {
        self.inner.start()
    }
    #[inline(always)]
    fn end(&self) -> u32 {
        self.inner.end()
    }
    #[inline(always)]
    fn chrom(&self) -> ChrRef<'static> {
        self.inner.chrom()
    }
}

impl<'a, T: Clone> Scored<T> for Bed6<'a, T> {
    #[inline(always)]
    fn score(&self) -> Option<T> {
        self.inner.score()
    }
}

impl<'a, T> Stranded for Bed6<'a, T> {
    fn strand(&self) -> Strand {
        self.strand
    }
}

impl<'a, T> Named<'a> for Bed6<'a, T> {
    fn name(&self) -> &str {
        self.inner.name()
    }
    fn to_cow(&self) -> RcCowString<'a> {
        self.inner.to_cow()
    }
}

impl <'a> ToSelfContained for Bed6<'a> {
    type SelfContained = Bed6<'static>;
    fn to_self_contained(&self) -> Self::SelfContained {
        Bed6 {
            inner: self.inner.to_self_contained(),
            strand: self.strand,
        }
    }
}