This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Sbicolor_255_v2.1.gene_exons.gff3 to Sbicolor_255_v2.1.gene_exons.gtf conversion

Hello All,

I have trouble converting Sbicolor_255_v2.1.gene_exons.gff3 file to Sbicolor_255_v2.1.gene_exons.gtf. I have used gffread of cufflinks to do the conversion. During the conversion process, gffread changed all the genes to exons (example is listed below). As a result, HTSeq couldn't count the transcripts as the gene information is missing. Please suggest a best way to fix this problem.

Sbicolor_255_v2.1.gene_exons.gff3

Chr01    phytozomev10    gene    1619    2809    .    +    .    ID=Sobic.001G000100.v2.1;Name=Sobic.001G000100
Chr01    phytozomev10    mRNA    1619    2809    .    +    .    ID=Sobic.001G000100.1.v2.1;Name=Sobic.001G000100.1;pacid=28395113;longest=1;Parent=Sobic.001G000100.v2.1

Sbicolor_255_v2.1.gene_exons.gtf

Chr01    phytozomev10    exon    1619    2809    .    +    .    transcript_id "Sobic.001G000100.1.v2.1"; gene_id "Sobic.001G000100.v2.1"; gene_name "Sobic.001G000100";
Chr01    phytozomev10    CDS    1951    2616    .    +    0    transcript_id "Sobic.001G000100.1.v2.1"; gene_id "Sobic.001G000100.v2.1"; gene_name "Sobic.001G000100";

Thanks

rna-seq
• 1,643 views • link
updated by Ram Ram • written by sandeepmarla • sandeepmarla

0 answers

No answers yet.

Log in to answer this question.

More posts like this