Question: how to combine junction.bed files with different rows into one table
0
gravatar for shirley0818
4.0 years ago by
shirley081890
United States
shirley081890 wrote:

Hi All,

I am working on ~100 samples for detecting alternative splicing. TopHat generates a junction.bed file for each sample. However, each of these bed files has different number of rows, and the coordinates of each junction is not same across samples. I think this junction.bed file includes known and novel junctions.

Since I am only interested in known junctions in Ensembl annotation database, how can I map these 100 junction.bed files to Ensembl gtf file and obtain a table matrix with the row as known-junction and column as sampleID?

Or do I need to create a exon-exon junction annotation bed file from Ensembl, then apply RSeQC to obtain reads for each junction against mapped .bam files?

Many thanks,

Shirley

Shirley

 

 

 

rna-seq • 1.4k views
ADD COMMENTlink written 4.0 years ago by shirley081890
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: 1073 users visited in the last hour