Beefy Boxes and Bandwidth Generously Provided by pair Networks
Clear questions and runnable code
get the best and fastest answer
 
PerlMonks  

Re: translating multiple DNA sequence to protein sequence

by polypompholyx (Chaplain)
on Aug 19, 2013 at 11:16 UTC ( #1050014=note: print w/ replies, xml ) Need Help??


in reply to translating multiple DNA sequence to protein sequence

If you're planning on using Perl for bioinformatics, you might be better off installing BioPerl rather than hand-rolling FASTA parsers and translation codon tables.

use Bio::SeqIO; my $sequences = Bio::SeqIO->new( -file => "sequence.fasta", -format => "fasta", ); while ( my $dna = $sequences->next_seq ){ my $protein = $dna->translate( -codontable_id => 1, # standard genetic code -frame => 0, #reading-frame offset 0 ); print $dna->display_id, "\n"; print $protein->seq, "\n\n"; }
Having said that, installing BioPerl (1.6.901) on Windows seems to be more difficult than I was expecting: I had to resort to force with Strawberry and CPAN, having simply given up trying to get it to install with ActivePerl and PPM.


Comment on Re: translating multiple DNA sequence to protein sequence
Select or Download Code

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://1050014]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others cooling their heels in the Monastery: (8)
As of 2014-09-23 02:38 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    How do you remember the number of days in each month?











    Results (210 votes), past polls