Bioperl beginner's program

[cindylanzao]

!usr/bin/perl -w

use Bio:B::GenBank;
use Bio::Restriction::EnzymeCollection;
use Bio::Restriction::Analysis;
my $gb=new Bio:B::GenBank();
my $seq=$gb->getSeqbyacc('AF303112');
print $seq->seq,"\n";
my $allcollection=Bio::Restriction::EnzymeCollection;
my $myenzyme=$allcollection->getenzyme('EcoRI');
my $analysis=Bio::Restriction::Analysis->new(-seq=>$seq); @fragment=$analysis->fragments($myenzyme);
print "@fragment\n";
my $proobj=$seq->translate(-complete=>1);
print $proobj->seq;

I want to get a sequence by accession number, and then find if there is a EcoRI enzyme site in this sequence and finally translate this sequence into protein. But the output is only the Nuclear Acid sequence, what's wrong, can someone tell me? Thanks!

[krobison]

Note: if you use the advanced editor, you can wrap your code in tags which will make it much easier to read (I've added use strict, which you should use strictly!). My tweaks to your code are in bold.

!usr/bin/perl -w
use strict;
use Bio:B::GenBank; 
use Bio::Restriction::EnzymeCollection; 
use Bio::Restriction::Analysis; 
my $gb=new Bio:B::GenBank(); 
my $seq=$gb->getSeqbyacc('AF303112'); 
print $seq->seq,"\n"; 
my $allcollection=Bio::Restriction::EnzymeCollection; 
my $myenzyme=$allcollection->getenzyme('EcoRI'); 
my $analysis=Bio::Restriction::Analysis->new(-seq=>$seq); 
my @fragment=$analysis->fragments($myenzyme); 
print "@fragment\n"; 
my $proobj=$seq->translate(-complete=>1); 
print $proobj->seq;

My next suggestion is that your learn the Perl debugger. If you use EMACS, you can run the debugger very nicely in the editor & easily set breakpoints.

I'm afraid I don't see an obvious problem, but if you look at what is contained in @fragment it might give a hint. Are you sure this sequence has an EcoRI site?