needletail 0.7.3

FASTX parsing and k-mer methods
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

use needletail::parse_fastx_file;

const TEST_FILES: [&str; 4] = [
    "./tests/data/test.fa.gz",
    "./tests/data/test.fa.bz2",
    "./tests/data/test.fa.xz",
    "./tests/data/test.fa.zst",
];

#[cfg(feature = "compression")]
#[test]
fn can_read_compressed_files_automatically() {
    use needletail::parser::Format;
    for p in &TEST_FILES {
        let mut reader = parse_fastx_file(p).unwrap();
        let mut i = 0;
        while let Some(record) = reader.next() {
            let seq = record.unwrap();
            assert_eq!(seq.format(), Format::Fasta);

            match i {
                0 => {
                    assert_eq!(seq.id(), b"test");
                    assert_eq!(seq.raw_seq(), b"AGCTGATCGA");
                    assert_eq!(seq.qual(), None);
                }
                1 => {
                    assert_eq!(seq.id(), b"test2");
                    assert_eq!(seq.raw_seq(), b"TAGC");
                    assert_eq!(seq.qual(), None);
                }
                _ => unreachable!("Too many records"),
            }
            i += 1;
        }
        assert_eq!(i, 2);
    }
}

#[cfg(not(feature = "compression"))]
#[test]
fn errors_on_compressed_files() {
    for p in &TEST_FILES {
        assert!(parse_fastx_file(p).is_err());
    }
}