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()
正則匹配中的search和sub用法
最后編輯于 ：
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明：文章內容（如有圖片或視頻亦包括在內）由作者上傳并發布，文章內容僅代表作者本人觀點，簡書系信息發布平臺，僅提供信息存儲服務。
平臺聲明：文章內容（如有圖片或視頻亦包括在內）由作者上傳并發布，文章內容僅代表作者本人觀點，簡書系信息發布平臺，僅提供信息存儲服務。
- 文/潘曉璐 我一進店門，熙熙樓的掌柜王于貴愁眉苦臉地迎上來，“玉大人，你說我怎么就攤上這事。” “怎么了？”我有些...
- 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著，像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發上，一...
- 文/蒼蘭香墨 我猛地睜開眼，長吁一口氣：“原來是場噩夢啊……” “哼！你這毒婦竟也來了？” 一聲冷哼從身側響起，我...
