|
|
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ 感谢参与,应助指数 +1 wizardfan: 金币+10, BioEPI+1, 感谢热心的回答 2013-04-24 06:31:58
GO to "http://www.ncbi.nlm.nih.gov/genomes/OrganelleResource.cgi?opt=organelle&taxid=6544"
click "Download"-> "download the list of accessions"-> fileA
GO to "http://www.ncbi.nlm.nih.gov/sites/batchentrez"
upload the fileA and got the all GBk webpage
click "send to" to save file with "Genebank full " -> file B
install perl and bioperl package
run "perl GetCds.pl fileB FileYouWanted"
GetCds.pl :
#!/usr/bin/perl
use Bio::SeqFeatureI;
use Bio: erl;
use Bio::Seq;
my %species;
@seq=read_all_sequences("$ARGV[0]",'genbank');
$flag=1;
foreach $seqobj (@seq) {
@features=$seqobj->get_SeqFeatures();
foreach $feat (@features) {
if($feat->primary_tag eq 'CDS'){
print "$feat\n";
@name=$feat->get_tag_values('gene');
$name=join "",@name;
open OUT,">>$ARGV[1]";
print OUT '>'.$seqobj->desc().$name."\n";
print OUT $feat->seq->seq()."\n";
close OUT;
}
}
} |
|