擅长:python、mysql、java
<p>将<a href="https://biopython.org/" rel="noreferrer">Biopython</a>用于此目的将更加健壮和简单:</p>
<pre><code>from Bio import SeqIO
motif = 'GGAC'
for record in SeqIO.parse('input.fa', 'fasta'):
if motif in record.seq:
print(record.id)
</code></pre>