5.2. read length histograms pacbio data

QC pacbio example

First, let us get a data set example. Note the .bam extension

from sequana import sequana_data

dataset = sequana_data("test_pacbio_subreads.bam")

Create a sequana.pacbio.BAMPacbio instance

from sequana.pacbio import PacbioSubreads

qc = PacbioSubreads(dataset)

plot the histogram of read length

qc.hist_read_length()

0it [00:00, ?it/s]
130it [00:00, 19842.05it/s]

plot the histogram of the SNRs for each base

qc.hist_snr()