擅长:python、mysql、java
<p>使用Python:</p>
<pre><code>strs = "atcttgttcaatggccgatcaacggtcgacaatcaa"
with open("1.fq") as f:
for line in f:
if line[20:24] == "aacg":
line = line.strip()
mismatches = sum(x!=y for x, y in zip(strs, line))
if mismatches <= 3:
print line, mismatches
atcttgttcaatggccgatcaacggtcgacaatcac 1
ttcttgttcaatggccgatcaacggtcgacaatcac 2
tccttgttcaatggccgatcaacggtcgacaatcac 3
</code></pre>