How can I use featurecounts after generating a bam file using BWA?
1
0
Entering edit mode
5.5 years ago
a.rex ▴ 350

I have a reference transcriptome, to which I have aligned reads using BWA. My goal is to quantify the aligned transcripts - how can I do this? Can I use featurescounts - bearing in mind that I have aligned to a transcript fa file.

RNA-Seq features counts • 2.8k views
ADD COMMENT
0
Entering edit mode

I read yesterday about Salmon and Kallisto, which can generate count metrics from a transcripts.fa file and the FQ reads. Maybe they have an option to use an aligned BAM?

ADD REPLY
2
Entering edit mode

Salmon does, see the section of the manual :)

ADD REPLY
3
Entering edit mode
5.5 years ago
h.mon 35k

If you aligned to the transcriptome, you should not use featureCounts, particularly if you have a well annotated transcriptome, or assembled the transcriptome yourself. In both cases, the transcriptome should have a lot of isoforms, which would result in a lot of multi-mapped reads, which would be discarded by featureCounts. See Why you should use alignment-independent quantification for RNA-Seq for a discussion on this.

Use RSEM in this case, or as ATpoint noted above, salmon quant.

ADD COMMENT

Login before adding your answer.

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