Question: How can I use featurecounts after generating a bam file using BWA?
0
gravatar for a.rex
6 months ago by
a.rex190
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
2

Salmon does, see the section of the manual :)

ADD REPLYlink written 6 months ago by ATpoint21k
3
gravatar for h.mon
6 months ago by
h.mon27k
Brazil
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.

Help
Access

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