擅长:python、mysql、java
<p>正如<strong>Katrielex</strong><a href="https://stackoverflow.com/questions/14397504/how-to-count-amino-acids-in-fasta-formated-file/14397815#comment20031034_14397504">points out</a>,<a href="http://docs.python.org/2/library/collections.html#collections.Counter" rel="nofollow noreferrer">^{<cd1>}</a>非常适合这项任务:</p>
<pre><code>In [1]: from collections import Counter
In [2]: Counter('MRMRGRRLLPIIL')
Out[2]: Counter({'R': 4, 'L': 3, 'M': 2, 'I': 2, 'G': 1, 'P': 1})
</code></pre>
<p>您可以将其应用于代码中<code>sequences</code>dict的值。在</p>
<p>但是,我建议不要在现实生活中使用此代码。像BioPython这样的图书馆做得更好。例如,您展示的代码将生成相当庞大的数据结构。由于FASTA文件有时非常大,可能会耗尽内存。而且,它不能以最好的方式处理可能的异常。在</p>
<p>最后,使用库只会节省您的时间。在</p>
<p>BioPython示例代码:</p>
^{pr2}$