Santcroos · Sander Bollen · Sander Bollen
--- a/vtools/gcoverage.py
+++ b/vtools/gcoverage.py
        return [gq]*maxlen


+# Credit:
+# https://gigabaseorgigabyte.wordpress.com/2017/06/26/averaging-basecall-quality-scores-the-right-way/
+# https://git.lumc.nl/klinische-genetica/capture-lumc/vtools/issues/3
+def qualmean(quals):
+    return -10*math.log(sum([10**(q/-10) for q in quals]) / len(quals), 10)