Gviz sashimi read count
1
0
Entering edit mode
bernt.popp • 0
@berntpopp-11398
Last seen 7.7 years ago

Hey,

I would like to plot the count of spliced reads supporting the junction above the splice bezier-curve, similar to the original sashimi plot now integrated in IGV.

Is there an option in Gviz for this?

 

Best, Bernt

Gviz sashimi sashimi plots read counting bam • 2.3k views
ADD COMMENT
1
Entering edit mode
Robert Ivanek ▴ 750
@robert-ivanek-5892
Last seen 12 months ago
Switzerland

Hi Bernt,

You can set the parameter sashimiNumbers  to TRUE. Then it should work. Either  when you create the track or when using plotTracks function.

Best

Robert

 

ADD COMMENT
0
Entering edit mode

--Hi,

is it possible to store in a file the count of spliced reads displayed by the sashimi plot when i set the parameter sashimiNumbers  to TRUE ?

best,

 

laurent --

 

 

ADD REPLY
0
Entering edit mode

To get those numbers via Gviz seems to me too complicated. You can use package GenomicAlignments to get those numbers from BAM files directly: function readGAlignments or readGAlignmentPairs to read in the data from BAM file and summarizeJunctions to produce the counts for junctions.

ADD REPLY
0
Entering edit mode

ok thank you so much for this help.

 

ADD REPLY

Login before adding your answer.

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