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


"""

#####################################
# Get a data set example 
from sequana import sequana_data
dataset  = sequana_data("test.fastq", "testing")



###########################################""
# Create a FastQC instance
from sequana import FastQC
qc = FastQC(dataset, verbose=False)

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