biors 0.46.0

Command-line tools for bio-rs biological AI model input workflows.
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

use serde_json::Value;
use std::process::Output;

mod common;

#[test]
fn fasta_validate_kind_flag_outputs_kind_aware_report() {
    let output = run_with_stdin(
        ["fasta", "validate", "--kind", "dna", "-"],
        ">seq1\nACGNU\n",
    );

    let value: Value = serde_json::from_slice(&output.stdout).expect("valid JSON output");
    assert_eq!(value["data"]["kind_counts"]["dna"], 1);
    assert_eq!(value["data"]["sequences"][0]["kind"], "dna");
    assert_eq!(value["data"]["sequences"][0]["alphabet"], "dna-iupac");
    assert_eq!(
        value["data"]["sequences"][0]["warnings"][0]["code"],
        "ambiguous_symbol"
    );
    assert_eq!(
        value["data"]["sequences"][0]["errors"][0]["code"],
        "invalid_symbol"
    );
}

#[test]
fn fasta_validate_auto_detects_mixed_sequence_kinds() {
    let output = run_with_stdin(
        ["fasta", "validate", "--kind", "auto", "-"],
        ">dna\nACGN\n>rna\nACGU\n>protein\nMEEPQSDPSV\n",
    );

    assert_mixed_kind_counts(&output);
}

#[test]
fn seq_validate_defaults_to_auto_detected_sequence_kinds() {
    let output = run_with_stdin(
        ["seq", "validate", "-"],
        ">dna\nACGN\n>rna\nACGU\n>protein\nMEEPQSDPSV\n",
    );

    assert_mixed_kind_counts(&output);
}

#[test]
fn seq_validate_kind_override_uses_kind_specific_messages() {
    let output = run_with_stdin(
        ["seq", "validate", "--kind", "rna", "-"],
        ">dna-looking\nACGT\n",
    );

    let value: Value = serde_json::from_slice(&output.stdout).expect("valid JSON output");
    assert_eq!(value["data"]["sequences"][0]["kind"], "rna");
    assert_eq!(value["data"]["sequences"][0]["errors"][0]["symbol"], "T");
    assert!(value["data"]["sequences"][0]["errors"][0]["message"]
        .as_str()
        .expect("message")
        .contains("RNA"));
}

fn run_with_stdin<const N: usize>(args: [&str; N], input: &str) -> Output {
    common::run_biors_stdin(&args, input)
}

fn assert_mixed_kind_counts(output: &Output) {
    let value: Value = serde_json::from_slice(&output.stdout).expect("valid JSON output");
    assert_eq!(value["data"]["kind_counts"]["dna"], 1);
    assert_eq!(value["data"]["kind_counts"]["rna"], 1);
    assert_eq!(value["data"]["kind_counts"]["protein"], 1);
    assert_eq!(value["data"]["sequences"][0]["kind"], "dna");
    assert_eq!(value["data"]["sequences"][1]["kind"], "rna");
    assert_eq!(value["data"]["sequences"][2]["kind"], "protein");
}