"""
Quality histogram a la fastQC
=================================


"""

#####################################
# Get a data set example
from sequana import sequana_data

dataset = sequana_data("test.fastq.gz")


###########################################""
# Create a FastQC instance
from sequana import FastQC

qc = FastQC(dataset, verbose=False)

####################################
# plot the histogram
qc.boxplot_quality()