Question: HT-seq-count default parameter to "score"secondary alignments contradicts its documentation that it ignores multi mappers by default. Not clear!
gravatar for Zee_S
9 weeks ago by
Zee_S50 wrote:

Dear Biostars Community,

Im trying to quantify RNAseq data using HTseq-count to get number of reads per gene. In the documentation it says that HTseq by default ignores multimapping reads. But in the manual ( the default for secondary alignments is set to "score" and not "ignore". I don't understand what this means. If it ignores multi mappers by default why is the default parameter for counting multimappers set to "score?"

Also, if I'm not using htseq-count for differential expression analyses but just to look at the expression level of all genes in my model organism, should I include multimapping reads in the gene counts? or not?

Im really confused and would be grateful for your insights on this issue.

Thanks a lot for your help.

ADD COMMENTlink written 9 weeks ago by Zee_S50

I guess it's only "scored" if you set the parameter "--nonunique" to "all". Afterwards, you can count the secondary and supplementary alignments.

ADD REPLYlink written 9 weeks ago by michael.ante3.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: 1705 users visited in the last hour