fastseq 0.1.4

fastseq is a Rust crate for parsing FASTA/Q sequences
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

## fastseq

fastseq is a high-performance Rust crate for parsing FASTA/Q sequences, for learning purposes.

It seemlessly parses FASTA/Q records from either plain files, (gzip, xz, bzip2, zstd, and lz4) compressed files or STDIN.

## Examples

```rust
use std::error::Error;

use fastseq::Reader;

fn main() -> Result<(), Box<dyn Error>> {
    // Read FASTA/Q records from stdin and write the number of records, total bases

    let in_file = "-";

    let mut reader = Reader::new(in_file)
        .map_err(|e| format!("failed to parse input file: {}: {}", in_file, e))?;

    let mut total_records: u64 = 0;
    let mut total_bases: u64 = 0;

    while let Some(res) = reader.next() {
        let seq = res?;

        total_records += 1;
        total_bases += seq.seq.len() as u64;
    }

    print!("{}\t{}\n", total_records, total_bases);

    Ok(())
}
```


## Benchmarking

see [fastx-tools](https://github.com/shenwei356/fastx-tools)


## Reference

- [needletail](https://github.com/onecodex/needletail)

## License

[MIT License](https://github.com/shenwei356/fastseq/blob/master/LICENSE)