Strand-specific coverage track for a BAM file#

Create a StrandSpecificCoverageTrack to show the strand-specific coverage information in a BAM file.

chr22
[<Axes: title={'center': 'chr22'}, ylabel='RNA-seq'>]

from pygv.viewer import GenomeViewer
from pygv.tracks.bam_track import StrandSpecificCoverageTrack
gv = GenomeViewer()
track = StrandSpecificCoverageTrack(
    "../examples/data/s03.chr22.bam",
    name="RNA-seq", y_label_rotation="vertical", y_label_ha="center")
gv.add_track(track)
gv.plot("chr22", 40345142, 40367239)

Total running time of the script: (0 minutes 0.115 seconds)

Gallery generated by Sphinx-Gallery