- bf = pysam.AlignmentFile(sys.argv[1], 'rb')
- results = pysam.AlignmentFile("deduplex.bam","wb",template=bf)
- readed = []
- for r in bf:
- name = r.query_name.split('|')[1]
- UMI1 = re.findall(r'.{8}', name)[0]
- UMI2 = re.findall(r'.{8}', name)[1]
- tag1 = UMI1,UMI2,r.pos,r.pnext
- tag2 = UMI2,UMI1,r.pos,r.pnext
- if tag1 not in readed and tag2 not in readed:
- results.write(r)
- readed.append(tag1)
- else:
- continue
- results.close()
- bf.close()
用append,列表是[(),()]这种格式
用extend,就全部连在一起了
请问一下有没有什么方法可以提高运行效率啊?


雷达卡




京公网安备 11010802022788号







