擅长:python、mysql、java
<p>首先是基因名,例如:ATK1</p>
<pre><code>item = 'ATK1'
animal = 'Homo sapien'
search_string = item+"[Gene] AND "+animal+"[Organism] AND mRNA[Filter] AND RefSeq[Filter]"
</code></pre>
<p>现在我们有一个搜索字符串来搜索id</p>
^{pr2}$
<p>如果没有找到id,则返回id作为一个列表。现在假设它返回列表中的1项。在</p>
<pre><code>seq_id = ids[0] #you must implement an if to deal with <0 or >1 cases
handle = Entrez.efetch(db="nucleotide", id=seq_id, rettype="fasta", retmode="text")
record = handleA.read()
</code></pre>
<p>这将为您提供一个fasta字符串,您可以将其保存到文件中</p>
<pre><code>out_handle = open('myfasta.fasta', 'w')
out_handle.write(record.rstrip('\n'))
</code></pre>