Rescaled qualities.

fastq-demultiplexing/Python/fq_demux.py

@@ -76,8 +76,8 @@ def main():
     bc_start = args.barcode_start - 1
 
     def bc_prob_profile(seq, qual):
-        # return [(letter, 10 ** (-quality / 10))
-        return [(letter, 10 ** (-ord(quality) / 10))
+        # return [(letter, 10 ** ((33 - quality) / 10))
+        return [(letter, 10 ** ((33 - ord(quality)) / 10))
                 for letter, quality
                 in islice(zip(seq, qual), bc_start, bc_start + bc_len)]