[Bioperl-l] Frame translation still gets an extra aa

Christine Hoogland christine.hoogland at gmail.com
Thu Jan 30 02:08:39 UTC 2014


Hi,

Back in 2011 there was a question/comment about frame translation behaviour:
http://bioperl.org/pipermail/bioperl-l/2011-January/034422.html

I have perl v5.14.2 and BioPerl v.1.6.910, and still get this extra AA 
using or not the complete_codons option:

$ perl -MBio::Seq -e 'print "$Bio::Seq::VERSION\n";';
1.005002102

$ perl -l -MBio::Seq -e 
'$x=Bio::Seq->new(-display_id=>"foo",-seq=>"AAACCCTTTGGG"); print 
$x->translate(-frame=>1, -complete_codons => 1)->seq'
NPLG

$ perl -l -MBio::Seq -e 
'$x=Bio::Seq->new(-display_id=>"foo",-seq=>"AAACCCTTTGGG"); print 
$x->translate(-frame=>1, -complete_codons => 0)->seq'
NPLG

I would expect NPL only, when asking for codon not completed 
(complete_codons 0, supposed to be the default):
a-aac-cct-ttg-gg
   N   P   L

What am I doing wrong?
Thanks,
Christine



More information about the Bioperl-l mailing list