正则匹配中的search和sub用法

import os
import re
import sys

try:
    infile = sys.argv[1]
    outfile = sys.argv[2]
except:
    print "Usage: python filter_blast.py infile outfile"

gff=file(infile)
newgff=file(outfile,"w")
j=0

hang=10000
genechr=''
num=1
for line in gff:
    j+=1
    line=line.strip().split("\t")
    if re.search('Chr', line[0]) and re.search('gene', line[2]):
        gene=line[8].split("=")[-1]
        #print gene
        chr=re.sub('Chr','Gm',line[0])
        chrn=re.sub('Chr','',line[0])
        newgene=chr+'g'+str(hang)
        hang+=1
        newgff.write(chr+"\t"+line[3]+
        "\t"+line[4]+"\t"+line[6]+"\t"+gene+"\t"+newgene+"\n")
print j,"records read"

gff.close()
newgff.close()
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 看个小软文感觉挺适合自己励志一下! 我经常在微信群里吹捧自己 昵称都是 贺静 帅无底限 残无人道! 有些男人就看不...
    e0c9b0c45626阅读 300评论 3 5
  • 如果没有十分喜欢我,就放手
    Amyuer阅读 181评论 0 0
  • 文/伊人若雨春雾蒙蒙雾里看花春风徐徐杨柳依依春雨洁净洗去烟尘沉睡的种子已苏醒在酥软的泥土里萌芽呵开大地之门召唤深冬...
    伊人若雨阅读 277评论 0 1
  • 今天几乎下了一天的雨,但仍然妨碍不了会员参加轰趴的热情。轰趴非常好的联系了会员们的情感,也挖掘很多人的咪霸潜能,本...
    米特Dragon阅读 253评论 1 0