[Biopython-dev] [Bug 2381] back-translate

[Peter]

Howard Salis wrote:
> 
> What does the current back-translate function output?
> 

Short example,

 >>> from Bio import Translate
 >>> from Bio.Seq import Seq
 >>> from Bio.Alphabet.IUPAC import unambiguous_dna
 >>> my_dna = Seq("GCCGCATGCATAGATAGATAG", unambiguous_dna)
 >>> my_prot = Translate.unambiguous_dna_by_id[11].translate(my_dna)
 >>> my_prot
Seq('AACIDR*', HasStopCodon(IUPACProtein(), '*'))
 >>> Translate.unambiguous_dna_by_id[11].back_translate(my_prot)
Seq('GCTGCTTGTATTGATCGTTAA', IUPACUnambiguousDNA())
 >>> my_dna
Seq('GCCGCATGCATAGATAGATAG', IUPACUnambiguousDNA())

i.e. The current back_translate picks one possible back translation.

Peter