Question: FeatureCounts count reads per gene (Including Introns)
0
gravatar for paulranum11
6 months ago by
paulranum1140
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 \
                  -M

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.

Thanks!

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
1
gravatar for ATpoint
6 months ago by
ATpoint28k
Germany
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 ensembl.org. 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.

Help
Access

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