1

I am writing some code intended to translate ambiguous DNA codes into possible amino acids and I am seeing some strange translation from the Biopython 1.56 package. It appears to be translating ambiguous DNA codes to 'J' which does not exist as a code for anything. I am running python 2.6.1 on Mac OS 10.6.6.

For example:

>>>from Bio.Seq import *
>>>translate('ARAWTAGKAMTA')
'XJXJ'

or

>>>from Bio.Seq import Seq
>>>c = Seq('ARAWTAGKAMTA')
>>>c.translate().tostring()
'XJXJ'

I have looked through the Bio.Data.CodonTable source and Bio.Seq source and I cannot find a reason why this would be happening. Any ideas?

Thanks!

Mark

2
  • The expected outcome should be 'XXXX'. I think the best answer came to my parallel post on biostar (biostar.stackexchange.com/questions/5688/…). Although it still does not natively solve the issue, it is a work around I can live with for now. Feb 24, 2011 at 20:50
  • +1 for the question, -1 because you could not find J is a well known letter to cope with L/I ambiguity in mass based sequencing methods. This equals zero.
    – joaquin
    Dec 3, 2011 at 7:18

1 Answer 1

4

J is Leucine (L) or Isoleucine (I), used in mass-spec (NMR).

See also http://biostar.stackexchange.com/questions/5688/biopython-translation-error

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.