Question: How can I use featurecounts after generating a bam file using BWA?
gravatar for a.rex
8 days ago by
a.rex170 wrote:

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 • 85 views
ADD COMMENTlink modified 8 days ago by h.mon23k • written 8 days ago by a.rex170

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 REPLYlink written 8 days ago by RamRS20k

Salmon does, see the section of the manual :)

ADD REPLYlink written 8 days ago by ATpoint13k
gravatar for h.mon
8 days ago by
h.mon23k wrote:

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 COMMENTlink written 8 days ago by h.mon23k
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: 1795 users visited in the last hour