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

Hello

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
0
gravatar for genomax
6 months ago by
genomax65k
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.

Help
Access

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