Extracting read coverage data from a Gviz AlignmentsTrack
1
0
Entering edit mode
@willbradshaw-19193
Last seen 2.1 years ago

I'm using Gviz to make read-coverage and Sashimi plots for a paper, and I'd like to extract information about how many reads contribute to each splice junction in the Sashimi plot. These numbers must be inferred at some point since Gviz is generating the AlignmentsTrack object directly from the BAM file, but it's not clear to me whether it's possible for a user to get access to the intermediate steps used to generate the Sashimi plot. Is there a way to do this?

gviz sashimi • 256 views
ADD COMMENT
0
Entering edit mode
Robert Ivanek ▴ 700
@robert-ivanek-5892
Last seen 16 hours ago
Switzerland

Hi, You can simply read in the part of your BAM file with GAlignments function and afterwards extract junction counts with summarizeJunctions, both function come from GenomicAlignments package.

ADD COMMENT

Login before adding your answer.

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