from pylab import hist
from sequana import sequana_data, BAM
b = BAM(sequana_data("test.bam"))
hist(b.get_mapped_read_length())