[Bioperl-l] Bio::DB::GenBank question (acc vs. version)

George Hartzell hartzell at alerce.com
Sun Sep 13 01:35:44 UTC 2009


It looks like get Bio::DB::GenBank::get_Seq_by_{version,acc} are
functionally identical.  They seem to trickle down to the same place
and walking through these two requests yields almost identical http
requests: 

  $db->get_Seq_by_version('J00522.1')
  GET http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?retmode=text&rettype=gbwithparts&db=nucleotide&tool=bioperl&id=J00522.1&usehistory=n

  $db->get_Seq_by_acc('J00522')
  GET http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?retmode=text&rettype=gbwithparts&db=nucleotide&tool=bioperl&id=J00522&usehistory=n

The only difference that I can see is that they index into different
secions of %PARAMSTRING defined in Bio::DB::GenBank, but those
sections contain the same information.

I'd like a general purpose tool that does The Right Thing whether
there's a .1 on the end of an identifier or not, and am just trying to
make sure I'm not doing something troublesome.

Am I correct about the above?

While I'm at it, I think that the comment

  # note that get_Stream_by_version is not implemented

in Bio::DB::GenBank was made obsolete by whoever commented out the

  $self->throw(...)

in get_Stream_by_version in Bio::WebDBSeqI.pm.

I'll happily commit the trivial doc fix if no one shoots down the
idea. (can't help big, might as well help small...).

Thanks,

g.



More information about the Bioperl-l mailing list