predictosaurus 0.3.0

Uncertainty aware haplotype based genomic variant effect prediction
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

use crate::transcription;
use crate::translation::amino_acids::AminoAcid;
use anyhow::Result;
use itertools::Itertools;

pub(crate) mod amino_acids;

/// Translates a DNA sequence to a protein sequence
pub(crate) fn dna_to_amino_acids(dna: &[u8]) -> Result<Vec<AminoAcid>> {
    let rna = transcription::transcribe_dna_to_rna(dna)?;
    let amino_acids = rna
        .iter()
        .tuples::<(_, _, _)>()
        .map(|codon| AminoAcid::from_codon(&[*codon.0, *codon.1, *codon.2]))
        .take_while(|result| result.is_ok())
        .collect::<Result<Vec<AminoAcid>>>()?;
    Ok(amino_acids)
}

#[cfg(test)]
mod tests {
    use super::*;
    #[test]
    fn test_dna_to_protein() -> Result<()> {
        assert_eq!(
            dna_to_amino_acids(b"ATGCGCGTA")?,
            vec![
                AminoAcid::Methionine,
                AminoAcid::Arginine,
                AminoAcid::Valine,
            ]
        );
        assert_eq!(
            dna_to_amino_acids(b"ATGCGCGTAAATGCGCGT")?,
            vec![
                AminoAcid::Methionine,
                AminoAcid::Arginine,
                AminoAcid::Valine,
                AminoAcid::Asparagine,
                AminoAcid::Alanine,
                AminoAcid::Arginine,
            ]
        );
        Ok(())
    }
}