How to count novel and known splice junction reads from BAM
1
0
Entering edit mode
9.5 years ago

Hi,

I have BAM file from TopHat 2 alignment with RefSeq genes annotation. Is there any way to count separately the exon junction reads coming from RefSeq transcripts and novel junctions (not from RefSeq).

Could somebody help me out?

RNA-Seq • 2.9k views
ADD COMMENT
0
Entering edit mode
9.5 years ago
Jordan ★ 1.3k

I have never tried this. But may be you can use bamtobed tool to convert the regions in bam to bed format. And then find common regions between this bed file and the gtf file.

You can may be try the following steps.

Step1: Convert bam to bed
bedtools bamtobed -i reads.bam > reads.bed

Step2: Convert gtf to bed
gtf2bed < refseq.gtf > refseq.bed

Step3: Compare the common regions between the two bed files
bedops --everything reads.bed refseq.bed
ADD COMMENT

Login before adding your answer.

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