擅长:python、mysql、java
<p>我的解决方案是:</p>
<p>将<code>TC</code>替换为<code>T C</code>,将<code>GC</code>替换为{<cd4>}(这是根据给定的索引完成的),然后根据空格字符拆分。。。。在</p>
<pre><code>def digest(seqs, enzymes):
res = []
for li in seqs:
for en in enzymes:
li = li.replace(en[0],en[0][:en[1]]+" " + en[0][en[1]:])
r = li.split()
res.append(r)
return res
seqs = ["AATTCCGGTCGGGGCTCGGGGG","AAAGCAAAATCAAAAAAGCAAAAAATC"]
enzymes = [["TC", 1],["GC",1]]
#enzymes = [["AAT", 2],["GC",1]]
print seqs
print digest(seqs, enzymes)
</code></pre>
<p><strong>结果是:</strong></p>
<p>对于<code>([["TC", 1],["GC",1]])</code></p>
^{pr2}$
<p>对于<code>([["AAT", 2],["GC",1]])</code></p>
<pre><code>['AATTCCGGTCGGGGCTCGGGGG', 'AAAGCAAAATCAAAAAAGCAAAAAATC']
[['AA', 'TTCCGGTCGGGG', 'CTCGGGGG'], ['AAAG', 'CAAAA', 'TCAAAAAAG', 'CAAAAAA', '
TC']]
</code></pre>