[Bioperl-l] Frame translation still gets an extra aa
Christine Hoogland
christine.hoogland at gmail.com
Thu Jan 30 03:18:16 UTC 2014
Hi Chris,
Thanks a lot for your fast answer. I updated BioPerl with
cpan -i 'Bio::Seq'
and I have now 1.006922
and the translation is working as expected!
with complete_codons: translate to NPLG
without (default): translate to NPL
Previously I installed following instructions at
http://www.bioperl.org/wiki/Installing_BioPerl_on_Ubuntu_Server
and didn't realised it was referring to current_core_unstable.tar.bz2...
probably not the latest version!!
Thanks again!
Christine
On 30/01/14 12:42, Fields, Christopher J wrote:
> The Bioperl version is wrong if you have v1.6.910; it should be something like 1.006910. Do you have an older version of bioperl somewhere in your perl library path? Try:
>
> perldoc -l Bio::Seq
>
> chris
>
> On Jan 29, 2014, at 8:08 PM, Christine Hoogland <christine.hoogland at gmail.com> wrote:
>
>> 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
>> _______________________________________________
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
--
Christine Hoogland
64/316 Long Street East
Graceville QLD 4075
Australia
Tel: +617 3129 9973
Mobile: +614 7891 9833
Email: christine.hoogland at gmail.com
LinkedIn: http://ch.linkedin.com/pub/christine-hoogland/23/778/333
Viadeo: http://www.viadeo.com/fr/profile/christine.hoogland
More information about the Bioperl-l
mailing list