This is a test version of Biostars. For the public version, visit https://www.biostars.org.
htseq gene atrribute error

$ htseq-count -q --additional-attr 534accepted_hits.bam 535accepted_hits.bam Rsolani.gff Error occured when processing GFF file (line 10 of file Rsolani.gff): Feature exon-AG1IA_00001-2 does not contain a 'gene_id' attribute [Exception type: ValueError, raised in count.py:79] what should I do now ?

software error next-gen rna-seq
• 1,748 views • link
updated by Strand NGS • Strand NGS • written by Balatheskulter • Balatheskulter

1 answer

The gff3 file you are passing doesn't have the 'gene_id' attribute. Please do check the file and see the additional options in htseq to pass the correct id.

Some of other posts here that you can refer to:

htseq-count error Feature does not contain a 'gene_id' attribute

htseq-count error: "does not contain 'gene_id'

and many more are there available here to help.

• 0 views • link
written by Strand NGS • Strand NGS

Log in to answer this question.

More posts like this
  • HTseq-count error
    written by sanatbhadsavle • sanatbhadsavle
    Hi, I am trying to count miRNA reads from a set of sam files. I am using HTseq-count to do so. I downloaded the mouse …
  • htseq-count returns : does not contain a 'gene' attribute
    written by ziziqolo • ziziqolo
    Dear BIOSTAR community, I'm trying to make count matrix with htseq-count, htseq-count -s yes -t gene -i gene 01.sorted.sam annotation_cattle.gff > 01.txt even with `--idattr=gene` …
  • cant find gene_id
    written by meisam.radfar • meisam.radfar
    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 …
  • Query on htseq count
    written by lavanyac790 • lavanyac790
    I am trying to run htseq-count for carrying out rna-seq analysis for solanum tuberosum and i used the following command: htseq-count --format bam --order pos …
  • Query on htseq count
    written by lavanyac790 • lavanyac790
    I am trying to run htseq-count for carrying out rna-seq analysis for solanum tuberosum and i used the following command: htseq-count --format bam --order pos …
  • HTSeq error about gene_id
    written by felipead66 felipead66
    ##gff-version 3 Pav_co3988703.1 augustus gene 1 195 0.7 - . ID=Pav_co3988703.1_g010.1.br;Note=PREDICTED: disease resistance protein RPM1-like Pav_co3988703.1 augustus exon 16 195 . - . ID=Pav_co3988703.1_g010.1.br:exon:1; Parent=Pav_co3988703.1_g010.1.br …
  • Error with HTseq RNAseq read count
    written by Bioinfonext Bioinfonext
    Hi, I am getting error while running HTseq, could you please suggest what is the issue; I am using this command $ htseq-count Oryza_indica.ASM465v1.44.chr.gff3 Leaf_T1_F_R10_S1_L001.sam …
  • htseq-count error Feature does not contain a 'gene_id' attribute
    written by madzayasodara • madzayasodara
    I'm trying to do read count using htseq-count > htseq-count --mode=union --stranded=no --idattr=gene_id -f sam > ./filteredfirstZBND1X_gff.sam ./filteredsecZBND1X_gff.sam > /home/madzays/data/finc > h_data/ZF_transcriptome/taeGut2_gff/taeGut2.gff > > ZBND1X_gff_counts.sam …
  • htseq-count error: "does not contain 'gene_id'
    written by lynne.m.a.harris • lynne.m.a.harris
    I'm attempting to use htseq-count to get read counts for some *Plasmodium *RNA-seq data. The command line looks as follows: htseq-count --format=bam --order=pos 3D7_1-1.bam PlasmoDB-28_Pfalciparum3D7.gff …
  • Converting gff file to gtf for htseq-count
    written by natsterbug • natsterbug
    After running TopHat2/2.1.0 on RNA-seq SE 50bp reads from *S.tuberosum*, I am now attempting to count reads mapping to each feature using htseq-count. Using the …