turns out that the boost second moment is not the second moment that is the variance. We now correctly calculate stddev for quality scores and k-mer variance, and we can plot them too!