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

use std::io;

use noodles_sam as sam;

/// A raw BAM record sequence.
#[derive(Debug, Eq, PartialEq)]
pub struct Sequence<'a> {
    src: &'a [u8],
    base_count: usize,
}

impl<'a> Sequence<'a> {
    pub(super) fn new(src: &'a [u8], base_count: usize) -> Self {
        Self { src, base_count }
    }

    /// Returns whether there are any bases.
    pub fn is_empty(&self) -> bool {
        self.src.is_empty()
    }

    /// Returns the number of bases in the sequence.
    ///
    /// This is _not_ the length of the buffer.
    pub fn len(&self) -> usize {
        self.base_count
    }
}

impl<'a> AsRef<[u8]> for Sequence<'a> {
    fn as_ref(&self) -> &[u8] {
        self.src
    }
}

impl<'a> TryFrom<Sequence<'a>> for sam::record::Sequence {
    type Error = io::Error;

    fn try_from(bam_sequence: Sequence<'a>) -> Result<Self, Self::Error> {
        use crate::reader::record::get_sequence;

        let mut src = bam_sequence.src;
        let mut sam_sequence = Self::default();
        let base_count = bam_sequence.base_count;
        get_sequence(&mut src, &mut sam_sequence, base_count)?;

        Ok(sam_sequence)
    }
}