!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!
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!
Comment