Question: HTseq-counts not counting correctly?
gravatar for vm.higareda
6 months ago by
vm.higareda20 wrote:


I am doing a DE analysis of Drosophila, after mapping my fasq file to Drosophila genome with hisat2. I am using htseq-count with the gtf file to get the count table.

In the count table, I have the identifier corresponding to 28S rRNA and 18S rRNA, but it has ZERO counts, that can not be possible.

Why I have zero counts? I know that I have ribosomal sequences in my treatments

rna-seq drosophila htseq • 228 views
ADD COMMENTlink written 6 months ago by vm.higareda20
gravatar for genomax
6 months ago by
United States
genomax65k wrote:

By default counting programs ignore multi-mapping reads which are bound to be the case with rRNA reads (see FAQ on this page).

If you want to count them then use this option --non-unique=all :

--nonunique=<nonunique mode>

    Mode to handle reads that align to or are assigned to more than one feature in the overlap <mode> of choice (see -m option). <nonunique mode> are none and all (default: none)
ADD COMMENTlink modified 6 months ago • written 6 months ago by genomax65k

and most common RNAseq protocols will likely do an ribosomal RNA depletion step as well (== to exclude those from the sequencing)

ADD REPLYlink written 6 months ago by lieven.sterck4.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1515 users visited in the last hour