Search
Question: Biostrings - translate gives wrong results for triplets TTG and CTG
0
gravatar for janas.sebastien
11 weeks ago by
janas.sebastien0 wrote:

All,

I also posted that message on Github, since I don't know what's the best place to report potential bugs.

Since at least Biostring 2.42, the function "translate" gives wrong results for triplets TTG and CTG.

translate(DNAString("TTG"))
1-letter "AAString" instance
seq: M

translate(DNAString("CTG"))
1-letter "AAString" instance
seq: M

=> Both should be L (leucine).

We think it's related to the erroneous default use of alternative genetic codes.

Cheers,

Seb

ADD COMMENTlink written 11 weeks ago by janas.sebastien0

I now understand that it's considered as an alternative start codon, that's why it's translated in M. My issue was that the change in the behavior of the code suddenly broke our unit tests. It would have be better to implement the new behavior as an option to translate. Because now, all people who are using translate to translate a single triplet will see their code broken.

ADD REPLYlink modified 11 weeks ago • written 11 weeks ago by janas.sebastien0

Moreover, the documentation does not say that the sequence must be a complete sequence (starting at the start codon).

ADD REPLYlink written 11 weeks ago by janas.sebastien0

It's fine to track this on github only.

ADD REPLYlink written 11 weeks ago by Martin Morgan ♦♦ 22k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 131 users visited in the last hour