fgumi 0.1.3

High-performance tools for UMI-tagged sequencing data: extraction, grouping, and consensus calling
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

pub mod array;
mod hex;
mod string;

use std::io;

use noodles::sam::alignment::record::data::field::Value;

use self::{array::write_array, hex::write_hex, string::write_string};
use crate::vendored::bam_codec::encoder::num::{
    write_f32_le, write_i8, write_i16_le, write_i32_le, write_u8, write_u16_le, write_u32_le,
};

pub fn write_value(dst: &mut Vec<u8>, value: &Value) -> io::Result<()> {
    match value {
        Value::Character(c) => write_u8(dst, *c),
        Value::Int8(n) => write_i8(dst, *n),
        Value::UInt8(n) => write_u8(dst, *n),
        Value::Int16(n) => write_i16_le(dst, *n),
        Value::UInt16(n) => write_u16_le(dst, *n),
        Value::Int32(n) => write_i32_le(dst, *n),
        Value::UInt32(n) => write_u32_le(dst, *n),
        Value::Float(n) => write_f32_le(dst, *n),
        Value::String(s) => write_string(dst, s)?,
        Value::Hex(s) => write_hex(dst, s)?,
        Value::Array(array) => write_array(dst, array)?,
    }

    Ok(())
}

#[cfg(test)]
mod tests {
    use bstr::ByteSlice;

    use super::*;

    #[test]
    fn test_write_value() -> Result<(), Box<dyn std::error::Error>> {
        use noodles::sam::alignment::record::data::field::value::Array;

        fn t(buf: &mut Vec<u8>, value: &Value, expected: &[u8]) -> io::Result<()> {
            buf.clear();
            write_value(buf, value)?;
            assert_eq!(buf, expected);
            Ok(())
        }

        use super::array::tests::T;

        let mut buf = Vec::new();

        t(&mut buf, &Value::Character(b'n'), b"n")?;
        t(&mut buf, &Value::Int8(1), &[0x01])?;
        t(&mut buf, &Value::UInt8(2), &[0x02])?;
        t(&mut buf, &Value::Int16(3), &[0x03, 0x00])?;
        t(&mut buf, &Value::UInt16(5), &[0x05, 0x00])?;
        t(&mut buf, &Value::Int32(8), &[0x08, 0x00, 0x00, 0x00])?;
        t(&mut buf, &Value::UInt32(13), &[0x0d, 0x00, 0x00, 0x00])?;
        t(&mut buf, &Value::Float(8.0), &[0x00, 0x00, 0x00, 0x41])?;

        t(&mut buf, &Value::String(b"ndls".as_bstr()), &[b'n', b'd', b'l', b's', 0x00])?;

        t(&mut buf, &Value::Hex(b"CAFE".as_bstr()), &[b'C', b'A', b'F', b'E', 0x00])?;

        t(
            &mut buf,
            &Value::Array(Array::UInt8(Box::new(T::new(&[0])))),
            &[
                b'C', // subtype = UInt8
                0x01, 0x00, 0x00, 0x00, // count = 1
                0x00, // values[0] = 0
            ],
        )?;

        Ok(())
    }
}