awk -f script.awk original.gtf >  genes.gtf
sed -i 's/gene_biotype/gene_type/' genes.gtf
sed -i 's/transcript_biotype/transcript_type/' genes.gtf
sed -i 's/three_prime_utr/UTR/' genes.gtf
sed -i 's/five_prime_utr/UTR/' genes.gtf
#awk '{print $0"  "}' genes.gtf >test.gtf && rm genes.gtf && mv test.gtf genes.gtf
grep  tRNA genes.gtf >trna.gtf


#removed exon
#I       PomBase exon    2562374 2562427 .       -       .       exon_number "1"; gene_id "SPSNRNA.06"; transcript_name "snu6"; transcript_source "PomBase"; gene_source "PomBase"; gene_name "snu6"; tag "basic"; exon_id "SPSNRNA.06.1:exon:1"; transcript_id "SPSNRNA.06.1"; gene_type "snRNA"; transcript_type "snRNA"
