orf 0.1.0

Translates nucleotide sequences to Amino Acid sequences. (RNA and DNA, shifts reading frames and reverses)
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

#![feature(question_mark)]
extern crate orf;
extern crate memmap;
extern crate num_cpus;

use std::fs::File;
use std::process::exit;
use memmap::{Mmap, Protection};

fn main() {
    if let Err(e) = transcribe() {
        println!("Error: {}", e);
        exit(1);
    }
}

fn transcribe() -> orf::Result {
    let mmap = Mmap::open_path("tests/test-nucleo.FASTA", Protection::Read)?;
    let input = unsafe { mmap.as_slice() };
    let mut output_file = File::create("results.txt")?;
    let n_threads = num_cpus::get() as u32;

    orf::start_parse(input, &mut output_file, n_threads)?;

    output_file.sync_all()?;

    Ok(())
}