Question: HTseq-counts not counting correctly?
gravatar for vm.higareda
10 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 • 341 views
ADD COMMENTlink written 10 months ago by vm.higareda20
gravatar for genomax
10 months ago by
United States
genomax70k 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 10 months ago • written 10 months ago by genomax70k

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

ADD REPLYlink written 10 months ago by lieven.sterck5.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: 1316 users visited in the last hour