擅长:python、mysql、java
<p>根据要求,Biopython解决方案:</p>
<pre><code>from Bio import Align
p = "ATCAT"
s = "ATGCAT"
s_annot = "135198"
aligner = Align.PairwiseAligner()
alignment = str(aligner.align(p, s)[0]).split()
middle = alignment.pop(1)
alignment.append("".join(c if m == "|" else m for c, m in zip(s_annot, middle)))
print("\n".join(alignment))
</code></pre>
<p>输出:</p>
<pre><code>AT-CAT
ATGCAT
13-198
</code></pre>