Question: FeatureCounts count reads per gene (Including Introns)
gravatar for paulranum11
6 months ago by
paulranum1140 wrote:

Hi Everyone,

I am using featureCounts to count the number of reads mapping to genes.

        featureCounts -a $myGTF \
                  -o counted_output \
                  -R BAM Aligned.sortedByCoord.out.bam \
                  -T $NUMCORES \

The above bash command works counting reads that map to the exonic sequences of each genes annotated in the $myGTF file.

My question is, how can I get featureCounts to also count reads mapping to intronic sequences? It is unclear to me how to achieve this from the documentation but it might be through some combination of -O and -f.


rna-seq featurecounts • 654 views
ADD COMMENTlink written 6 months ago by paulranum1140

If you can get/create an annotation file with intron features, you can then count using that file with -t intron option.

ADD REPLYlink modified 6 months ago • written 6 months ago by genomax76k
gravatar for ATpoint
6 months ago by
ATpoint28k wrote:

Alternatively make a custom reference in featureCount's SAF format (see manual for details) to include introns.

ADD COMMENTlink written 6 months ago by ATpoint28k

This is what i ended up doing. It was easy to make this SAF format file using biomart from For future readers the columns to include for this format are: GeneName, Chromosome, GeneStart, GeneEnd, Strand

Note that the column names header needs to be removed before the SAF file is used with featureCounts.

ADD REPLYlink written 6 months ago by paulranum1140
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: 827 users visited in the last hour