from sequana import FastQC, sequana_data
fastqc = FastQC(sequana_data("test.fastq"))
for x in 'ACGT':
    fastqc.get_actg_content()[x].hist(
        alpha=0.5, label=x, histtype='step', lw=3, bins=10)
from pylab import legend
legend()