Question: exon counting with featureCounts
0
gravatar for inah
22 months ago by
inah0
inah0 wrote:

Hi,

  I want to obtain read counts at the exon level using featureCounts. I run Rsubread and use these options:

annot.ext="/home/inah/RefGTF/GRCh38/annotation/Homo_sapiens.GRCh38.85.gtf",
isGTFAnnotationFile=TRUE,
GTF.featureType="exon",  GTF.attrType="exon_id",  useMetaFeatures=FALSE

The resulting count matrix has 1,182,163 rows but only 678,276 unique row identifiers, so there are multiple rows for the same exon-id.

Of course when I set useMetaFeatures=TRUE, then the resulting count matrix has 678,276 rows all with unique identifiers.

Sorry for the ignorant question, but why do I have so many more "exons" than "exon-ids"?

Thanks, Ina

rnaseq • 892 views
ADD COMMENTlink written 22 months ago by inah0

You should add "Rsubread" as a tag to your question, otherwise the Rsubread author might not be alerted to your question.

ADD REPLYlink written 22 months ago by Gordon Smyth37k
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 16.09
Traffic: 115 users visited in the last hour