5.3. BAM module example

Plot histogram of MAPQ values contained in a BAM file

5.3.1. Get a data set (BAM file) for testing

first import the relevant modules

from sequana import BAM, sequana_data

datatest = sequana_data("test.bam", "doc")

Use sequana.bamtools.BAM class to plot the MAPQ historgram

b = BAM(datatest)

# for some unknown reasons, we must call show in this example for this image to
# appear in the doc
from pylab import show

plot bam

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

Gallery generated by Sphinx-Gallery