擅长:python、mysql、java
<p>有没有理由不使用生物圈?在</p>
<pre><code>from Bio import AlignIO
alignment =AlignIO.read("alignment.fas", "fasta")
n=0
while n<len(alignment[0]):
A=alignment[:,n].count('A')
C=alignment[:,n].count('C')
G=alignment[:,n].count('G')
T=alignment[:,n].count('T')
gap=alignment[:,n].count('-')
print "at position %s there are %s A's, %s C's, %s G's, %s T's and %s gaps" % (n, A, C, G, T, gap)
n=n+1
</code></pre>
<p>确保您有一个真正的对齐方式(即序列长度相同)。<br/>
p、 我为印刷声明的丑陋格式道歉。。。在</p>
<p>这就回来了</p>
^{pr2}$