small RNA seq筛选reads长度
small RNA seq筛选reads长度,以筛选18~25 nt为例pythonfw = open('out', 'w')with open('inp', 'r') as fr:content = fr.readlines()for i in range((len(content) + 1)/4):site = 4 * iif ...
·
small RNA seq筛选reads长度,以筛选18~25 nt为例
python
fw = open('out', 'w')
with open('inp', 'r') as fr:
content = fr.readlines()
for i in range((len(content) + 1)/4):
site = 4 * i
if 16 < len(content[1 + site]) < 25: #保留长度大于17,小于26的reads
fw.write(content[0 + site] + content[1 + site] + content[2 + site] + content[3 + site].strip() + '\n')
更多推荐
已为社区贡献4条内容
所有评论(0)