Regarding writing readable Perl code, there are some automatic plug-ins for VIM, EMACs or JEDIT, or you can simply use PerlTidy: http://perltidy.sourceforge.net/
Bioperl: Finding Restriction Enzyme Sites
#!usr/bin/perl -w
use Bio::DB::GenBank;
use Bio::Restriction::EnzymeCollection;
use Bio::Restriction::Analysis;
my $gb=new Bio::DB::GenBank();
my $seq=$gb->get_Seq_by_acc('AF303112');
print $seq->seq,"\n";
my $all_collection=Bio::Restriction::EnzymeCollection;
my $my_enzyme=$all_collection->get_enzyme('EcoRI');
my $analysis=Bio::Restriction::Analysis->new(-seq=>$seq);
@fragment=$analysis->fragments($my_enzyme);
print "@fragment\n";
my $pro_obj=$seq->translate(-complete=>1);
print $pro_obj->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!
• 4,704 views
•
link
updated
by
Neilfws
•
written
by
cindylanzao •
•
written
by
cindylanzao •
1 answer
First: if you're going to declare your variables using my, you should:
use strict;
Second - now that we are using strict, you need to change @fragment to:
my @fragment = $analysis->fragments($my_enzyme);
Third: having done those things, you should now see the error:
Bareword "Bio::Restriction::EnzymeCollection" not allowed while "strict subs" in use at...
I suspect that what you want is:
my $all_collection = Bio::Restriction::EnzymeCollection->new;
Your code should now run and print the sequence, fragment and translation. In this example, there does not appear to be an EcoRI site since the sequence and fragment are the same.
Finally: Perl code is much more readable if you put in more spaces and align related lines. For example, compare this to what you typed:
my $all_collection = Bio::Restriction::EnzymeCollection->new;
my $my_enzyme = $all_collection->get_enzyme('EcoRI');
my $analysis = Bio::Restriction::Analysis->new(-seq => $seq);
my @fragment = $analysis->fragments($my_enzyme);
Log in to answer this question.
More posts like this
-
Bioperl - how can i print first result of result sequence per iteration?
written by Xylanaser
-
Perl script to retrieve species name from Swissprot
written by anicet.ebou
-
Get Fasta file with Protein Sequences given a file with Genbank Ids using Perl
written by Andrea •
-
Extracting Random Sequences From A Large Fasta File Using Bioperl
written by joshualevipayne
-
Getting Only Orthologs From Ensembl Database
written by jaip217 •
-
Parsing A Genbank Record With Bioperl
written by robjohn7000
-
Bioperl Retrieve Fasta Sequence With Colon In Header
written by jason.r.gallant •
-
Bio::Db::Genbank With Coordinates Reports A Strange Error
written by Puriney •
-
Bioperl Sirna Design
written by User 3992 •
-
Getting Genome Coordinates From Refseq Exon Mrna Position Data?
written by Krisr
(1) I've indented your lines of code with 4 spaces: note how this improves readability. (2) I think you wanted to post this under category "question", not "forum"; changed that too.
Thank you very much,
Changed the title too. "Bioperl" is not very informative. Nor were "beginner's" or "task" as tags. Good titles and tags make your question easier to find and answer.