This is a test version of Biostars. For the public version, visit https://www.biostars.org.
cant find gene_id

i want to use htseq-count for counting reads of my RNA-seq project

i use GFF3 file ( wheat ) from plantensemble website but when i write this command:

htseq-count -q -i gene_id -f bam X.bam Triticum_aestivum.IWGSC.49.gff3

i see this:

Feature transcript:ENSRNA050013875-T1 does not contain a ‘gene_id’ attribute [exception type: ValueError, raised in features.py.329]

i use ID=gene but this is not work too.

in my gff3 file i van see:

ID=gene:ENSRNA050013875 gene_id=ENSRNA050013875 ID=transcript:ENSRNA050013875-T1 parent=gene:ENSRNA050013875 transcript_id=ENSRNA050013875-T1 parent=transcript:ENSRNA050013875-T1 Name:ENSRNA050013875-E1 exon_id=ENSRNA050013875E1

just Name and exon_id work but i dont need them.

rna-seq software error

0 answers

No answers yet.

Log in to answer this question.