Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
mapping junctions.bed into exons/introns microsat RNA Sequencing 0 03-06-2014 09:18 AM
Tophat junctions.bed RockChalkJayhawk RNA Sequencing 7 12-12-2013 11:56 AM
tophat junctions.bed file upendra_35 RNA Sequencing 2 10-23-2012 08:17 AM
convert junctions.bed to .juncs xy6699 RNA Sequencing 3 05-22-2012 08:54 AM
tophat junctions.bed MerFer Bioinformatics 0 06-16-2010 03:57 AM

Thread Tools
Old 05-06-2014, 09:17 AM   #1
Location: Dresden

Join Date: Jul 2012
Posts: 49
Default Disparate numbers in IGV sashimi plot and junctions.bed

I trying to detect widespread splicing defects and amongst other things I ran tophat to find new splice junctions on my libraries (PE, 100 bases):

tophat -p 12 \
-r 200 --library-type fr-firststrand  \
--transcriptome-index=$bowtie_tindex \
-o $outputdir \
$bowtie_gindex \
$fastq_R1 $fastq_R2
Afterwards, the junction.bed files were merged and used as input for DESeq2 (condition vs control, 2 replicates). Then I tried to look at some of the significant events in IGV to get some nice sushi plots... and here lies my problem. This event, and there are more, has a clear pattern, being expressed in one condition and not in the other:

chr17_24721664_24721852_chr17_24721714_24721802 0 0 264 591
However, when I look at that junction in igv the two libraries that have no expression have a high count:

Is this something to do with the way IGV counts reads in that junction, am I missing something, or is it all wrong?
krespim is offline   Reply With Quote

igvtools, junctions, sashimi, tophat

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 03:45 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO