擅长:python、mysql、java
<p>这将同时迭代每个索引,并为每个核苷酸出现添加一个索引</p>
<pre><code>result = {'A': 13*[0], 'G': 13*[0], 'T': 13*[0], 'C': 13*[0]}
for index, sequence in enumerate(zip(*sequences)):
for nucleotide in sequence:
result[nucleotide][index] += 1
</code></pre>
<p>输出:</p>
<pre><code>{'A': [1, 4, 1, 0, 0, 3, 4, 1, 1, 3, 0, 2, 0], 'G': [1, 0, 4, 6, 0, 0, 1, 3, 1, 0, 0, 1, 2], 'T': [3, 0, 0, 0, 6, 1, 0, 2, 3, 3, 2, 3, 0], 'C': [1, 2, 1, 0, 0, 2, 1, 0, 1, 0, 4, 0, 4]}
</code></pre>