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.
0 answers
No answers yet.
Log in to answer this question.