Receiving more reads from exons than whole genes using featurecounts
Entering edit mode
18 months ago


In order to receive intronic counts for my rnaseq data, I wanted to substract the reads for exons from the total reads of the respective genes. For this, I am using featureCounts from the RSubread package in R. However, when I substracted the reads, some genes were returning negative counts (meaning there are more reads in the exons than in the whole gene). As a .gtf file I am using the standard mm10 annotation from ensembl with the standard options with GTF.featureType set to "gene" or "exon" respectively, except that I set fracOverlap to 1 (to remove any intronic reads from the exon reads). Does anybody have an idea how this can happen?

Thanks in advance

rna-seq featurecounts RSubread • 373 views

Login before adding your answer.

Traffic: 1577 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6