0: #!/usr/bin/perl -w
1:
2: use strict;
3: use Bio::DB::GenBank;
4: use Bio::SeqIO;
5: use File::Basename;
6:
7: # configuration information
8: my $outfile = 'outfile';
9: my $progname = basename $0;
10:
11: # usage information
12: my $usage = "usage: $0 accession1 [accession2 ...]\n";
13: die $usage unless @ARGV;
14:
15: # start program flow
16: my $gb = new Bio::DB::GenBank(
17: -retrievaltype => 'tempfile',
18: -format => 'Fasta'
19: );
20:
21: my $seqout = Bio::SeqIO->new( -file => ">$outfile", -format => 'FASTA' );
22:
23: my $seqio = $gb->get_Stream_by_acc( @ARGV );
24: while( my $sequence = $seqio->next_seq ) {
25: $seqout->write_seq($sequence);
26: } In reply to Fetching Biological Sequences from the NCBI by biosysadmin
| For: | Use: | ||
| & | & | ||
| < | < | ||
| > | > | ||
| [ | [ | ||
| ] | ] |