How to calculate Exon-Intron junction perbase coverage?
Hi all,
I want to extract the Exon-Intron junction reads from the paired end RNA-Seq data. I’m only considering the reads that are crossing Exon-Intron boundary. For that I constructed a bed file flanking 100 nt up and downstream of the Exon-Intron junction. But I wanted to make sure that I have enough reads are supporting the junction. So I wanted to count the per-base frequency for the junctions. I’m considering only if I have >5x coverage. How can i calculate the per-base frequency? Is there any flaws in my strategy. if yes, please fell free to address.
Thanks.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
is there any way to normalize the counts and compare it for the splicing quantification?
written by PK •
-
Identify reads that span intron-exon junctions in RNA-Seq
written by oriolebaltimore
-
STAR using splice junctions from SJ.out to find flanking exons (GRCh37)
written by stefanos.bamopoulos •
-
How to extract reads supporting protein-coding gene splice junctions from STAR's SJ.out.tab output
written by Nicolas Rosewick
-
how to combine junction.bed files with different rows into one table
written by shirley0818
-
Count splice junction reads from SAM/BAM files
written by juncheng
-
Selecting and de-selecting exon-exon junctions based on depth or coverage from RNA-seq mapping
written by mjoyraj
-
Splice Junction file intersection with genome annotation
written by ruchiksy •
-
Question About Minor Isoform Filtering Of Tophat
written by Fan •
-
Intron Retention And Alternative 5'3' Splice Site Identification
written by Arun
