from sequana import *
from pylab import linspace, plot, grid, xlabel, ylabel

b = BAM(sequana_data("measles.fa.sorted.bam"))
b.plot_insert_size()