Question: How can I use featurecounts after generating a bam file using BWA?
gravatar for a.rex
6 months ago by
a.rex190 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 • 243 views
ADD COMMENTlink modified 6 months ago by h.mon27k • written 6 months ago by a.rex190

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 6 months ago by RamRS23k

Salmon does, see the section of the manual :)

ADD REPLYlink written 6 months ago by ATpoint21k
gravatar for h.mon
6 months ago by
h.mon27k 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 6 months ago by h.mon27k
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: 616 users visited in the last hour