mmft 0.2.1

A minimal fasta 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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47

use crate::utils::{error, stdin};
use anyhow::{bail, Result};
use bio::io::fasta;
use std::io;

pub fn get_number_seq_bases(matches: &clap::ArgMatches) -> Result<()> {
    let input_file = crate::get_fasta_files(matches);

    match input_file {
        // read directly from files
        Some(f) => {
            for el in f.iter() {
                let basename = crate::get_basename_from_pathbuf(el)?;

                let mut reader = fasta::Reader::from_file(el)?.records();
                let mut nb_reads = 0;
                let mut nb_bases = 0;

                while let Some(Ok(record)) = reader.next() {
                    nb_reads += 1;
                    nb_bases += record.seq().len();
                }
                // to stdout
                println!("{}\t{}\t{}", basename, nb_reads, nb_bases)
            }
        }
        // read from stdin
        None => match stdin::is_stdin() {
            true => {
                let mut records = fasta::Reader::new(io::stdin()).records();
                let mut nb_reads = 0;
                let mut nb_bases = 0;

                while let Some(Ok(record)) = records.next() {
                    nb_reads += 1;
                    nb_bases += record.seq().len();
                }
                // to stdout
                println!("{}\t{}", nb_reads, nb_bases)
            }
            false => {
                bail!(error::StdinError::NoSequence)
            }
        },
    }
    Ok(())
}