This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get The Idea Of Splicing From Reads Mapped In Rna-Seq

I've got a set of 100 bam files from a public experiment, I want to have an idea of splicing in each of them regarding three exons,without entering in some kind of depth-level procedure like Cufflinks or DEXSeq,

Lets say that my exons are named 1,2 and 3, and I want to know in how many samples I have a splicing event of the number two, so i was looking in the threads and I found that using coverageBed with my bed file of the three exons I could get some kind of idea per bam file

coverageBed -split -abam my_alignment -b exons_to.bed

Am I correct?

I was also thinking of getting the reads mapped in flanking end positions of read 1 and start of read 3 with samtools

What do you think about it? Any idea will be kindly appreciated

Thanks in advance!

rnaseq splicing bedtools samtools
• 3,280 views • link
updated by Charles Warden Charles Warden • written by antgomo • antgomo

2 answers

Use bwa mem then sort the reads by name and you can easily get reads that map across junctions.

• 0 views • link
written by Istvan Albert Istvan Albert

I would recommend MATS for identifying splicing events (for example, exon skipping - for exon #2 in your case, I believe):

http://rnaseq-mats.sourceforge.net/

MISO is probably a more popular option, but I think MATS is more specific (and the result table is easier to read):

http://genes.mit.edu/burgelab/miso/

If you have long and/or paired end reads, I would also recommend viewing splicing junctions with a sashimi plot in IGV

• 135 views • link
written by Charles Warden Charles Warden

Log in to answer this question.

More posts like this
  • Java Samtools Api: Write Bam To Stdout
    written by dariober dariober
    <p>Hello,</p> <p>I'm trying to use the Java API for picard/samtools to write a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file to standard output.</p> <p>With the code below I can …
  • Pcr Free Variant Calls
    written by tommy.carstensen • tommy.carstensen
    <p>If I have a PCR free high coverage <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> and a low coverage <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> and I downsample the PCR free high coverage <a …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • Samtools View Of A Sam File Entire Chromosome
    written by dfernan dfernan
    <p>Hi, </p> <p>I have (what I think) is a valid <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file.</p> <p>When I do:</p> <pre><code>samtools view -S F1i1F.C57.2.Aligned.out.sam chr19_C57 </code></pre> <p>I get the …
  • Length Of Contigs In Bam File Off By 1 Versus Reference Fasta
    written by Tancata • Tancata
    <p>I have been trying to call some variants with GATK. It complains, rightly, that my the length of the contigs in my <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file …
  • Reduce Bam File With Specific Loci In Bedfile
    written by marcodpc • marcodpc
    <p>I want to reduce my <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file, obtaining a new <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file that contain </p> <p>only loci in the bed file in input …
  • How To Interpret The Flag Field In The Sam Specification
    written by deepthithomaskannan deepthithomaskannan
    <p>Hi all ,</p> <p>I created a <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and converted it to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file using <a href='http://samtools.sourceforge.net/'>samtools</a>.</p> <p>One of the record in <a …
  • Problem With Sorting Bam File: Samtools
    written by lyfsa • lyfsa
    <p>Hi!</p> <p>I have three aligned pair end reads in <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> format. Two of them are of same sample run in different lanes. I first …
  • Sam To Bam Conversion Reference Fasta File
    written by Atom Smasher • Atom Smasher
    <p>Hello,</p> <p>I had some data in "eland_export.txt" format and I want to convert it to the <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> format first and then later to the …
  • Ucsc/Ensembl Bam Converter?
    written by 2184687-1231-83- 2184687-1231-83-
    <p>Is there any tool to convert <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files that have been produced against Ensembl genome references to UCSC-compatible <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files? And viceversa? AFAICS, …