nucgen 0.2.0

A simple tool and library for generating random nucleotide 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
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74

mod cli;

use cli::Args;

use anyhow::Result;
use clap::Parser;
use nucgen::{write_fasta, write_fastq, Format, Sequence};
use std::io::Write;

fn gen_single_end(args: Args) -> Result<()> {
    let mut rng = args.get_rng();
    let mut seq = Sequence::new();
    let qual = vec![b'?'; args.slen];

    let mut out = args.output_handle_single()?;
    let format = args.format();

    // Write the records
    for idx in 0..args.num_records {
        seq.fill_buffer(&mut rng, args.slen);
        match format {
            Format::Fasta => write_fasta(&mut out, idx, seq.bytes())?,
            Format::Fastq => write_fastq(&mut out, idx, seq.bytes(), &qual)?,
        }
    }

    // Flush the output buffer
    out.flush()?;

    Ok(())
}

fn gen_paired_end(args: Args) -> Result<()> {
    let mut rng = args.get_rng();
    let mut s1 = Sequence::new();
    let mut s2 = Sequence::new();
    let q1 = vec![b'?'; args.slen];
    let q2 = vec![b'?'; args.xlen];

    let (mut out_r1, mut out_r2) = args.output_handle_paired()?;
    let format = args.format();

    // Write the records
    for idx in 0..args.num_records {
        s1.fill_buffer(&mut rng, args.slen);
        s2.fill_buffer(&mut rng, args.xlen);
        match format {
            Format::Fasta => {
                write_fasta(&mut out_r1, idx, s1.bytes())?;
                write_fasta(&mut out_r2, idx, s2.bytes())?;
            }
            Format::Fastq => {
                write_fastq(&mut out_r1, idx, s1.bytes(), &q1)?;
                write_fastq(&mut out_r2, idx, s2.bytes(), &q2)?;
            }
        }
    }

    // Flush the output buffers
    out_r1.flush()?;
    out_r2.flush()?;

    Ok(())
}

fn main() -> Result<()> {
    let args = Args::parse();

    if args.paired() {
        gen_paired_end(args)
    } else {
        gen_single_end(args)
    }
}