Sashimi Plot
Sashimi plots visualize splice junctions from aligned RNA-seq data and a gene annotation track. IGV displays the Sashimi plot in a separate window and allows for more manipulations of the plots than the junctions track.
- To view a Sashimi plot of your alignment data, first zoom out the view to contain the entire region of interest as scrolling and zooming in the Sashimi plot will be limited to this initial region.
- Right click on the alignment track to bring up the pop-up menu, and select Sashimi Plot.
- Select one feature track to serve as the annotation.
- If there is only one possible feature track, e.g., the default RefSeq genes track loaded with the reference genome, then it is automatically loaded.
- If you loaded additional feature tracks, IGV presents a dialog for you to select one for the new plot.
- IGV prompts you to select which alignment tracks you would like to view as Sashimi plots. Select any number and press OK.
The Sashimi plot is displayed in a separate window. The coverage for each alignment track is plotted as a bar graph. Arcs representing splice junctions connect exons. Arcs display the number of reads split across the junction (junction depth). Genomic coordinates and the gene annotation track are shown below the junction tracks.
- Hovering the mouse over each of the exons in the feature annotation track displays additional information in a yellow tooltip.
- Zoom in using the + button at top, and scroll by click-dragging the panel.
- To view only those junctions which overlap a particular exon, select that exon by clicking on it.
- Multiple exons can be selected using ctrl + <click> and they will be highlighted as white boxes.
- To clear selections, click on a blank area of the annotations section of the panel.
Static images of Sashimi plots can also be generated outside IGV with sashimi_plot, a Python tool which is part of the MISO package. Read more about sashimi_plot here.
Popup Menu Options
Command | Description |
Set Exon Coverage Max |
|
Set Junction Coverage Min |
|
Set Junction Coverage Max |
|
Set Color |
|
Show Exon Coverage Data |
|
Text |
|
Combine Strands Forward Strand Reverse Strand | A junction's strandedness is determined by the BAM file XS tag value for the split read. How you assigned the XS tag values to the reads determines whether you potentially display novel junctions or display junctions reflecting previously determined junction annotations. See the Splice Junctions page for more details.
|
Save Image |
|
Source: https://software.broadinstitute.org/software/igv/Sashimi
Posted by: dorothyyoust1959.blogspot.com
Posting Komentar