Question: HTseq-counts not counting correctly?
gravatar for vm.higareda
2.4 years 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 • 802 views
ADD COMMENTlink written 2.4 years ago by vm.higareda20
gravatar for GenoMax
2.4 years ago by
United States
GenoMax96k 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 2.4 years ago • written 2.4 years ago by GenoMax96k

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

ADD REPLYlink written 2.4 years ago by lieven.sterck10k
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: 1692 users visited in the last hour