[BioPython] problem accessing ncbi through	GenBank.NCBIDictionary
    Peter 
    biopython at maubp.freeserve.co.uk
       
    Tue Sep  4 13:05:21 UTC 2007
    
    
  
Miguel Ortiz-Lombardía wrote:
> Hello everyone.
> 
> I'm trying to retrieve from NCBI a series of GeneBank records from a list
> read from a file.
How many GenBenk identifiers are we talking about? Just trying to get an 
idea of the scale of the problem.  It certainly sounds like either 
network failures or timeouts.  Have you try something like this?
from Bio import GenBank
from urllib2 import HTTPError
ncbi_dict = GenBank.NCBIDictionary("protein", "genbank")
ids = ['14598510', '16904191']
output = open('saved.gb','w')
for gbid in ids:
     print "Fetching %s" % gbid
     try :
         gb_record = ncbi_dict[gbid]
     except HTTPError, e :
         #Check error code?
         print str(e)
         print "Re-trying %s" % gbid
         gb_record = ncbi_dict[gbid]
     output.write(gb_record)
output.close()
print "Done"
Peter
    
    
More information about the Biopython
mailing list