fxsplit 0.0.3

split FASTX into N chunks/files/headers
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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100

mod common;

use common::*;

#[test]
fn fasta_chunks_mean_records_per_file() {
    let t = tempdir();
    let input = join(t.path(), "input.fasta");
    let outdir = join(t.path(), "out_chunks");

    write_text(
        &input,
        ">r1\nAAAA\n>r2\nCCCC\n>r3\nGGGG\n>r4\nTTTT\n>r5\nNNNN\n",
    );

    run_split(vec![
        "--file".to_string(),
        path_str(&input),
        "--chunks".to_string(),
        "2".to_string(),
        "--threads".to_string(),
        "2".to_string(),
        "--outdir".to_string(),
        path_str(&outdir),
    ]);

    let files = list_files(&outdir);
    assert_eq!(files.len(), 3);

    let counts: Vec<usize> = files
        .iter()
        .map(|name| count_fasta_records(&join(&outdir, name), false))
        .collect();
    assert_eq!(counts, vec![2, 2, 1]);
}

#[test]
fn fasta_files_mean_number_of_output_files() {
    let t = tempdir();
    let input = join(t.path(), "input.fasta");
    let outdir = join(t.path(), "out_files");

    write_text(
        &input,
        ">r1\nAAAA\n>r2\nCCCC\n>r3\nGGGG\n>r4\nTTTT\n>r5\nNNNN\n",
    );

    run_split(vec![
        "--file".to_string(),
        path_str(&input),
        "--files".to_string(),
        "3".to_string(),
        "--threads".to_string(),
        "2".to_string(),
        "--outdir".to_string(),
        path_str(&outdir),
    ]);

    let files = list_files(&outdir);
    assert_eq!(files.len(), 3);

    let counts: Vec<usize> = files
        .iter()
        .map(|name| count_fasta_records(&join(&outdir, name), false))
        .collect();
    assert_eq!(counts, vec![2, 2, 1]);
}

#[test]
fn fasta_headers_use_sanitized_ids() {
    let t = tempdir();
    let input = join(t.path(), "input.fasta");
    let outdir = join(t.path(), "out_headers");

    write_text(&input, ">seq one\nAAAA\n>seq/two\nCCCC\n>seq/two\nGGGG\n");

    run_split(vec![
        "--file".to_string(),
        path_str(&input),
        "--headers".to_string(),
        "--threads".to_string(),
        "2".to_string(),
        "--outdir".to_string(),
        path_str(&outdir),
    ]);

    let files = list_files(&outdir);
    assert_eq!(
        files,
        vec![
            "seq.fasta".to_string(),
            "seq_two.fasta".to_string(),
            "seq_two_1.fasta".to_string()
        ]
    );

    for file in files {
        assert_eq!(count_fasta_records(&join(&outdir, &file), false), 1);
    }
}