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)
b.plot_bar_mapq()


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

show()