gff3格式转换

使用gmap软件输出的gff3文件并不能直接用来作为gffcompare的输入文件
故此，需要进行转换，转换也简单就是讲target字段删除即可。

---

#!/usr/bin/env python
# -*- coding: utf-8 -*- 
__author__ = 'shengwei ma'
__author_email__ = '[email protected]'

with open('LYL_5A.gff3', 'r') as f:
    for line in f:
        line1 = line.strip().split('\t')
        if line1[2] == 'gene':
            print line,
        if line1[2] == 'mRNA':
            print line,
        if line1[2] == 'exon':
            new = line1[8].split(';')
            print '%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s;%s;%s\n' %(line1[0], line1[1], line1[2],line1[3],line1[4],line1[5],
                                                           line1[6], line1[7], new[0], new[2], new[2]),