使用gmap软件输出的gff3文件并不能直接用来作为gffcompare的输入文件
故此,需要进行转换,转换也简单就是讲target字段删除即可。
#!/usr/bin/env python
# -*- coding: utf-8 -*-
__author__ = 'shengwei ma'
__author_email__ = 'shengweima@icloud.com'
with open('LYL_5A.gff3', 'r') as f:
for line in f:
line1 = line.strip().split('\t')
if lin