This is a test version of Biostars. For the public version, visit https://www.biostars.org.
htseq-count overcounting non-unique alignments

Hi,

I have some reads that I aligned with STAR aligner version 2.4.2a and then ran through htseq-count version 0.6.1.

I'm noticing a discrepancy between what STAR reports as "Number of reads mapped to multiple loci" (7.1M, ~11% of reads) and what htseq-count reports as __alignment_not_unique (22.3M, ~31% of reads). This is a pretty big discrepancy. I ran the same thing on multiple files, both single-end and paired-end, and while the size of the discrepancy varies, it looks like it's always there.

I was wondering if anyone could help me figure out what's causing this. Shouldn't HTSeq be taking that count directly from the bam file?

Thanks,

Jake

htseq-count rna-seq star
• 7,139 views • link
updated by Ram Ram • written by jc.szamosi • jc.szamosi

What about "% of reads mapped to too many loci"? This is different from "% of reads mapped to multiple loci".

• 0 views • link
updated by Ram Ram • written by Ashutosh Pandey Ashutosh Pandey

1 answer

Hi Jake,

The __alignment_not_unique is rather an alignment count than an read count. If you sum up the gene counts, no features, and ambiguous counts, you'll get very very close to your uniquely mapping reads form the STAR report.

Cheers,
Michael

• 0 views • link
updated by Ram Ram • written by michael.ante michael.ante

So what you're saying is that STAR counts each multiply-mapped read once, but HTSeq counts them as many times as they're mapped? That makes sense. Thank you!

• 0 views • link
written by jc.szamosi • jc.szamosi

Exactly, that's my experience.

• 0 views • link
written by michael.ante michael.ante

Log in to answer this question.

More posts like this