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

Hi,

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
ADD COMMENT

Login before adding your answer.

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

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

Powered by the version 2.3.6