hyper-gen 0.2.2

HyperGen is a high-performance Rust library to sketch genomics files into hypervectors and realize fast Average Nucleotide Identity (ANI) approximation.
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

use std::fs::File;
use std::io::{BufRead, BufReader};
use std::path::PathBuf;

// Read merged sequences from a genome file into single u8 vector
pub fn read_merge_seq(file_name: &PathBuf) -> Vec<u8> {
    let mut fna_seqs = Vec::<u8>::new();

    let file = File::open(file_name).unwrap();
    let mut reader = BufReader::new(file);

    let mut buf = String::new();
    while reader.read_line(&mut buf).unwrap() > 0 {
        if !buf.starts_with('>') {
            if buf.ends_with('\n') {
                buf.pop(); // remove \n character
            }
            if buf.ends_with('\r') {
                buf.pop();
            }
            fna_seqs.extend_from_slice(buf.as_bytes());
        } else {
            fna_seqs.push(b'N');
        }
        buf.clear();
    }

    fna_seqs
}