Fix bug in previous commit

- in EVC, pileup is counting ACGTIDN, where "N" is treated as "N" not "A"
KewinOgink · Dec 10, 2019 · b2fd869 · b2fd869
1 parent 98f1f7e
commit b2fd869
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/dataPrepScripts/ExtractVariantCandidates.py b/dataPrepScripts/ExtractVariantCandidates.py
@@ -287,7 +287,8 @@ def make_candidates(args):
 
                 elif c == "M" or c == "=" or c == "X":
                     for _ in range(advance):
-                        base = BASE2ACGT[SEQ[query_position]]
+                        query_base = SEQ[query_position]
+                        base = BASE2ACGT[query_base] if query_base != "N" else query_base
                         pileup[reference_position][base] += 1
 
                         # those CIGAR operations consumes query and reference