RSEM BAM outputs, which one to use for RSeQC?
0
0
Entering edit mode
2.8 years ago
Freek ▴ 10

Hi,

I have used RSEM to estimate Transcript abundances in a sample. RSEM uses STAR to do the mapping (if asked) and optionally can output a BAM file where the reads are aligned against the genome and not the transcriptome (this genome file ends in .STAR.genome.bam and is generated with the option "--star-output-genome-bam "). Now I was wondering, I want to use RSeQC to look at various quality related aspects of my BAM file, but which one should I use as input? The one aligned against the transcriptome (.transcript.bam) of the one aligned against the genome (.STAR.genome.bam)?

Highest regards,

Freek.

RSEM STAR RSeQC BAM • 1.2k views
ADD COMMENT
1
Entering edit mode

Did you ever find the answer to this question? I have been used genomic BAM files when processing data through RSeQC, but I wanted to double check whether this was the correct approach.

ADD REPLY
0
Entering edit mode

Hi, thanx for coming back at this. 16 months later, 16 months wiser :)

I use the genome BAM file, although I never tested the transcriptome BAM file, I'm assuming, since the transcriptome BAM was created with the transcriptome as a reference, that there are no reads mapped to anything other than exons (valid transcripts) in said BAM file, this would mean that RSeQC can not determine the number of intronic reads or reads on intergenic regions etc. Since RSeQC uses it's own gene/transcript definitions (I downloaded the BED file they offer), I think it only makes sense to apply RSeQC to a genome BAM.

Still, I might just feed RSeQC the transciptome BAM once and see what happens, in the days to come :)

ADD REPLY

Login before adding your answer.

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