fastq_rs 0.0.2

Multi purpose fastq toolkit.
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

use crate::common::{general_bufwriter, needletail_fastq_reader};
use anyhow::Result;
use regex::Regex;
use std::path::PathBuf;

pub fn fastq_grep(fastq: Option<PathBuf>, pattern: String, outfile: Option<PathBuf>) -> Result<()> {
    let mut reader = needletail_fastq_reader(fastq)?;
    let mut writer = general_bufwriter(outfile)?;

    let pattern = Regex::new(pattern.as_str())?;

    while let Some(record) = reader.next() {
        let record = match record {
            Ok(record) => record,
            Err(_) => continue,
        };

        match pattern.captures(std::str::from_utf8(record.id())?) {
            Some(_) => record.write(&mut writer, None)?,
            None => continue,
        }
    }

    Ok(())
}