擅长:python、mysql、java
<p>读过你重复的问题<a href="http://www.biostars.org/p/94164/" rel="nofollow">http://www.biostars.org/p/94164/</a>(请不要像这样重复发帖),下面是一个最小的Biopython答案:</p>
<pre><code>import sys
from Bio import SeqIO
filename = sys.argv[1] # Takes first command line argument input filename
for record in SeqIO.parse(filename, "genbank"):
for feature in record.features:
if feature.type == "CDS":
locus_tag = feature.qualifiers.get("locus_tag", ["???"])[0]
product = feature.qualifiers.get("product", ["???"])[0]
print("%s\t%s" % (locus_tag, product))
</code></pre>
<p>只需稍作修改,您就可以将其写入文件。在</p>