splitname

GSVIVT01000001001 pacid=17816625 transcript=GSVIVT01000001001 locus=GSVIVG01000001001 ID=GSVIVT01000001001.Genoscope12X annot-version=Genoscope.12X
变为
GSVIVT01000001001

##splitname.py for python2
#!/usr/lib/python
#-*- coding:utf-8 --*-

import sys, os
def splitname(infile, outfile):
        ini = open(infile, 'r')
        out = open(outfile, 'w')

        for line in ini:
                if line[0]=='>':
                        lineL=line.split(' ')
                        print >>out, lineL[0]
                else:
                        print >>out, line.strip()
if __name__=="__main__":
        infile = sys.argv[1]
        outfile = sys.argv[2]
        splitname(infile, outfile)
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容