[Bioperl-l] clone ID to Accession ID genbank
Iris Buiton
irisbuiton at gmail.com
Thu Feb 4 16:06:25 UTC 2010
Hi members of Bioperl,
I am trying to run some code that given a clone ID (for instance,
CH255-149A8) get the accesion ID (AC238489.2).
The problem is that using this ID, it works perfectly, but using this
other ID (CH251-558O8) I do not get anything (not $seq object is
obtained). Checking to Genbank, this clone exist, there is seqeunce,
and the accession ID should be AC147515.1
This is the code:
my $name_file="CH251-558O8";
my $query = Bio::DB::Query::GenBank->new
(-query => $name_file,
-reldate => '30',
-db => 'nucleotide');
#eval 'my $seqio = $gb->get_Stream_by_query($query)';
print $query->count,"\n";
my $gb = Bio::DB::GenBank->new();
my $seqio = $gb->get_Stream_by_query($query);
print "Running $name_file.....\n";
while( my $seq = $seqio->next_seq ) {
if ($seq->accession_number)
{
print "$name_file = ", $seq->accession_number," = ", $seq-
>length,"\n" ;
}else
{
print "$name_file = #N/A \n" ;
}
}
Any ideas what am I doing wrong?
Thanks!
Iris
More information about the Bioperl-l
mailing list