[Biopython] Searching for and downloading sequences using the history

Carlos Javier Borroto carlos.borroto at gmail.com
Fri Sep 18 12:59:58 EDT 2009


Hi all,

I'm trying to download all of the EST from a specie, I'm following the
example on the tutorial which seems to be exactly what I need. But I
running into this problem:

>>> from Bio import Entrez
>>> Entrez.email = "carlos.borroto at gmail.com"
>>> dbname = "nucest"
>>> query_term = "Genus specie"
>>> search_handle = Entrez.esearch(db=dbname,term=query_term,usehistory="y")
>>> search_results = Entrez.read(search_handle)
>>> search_handle.close()
>>> len(search_results["IdList"])
20
>>> print search_results["Count"]
193951

So the assert statement if failing:
>>> gi_list = search_results["IdList"]
>>> count = int(search_results["Count"])
>>> assert count == len(gi_list)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
AssertionError

And most important I'm not getting all of the ids.

Did someone knows what I'm doing wrong?

thanks in advance
-- 
Carlos Javier Borroto
Baltimore, MD
Phone: (410) 929 4020


More information about the Biopython mailing list