Fix bug in bam parsing that misprocessed low qual mutants (resolves #31)

resolves #31
arkal · Nov 7, 2017 · 7b9874c · 7b9874c
1 parent 03eca0f
commit 7b9874c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/transgene.py b/transgene.py
@@ -220,7 +220,7 @@ def get_alignment_info(rna_bam, dna_bam, vcf_record):
             for read in pileup.pileups:
                 if not (read.is_refskip or read.is_del):
                     base = read.alignment.seq[read.query_position]
-                    if ord(read.alignment.qual[read.query_position]) >= 30:
+                    if ord(read.alignment.qual[read.query_position]) >= 63:  # (33PHRED + 30 QUAL)
                         output_counts_[base]['counts'] += 1
                         if read.alignment.is_read2:
                             output_counts_[base]['r'] += 1