binseq 0.11.0

A high efficiency binary format for sequencing data
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
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
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
//! VBQ file and block header definitions.
//!
//! `FileHeader` opens the file; a `BlockHeader` precedes each record block.
//! Both are fixed-size (32 bytes) and carry magic numbers for validation.

use std::io::{Read, Write};

use bitnuc_deprec::BitSize;

use crate::{
    error::{HeaderError, ReadError, Result},
    utils::{read_u32_le, read_u64_le},
};

/// Magic bytes at the start of a VBQ file on disk
pub const FILE_MAGIC: [u8; 4] = *b"VSEQ";
const MAGIC: u32 = u32::from_le_bytes(FILE_MAGIC);

/// Block magic number: "BLOCKSEQ" in ASCII (0x5145534B434F4C42)
#[allow(clippy::unreadable_literal)]
const BLOCK_MAGIC: u64 = 0x5145534B434F4C42;

/// Current format version number
const FORMAT: u8 = 1;

/// Size of the file header in bytes
pub const SIZE_HEADER: usize = 32;

/// Size of each block header in bytes
pub const SIZE_BLOCK_HEADER: usize = 32;

/// Default virtual block size in bytes (128KB)
pub const BLOCK_SIZE: u64 = 128 * 1024;

/// Reserved bytes in the file header (placeholder value 42)
pub const RESERVED_BYTES: [u8; 13] = [42; 13];

/// Reserved bytes in block headers (placeholder value 42)
pub const RESERVED_BYTES_BLOCK: [u8; 12] = [42; 12];

#[derive(Default, Debug, Clone, Copy)]
pub struct FileHeaderBuilder {
    qual: Option<bool>,
    block: Option<u64>,
    compressed: Option<bool>,
    paired: Option<bool>,
    bitsize: Option<BitSize>,
    headers: Option<bool>,
    flags: Option<bool>,
}
impl FileHeaderBuilder {
    #[must_use]
    pub fn new() -> Self {
        Self::default()
    }
    #[must_use]
    pub fn qual(mut self, qual: bool) -> Self {
        self.qual = Some(qual);
        self
    }
    #[must_use]
    pub fn block(mut self, block: u64) -> Self {
        self.block = Some(block);
        self
    }
    #[must_use]
    pub fn compressed(mut self, compressed: bool) -> Self {
        self.compressed = Some(compressed);
        self
    }
    #[must_use]
    pub fn paired(mut self, paired: bool) -> Self {
        self.paired = Some(paired);
        self
    }
    #[must_use]
    pub fn bitsize(mut self, bitsize: BitSize) -> Self {
        self.bitsize = Some(bitsize);
        self
    }
    #[must_use]
    pub fn headers(mut self, headers: bool) -> Self {
        self.headers = Some(headers);
        self
    }
    #[must_use]
    pub fn flags(mut self, flags: bool) -> Self {
        self.flags = Some(flags);
        self
    }
    #[must_use]
    pub fn build(self) -> FileHeader {
        FileHeader {
            block: self.block.unwrap_or(BLOCK_SIZE),
            qual: self.qual.unwrap_or(false),
            compressed: self.compressed.unwrap_or(false),
            paired: self.paired.unwrap_or(false),
            bits: self.bitsize.unwrap_or_default(),
            headers: self.headers.unwrap_or(false),
            flags: self.flags.unwrap_or(false),
            ..FileHeader::default()
        }
    }
}

/// 32-byte header at the start of every VBQ file.
#[derive(Clone, Copy, Debug, PartialEq)]
pub struct FileHeader {
    /// Magic number "VSEQ" (4 bytes)
    pub magic: u32,

    /// Format version (1 byte)
    pub format: u8,

    /// Virtual (uncompressed) block size in bytes (8 bytes)
    pub block: u64,

    /// Whether quality scores are included (1 byte)
    pub qual: bool,

    /// Whether blocks are ZSTD compressed (1 byte)
    pub compressed: bool,

    /// Whether records contain paired sequences (1 byte)
    pub paired: bool,

    /// Bits per nucleotide: 2-bit standard or 4-bit ambiguous (1 byte)
    pub bits: BitSize,

    /// Whether length-prefixed sequence headers are included (1 byte)
    pub headers: bool,

    /// Whether per-record flags are included (1 byte)
    pub flags: bool,

    /// Reserved bytes for future extensions (13 bytes)
    pub reserved: [u8; 13],
}
impl Default for FileHeader {
    /// Default block size (128KB), 2-bit encoding, all features disabled.
    fn default() -> Self {
        Self {
            magic: MAGIC,
            format: FORMAT,
            block: BLOCK_SIZE,
            qual: false,
            compressed: false,
            paired: false,
            headers: false,
            flags: false,
            bits: BitSize::default(),
            reserved: RESERVED_BYTES,
        }
    }
}
impl FileHeader {
    /// Parses a header from a 32-byte buffer, validating magic and version.
    pub fn from_bytes(buffer: &[u8; SIZE_HEADER]) -> Result<Self> {
        let magic = read_u32_le(&buffer[0..4]);
        if magic != MAGIC {
            return Err(HeaderError::InvalidMagicNumber(magic).into());
        }
        let format = buffer[4];
        if format != FORMAT {
            return Err(HeaderError::InvalidFormatVersion(format).into());
        }
        let block = read_u64_le(&buffer[5..13]);
        let qual = buffer[13] != 0;
        let compressed = buffer[14] != 0;
        let paired = buffer[15] != 0;
        let bits = match buffer[16] {
            0 | 2 | 42 => BitSize::Two,
            4 => BitSize::Four,
            x => return Err(HeaderError::InvalidBitSize(x).into()),
        };
        let headers = match buffer[17] {
            0 | 42 => false, // backwards compatibility
            _ => true,
        };
        let flags = buffer[18] != 0;
        let Ok(reserved) = buffer[19..32].try_into() else {
            return Err(HeaderError::InvalidReservedBytes.into());
        };
        Ok(Self {
            magic,
            format,
            block,
            qual,
            compressed,
            paired,
            bits,
            headers,
            flags,
            reserved,
        })
    }

    /// Serializes the header as 32 bytes to a writer.
    pub fn write_bytes<W: Write>(&self, writer: &mut W) -> Result<()> {
        let mut buffer = [0u8; SIZE_HEADER];
        buffer[0..4].copy_from_slice(&self.magic.to_le_bytes());
        buffer[4] = self.format;
        buffer[5..13].copy_from_slice(&self.block.to_le_bytes());
        buffer[13] = self.qual.into();
        buffer[14] = self.compressed.into();
        buffer[15] = self.paired.into();
        buffer[16] = self.bits.into();
        buffer[17] = self.headers.into();
        buffer[18] = self.flags.into();
        buffer[19..32].copy_from_slice(&self.reserved);
        writer.write_all(&buffer)?;
        Ok(())
    }

    /// Reads and parses a 32-byte header from a reader.
    pub fn from_reader<R: Read>(reader: &mut R) -> Result<Self> {
        let mut buffer = [0u8; SIZE_HEADER];
        reader.read_exact(&mut buffer)?;
        Self::from_bytes(&buffer)
    }

    #[must_use]
    pub fn is_paired(&self) -> bool {
        self.paired
    }
}

/// 32-byte header preceding each block in a VBQ file.
#[derive(Clone, Copy, Debug)]
pub struct BlockHeader {
    /// Magic number "BLOCKSEQ" (8 bytes)
    pub magic: u64,

    /// Actual on-disk block size in bytes; differs from the virtual size when compressed (8 bytes)
    pub size: u64,

    /// Number of records in this block (4 bytes)
    pub records: u32,

    /// Reserved bytes for future extensions (12 bytes)
    pub reserved: [u8; 12],
}
impl BlockHeader {
    /// Creates a new block header with the given on-disk size and record count.
    #[must_use]
    pub fn new(size: u64, records: u32) -> Self {
        Self {
            magic: BLOCK_MAGIC,
            size,
            records,
            reserved: RESERVED_BYTES_BLOCK,
        }
    }

    #[must_use]
    pub fn empty() -> Self {
        Self {
            magic: BLOCK_MAGIC,
            size: 0,
            records: 0,
            reserved: RESERVED_BYTES_BLOCK,
        }
    }

    #[must_use]
    pub fn is_empty(&self) -> bool {
        self.size == 0 && self.records == 0
    }

    /// Serializes the block header as 32 bytes to a writer.
    pub fn write_bytes<W: Write>(&self, writer: &mut W) -> Result<()> {
        let mut buffer = [0u8; SIZE_BLOCK_HEADER];
        buffer[0..8].copy_from_slice(&self.magic.to_le_bytes());
        buffer[8..16].copy_from_slice(&self.size.to_le_bytes());
        buffer[16..20].copy_from_slice(&self.records.to_le_bytes());
        buffer[20..].copy_from_slice(&self.reserved);
        writer.write_all(&buffer)?;
        Ok(())
    }

    /// Parses a block header from a 32-byte buffer, validating the magic number.
    pub fn from_bytes(buffer: &[u8; SIZE_BLOCK_HEADER]) -> Result<Self> {
        let magic = read_u64_le(&buffer[0..8]);
        if magic != BLOCK_MAGIC {
            return Err(ReadError::InvalidBlockMagicNumber(magic, 0).into());
        }
        let size = read_u64_le(&buffer[8..16]);
        let records = read_u32_le(&buffer[16..20]);
        Ok(Self::new(size, records))
    }

    #[must_use]
    pub fn size_with_header(&self) -> usize {
        self.size as usize + SIZE_BLOCK_HEADER
    }
}

#[cfg(test)]
mod tests {
    use super::*;

    // ==================== FileHeaderBuilder Tests ====================

    #[test]
    fn test_builder_block_and_bitsize() {
        let header = FileHeaderBuilder::new()
            .block(4096)
            .bitsize(BitSize::Four)
            .build();
        assert_eq!(header.block, 4096);
        assert_eq!(header.bits, BitSize::Four);
    }

    #[test]
    fn test_builder_defaults() {
        let header = FileHeaderBuilder::new().build();
        assert_eq!(header.block, BLOCK_SIZE);
        assert!(!header.qual);
        assert!(!header.compressed);
        assert!(!header.paired);
        assert!(!header.headers);
        assert!(!header.flags);
    }

    // ==================== FileHeader Constructor Tests ====================

    #[test]
    fn test_file_header_default() {
        let header = FileHeader::default();
        assert_eq!(header.block, BLOCK_SIZE);
        assert!(!header.is_paired());
    }

    // ==================== FileHeader from_bytes/from_reader Tests ====================

    #[test]
    fn test_file_header_roundtrip() {
        let header = FileHeaderBuilder::new()
            .qual(true)
            .paired(true)
            .headers(true)
            .flags(true)
            .build();
        let mut buffer = Vec::new();
        header.write_bytes(&mut buffer).unwrap();
        let mut cursor = std::io::Cursor::new(buffer);
        let parsed = FileHeader::from_reader(&mut cursor).unwrap();
        assert_eq!(parsed, header);
    }

    #[test]
    fn test_file_header_from_bytes_four_bit() {
        let header = FileHeaderBuilder::new().bitsize(BitSize::Four).build();
        let mut buffer = [0u8; SIZE_HEADER];
        {
            let mut cursor = std::io::Cursor::new(&mut buffer[..]);
            header.write_bytes(&mut cursor).unwrap();
        }
        let parsed = FileHeader::from_bytes(&buffer).unwrap();
        assert_eq!(parsed.bits, BitSize::Four);
    }

    #[test]
    fn test_file_header_from_bytes_invalid_magic() {
        let buffer = [0u8; SIZE_HEADER];
        let result = FileHeader::from_bytes(&buffer);
        assert!(result.is_err());
    }

    #[test]
    fn test_file_header_from_bytes_invalid_format() {
        let header = FileHeader::default();
        let mut buffer = [0u8; SIZE_HEADER];
        {
            let mut cursor = std::io::Cursor::new(&mut buffer[..]);
            header.write_bytes(&mut cursor).unwrap();
        }
        buffer[4] = 99;
        let result = FileHeader::from_bytes(&buffer);
        assert!(result.is_err());
    }

    #[test]
    fn test_file_header_from_bytes_invalid_bitsize() {
        let header = FileHeader::default();
        let mut buffer = [0u8; SIZE_HEADER];
        {
            let mut cursor = std::io::Cursor::new(&mut buffer[..]);
            header.write_bytes(&mut cursor).unwrap();
        }
        buffer[16] = 99;
        let result = FileHeader::from_bytes(&buffer);
        assert!(result.is_err());
    }

    #[test]
    fn test_file_header_from_reader_truncated() {
        let mut cursor = std::io::Cursor::new(vec![0u8; 5]);
        let result = FileHeader::from_reader(&mut cursor);
        assert!(result.is_err());
    }

    // ==================== BlockHeader Tests ====================

    #[test]
    fn test_block_header_from_bytes_invalid_magic() {
        let buffer = [0u8; SIZE_BLOCK_HEADER];
        let result = BlockHeader::from_bytes(&buffer);
        assert!(result.is_err());
    }

    #[test]
    fn test_block_header_roundtrip() {
        let header = BlockHeader::new(2048, 42);
        let mut buffer = Vec::new();
        header.write_bytes(&mut buffer).unwrap();
        let mut fixed = [0u8; SIZE_BLOCK_HEADER];
        fixed.copy_from_slice(&buffer);
        let parsed = BlockHeader::from_bytes(&fixed).unwrap();
        assert_eq!(parsed.size, 2048);
        assert_eq!(parsed.records, 42);
        assert!(!parsed.is_empty());
    }
}