from sequana import sequana_data, BAM
b = BAM(sequana_data("test.bam"))
b.plot_read_length()