python7

2021-05-23  本文已影响0人  rong酱
# -*- coding: utf-8 -*-

import os
import fileinput

keylist = []
bamlines = open("n1.filter.sam","rb").readlines()
for lines in bamlines:
    line = lines.strip().split("\t")
    for liness in line:
        if 'CR' in liness:
            barcode = liness.strip().split(":")[-1]
            name = "n1.filter.TAG_CB_" + str(barcode)
            if barcode not in keylist:
                keylist.append(barcode)
                os.system('touch %s'%(name))
                content = open(name,"w")
                content.write(lines)
            else:
                content.write(lines)

不念往事,立足现在

上一篇 下一篇

猜你喜欢

热点阅读