我想在pubmed数据库中搜索一个基因符号列表(下面命名为t),以便(最终)检索相关基因的DNA序列。我想限制我的搜索只限于人类,但我现在的代码给了我除了人类以外的有机体。
from Bio import Entrez
Entrez.email = '...' #my email: always tell Entrez who you are
t = ['FOXO3']
for i in range(len(t)):
search = 'human[orgn]'+t[i]
handle = Entrez.esearch(db='gene',term=search)
record = Entrez.read(handle)
t = record[u'IdList']
handle = Entrez.efetch('nucleotide',id=t[0],rettype='gb',retmode='text')
print handle.read()
有人可以看到我出错的地方吗?
答案 0 :(得分:0)
你弄乱了数据库。在您使用db=gene
的搜索中,但在efetch中,您将其更改为db=nucleotide
。它们是不同的东西: