featureCounts parameters for RNA-Seq
1
0
Entering edit mode
22 months ago

Dear Community,

I have been running featureCounts to count mapped reads in RNA-Seq paired-end data with the following parameters

featureCounts -T 10  -t exon -s 2 -p -g gene_id -a <path_to_annotation.gtf> -o <path_to_outoput_directory> mapped.bam

as mentioned in http://bioinf.wehi.edu.au/featureCounts/

However, I came across two posts where the parameters

and

have been used.

This got me a little confused and will like to to know if and when I should use these parameters as well as what will be the difference in the output from the default parameters?

Many thanks in advance

featurecounts RNA-Seq Subread ReadCounts • 2.0k views
ADD COMMENT
0
Entering edit mode

Hey do you get some answers for this? I am also using featurecounts now, paired ends RNAseq, and I am not sure if i have to use specify strandness? Could you please help me ?

ADD REPLY
0
Entering edit mode

If your sequencing is stranded then you need to specify to featureCounts which strand has been sequenced. If your sequencing is not stranded then you don't. Pretty logical!

ADD REPLY
0
Entering edit mode

Thanks for your answer. May I ask do you reccomand --primary to be set on?

ADD REPLY
0
Entering edit mode

Generally, we recommend that featureCounts be run with default settings, which does not include --primary. Just follow the documentation. The specific use depends on your intended application.

ADD REPLY
0
Entering edit mode

Thanks so much Gordon Smyth. In default setting -t exon is set on. So when is -t gene reccomanded? In my bacteria gff file, I do not see exon, but in my fungi gff file, I do. I wonder can I just set -t gene for both bacteria and fungi? I really appreciate your help. Best regards, Yanfang

ADD REPLY
0
Entering edit mode
12 weeks ago

Honestly, if you aren't using something like Kallisto, which outputs counts more or less natively, I'd be using RSEM, not featureCounts RSEM is smarter about dealing with ambiguously assigned reads.

ADD COMMENT

Login before adding your answer.

Traffic: 1449 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