python练习
2020-05-26 本文已影响0人
M78_a
image.png
import sys
from Bio import SeqIO
with open(sys.argv[1],'r') as f:
for s in SeqIO.parse(f, "fasta"):
#print(s.seq)
if s.seq[0:3] == "ATG":
print(s.id)
