fq文件随机提取sys.argv[1]条reads python

来源:互联网 发布:centos pcre 下载 编辑:程序博客网 时间:2024/05/18 13:08
import random
import sys
fin=open('reads.fq','r').readlines()
fout=open('out.fq','w')

a=fin[::4]
b=fin[1::4]
c=fin[2::4]
d=fin[3::4]
lt=list(zip(a,b,c,d))
n=[i for i in range(len(lt))]
m=sorted(random.sample(n,int(sys.argv[1])))
for i in range(len(m)):
    fout.write(''.join(lt[m[i]]))
原创粉丝点击