Question: Use BAM file created from subjunc function (Rsubread package) for feature counts?
gravatar for Stephen Piccolo
3.3 years ago by
United States
Stephen Piccolo530 wrote:

I'm able to process RNA-Seq data using Rsubread to obtain feature counts after using the align function. I can also identify exon-exon junctions and fusions after aligning using the subjunc function. Is it acceptable (to save processing time) to apply featureCounts() to the subjunc BAM files? Will this produce different feature counts than I would get if I used the align function instead?

ADD COMMENTlink modified 3.3 years ago by Wei Shi2.9k • written 3.3 years ago by Stephen Piccolo530
gravatar for Wei Shi
3.3 years ago by
Wei Shi2.9k
Wei Shi2.9k wrote:

You will get slightly different counts. Subjunc is more stringent in read mapping in order to control errors in calling exon-exon junctions. If your analysis requires the identification of exon-exon junctions, then you should use Subjunc.

We typically use Subread aligner (align()) in our gene expression analysis pipeline. But there should be only minor difference in the list of differentially expressed genes identified from using the two aligners. The two aligners use the same mapping paradigm (seed and vote) to align reads.

ADD COMMENTlink written 3.3 years ago by Wei Shi2.9k
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.2.0
Traffic: 205 users visited in the last hour