pysam筛选reads写bam
使用pysam从bam中筛选readsbamfile = 'Final.bam'bam = pysam.AlignmentFile(bamfile, "r")outfile = pysam.AlignmentFile('out.bam', "wb", template=bam)ALL = set()for reads in bam: status = '''some filters''' if status: outfile.write(reads) print
原创
2021-07-10 16:37:13 ·
923 阅读 ·
0 评论