seq_rev_translate {bioseq}

Reverse translate amino acid sequences

Description

The function perform reverse translation of amino acid sequences. Such operation does not exist in nature but is provided for completeness. Because of codon degeneracy it is expected to produce many ambiguous nucleotides.

Usage

seq_rev_translate(x, code = 1)

Arguments

x

an amino acid sequence (bioseq_aa)

code

an integer indicating the genetic code to use for reverse translation (default 1 uses the Standard genetic code). See Details.

Details

Gaps (-) are interpreted as unknown amino acids (X) but can be removed prior to the translation with the function seq_remove_gap.

Value

a vector of DNA sequences.

See Also

Other biological operations: rev_complement, seq_translate(), transcription

Examples

x <- dna("ACTTTGGCTAAG")
y <- seq_translate(x)
z <- seq_rev_translate(y)
z
# There is a loss of information during the reverse translation
all.equal(x, z)