[Biopython] Entrez.efetch from gene

Michael Fahy fahy at chapman.edu
Tue Jan 25 20:34:50 EST 2011


Trying to use Entrez.efetch() to query the gene database.

The efetch help at
http://www.ncbi.nlm.nih.gov/entrez/query/static/efetchseq_help.html says
there are no retrieval types supported by the gene database.  If I do an
efetch query without specifying a value for rettype, it returns html.  Is
there a way in Biopython to parse this html?  Or is there another way to
query the gene database so it will return data that can be parsed?

Sample code:

from Bio import Entrez
Entrez.email = 'email at chapman.edu'

search_database = 'gene'
search_term = 'YIL065C'

handle = Entrez.esearch(db='gene',term=search_term)
record = Entrez.read(handle)
reclist = record['IdList']  

handle = Entrez.efetch(db=search_database,  id =reclist[0])

myrecord = handle.read()
print myrecord

----------------------------------------------------------
Michael A. Fahy
fahy at chapman.edu




More information about the Biopython mailing list