qq_26012913的博客

首先我们要把gtf文件中的exon抓取出来grep "exon" genome.gtf > genome_exon.gtf然后提取genome_exon.gtf文件中的gene的exon的长度和得到我们想要的gene的长度python count_genelen_from_gft.py genome_exon.gtf gene.len这其中count_genelen_from_gft.py的代码如下：import sys,refile1 = sys.argv[1]file2 = sy