[Bioperl-l] How do you change the default on Bio::DB::GenBank?
    Mark Wilkinson 
    markw at illuminae.com
       
    Wed Jul  2 23:31:54 EDT 2003
    
    
  
Hi BioPerlers, 
Can someone tell me how to change the default db name in
Bio::DB::GenBank?  I've spent the past two hours in debug mode and can't
see any point at which an argument for the database name is even tested
- it is coded to default to 'protein' which isn't very useful to me. 
Such a switch doesn't appear anywhere in the documentation either as far
as I can tell.  I took a guess and thought the following would work:
use Bio::DB::GenBank; 
my $gb = Bio::DB::GenBank->new(-db=>"nucleotide");
print $gb->get_Seq_by_acc("AF005495");
But that still sends out a GET request to the protein database:
http://www.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?retmode=text&rettype=gb&db=protein&tool=bioperl&id=J00522&usehistory=n
Whereas what I want it to do is:
http://www.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?retmode=text&rettype=gb&db=nucleotide&tool=bioperl&id=J00522&usehistory=n
Help please :-)
Cheers all!
M
-- 
Mark Wilkinson <markw at illuminae.com>
Illuminae
    
    
More information about the Bioperl-l
mailing list