擅长:python、mysql、java
<p>当您读取fasta文件时,您可以创建一个类似于掩码的列表:</p>
<pre><code>with open('mask.txt') as mf:
mask = [ s.strip() == '1' for s in mf.readlines() ]
</code></pre>
<p>然后:</p>
^{pr2}$
<p>或者:</p>
<pre><code>from itertools import izip
for b, line in izip(open(mask_file), open(seq_file)):
if b.strip() == '1':
*something* line
</code></pre>