将e-utils命令转换为Bio.Entrez(BioPython)中的等效命令

时间:2019-03-19 12:21:53

标签: python-3.x biopython

使用BioPython连接每个命令时遇到麻烦。

有人可以帮助我使用BioPython将此命令行转换为等效命令行吗?

esearch -db assembly -query "GCF_002514765.1" | elink -target biosample | efetch -format xml

谢谢

1 个答案:

答案 0 :(得分:0)

摘要可能会解决。

from Bio import Entrez

acc = "GCF_002514765.1"
search = Entrez.read(Entrez.esearch(db="assembly", term=acc))
summary = Entrez.read(Entrez.esummary(db="assembly",id=idsearch['IdList'],report="full"))