1 Jun 10:06
Extract genomic data from GenBank
Will Bryant <w.bryant <at> ucl.ac.uk>
2009-06-01 08:06:58 GMT
2009-06-01 08:06:58 GMT
I'm trying to retrieve the complete GenBank format sequence file for a
specified bacterium using get_Seq_by_gi, but I keep getting 'gi does not
exist' errors, even when trying the example gi '405830'. The script was
running fine September last year, but when I came back to it this week
it wasn't working. Am I missing something obvious?
In case it's important, I'm using ActivePerl 5.10.0, bioperl 1.5.2_100
Code:
#!/usr/bin/perl -w
use strict;
use Bio::Perl;
use Bio::DB::GenBank;
my $gb = new Bio::DB::GenBank(-db => 'genome', -format => 'genbank');
my $straincomp = $gb->get_Seq_by_gi('405830');
my $seqout = 0;
#my $set_output_file = '$seqout = Bio::SeqIO->new( -format =>
\'genbank\', -file =>
\'>c:\\phd\\modelling\\working\\gi'.$ARGV[0].'_data.gb\');';
#print $set_output_file;
eval ($set_output_file);
$seqout -> write_seq($straincomp);
(Continue reading)
RSS Feed