使用GFF3文件分析外显子及其数量和内含子的大小

    with open('TSG.gff3', 'r') as f:
    for line in f:
        lin = line.strip().split('\t')
        name = lin[8].split(';')[1].split('=')[-1]
        if lin[2] == 'gene':
            print '\n',
            print name, lin[0], lin[3], lin[4],lin[6],
        if lin[2] == 'exon':
            print lin[3], lin[4],
    
    # 上述程序统计每个基因的外显子起始与结束的位置，保存为1.txt,注意需要打开1.txt编辑删除第一行的空行
    with open('1.txt', 'r') as f:
    for line in f:
        lin = line.strip().split()
        a = len(lin)
        for i in range(6, a, 2):
            exon = int(lin[i]) - int(li