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

Read length    Mean length : 1405.68

plot the histogram of the SNRs for each base

plot qc pacbio

Total running time of the script: (0 minutes 0.686 seconds)

Gallery generated by Sphinx-Gallery