Merge pull request #98 from dufeiyu/bugfix

Fix another bug of division by zero
genome · Oct 13, 2023 · 370dc54 · 370dc54
2 parents 9e5c135 + 1edb108
commit 370dc54
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/dockerfiles/docker-myeloseq/filterHaloplex.py b/dockerfiles/docker-myeloseq/filterHaloplex.py
@@ -503,7 +503,12 @@ def revcomp(seq):
 
             ro = rec.samples[0]['AD'][0]
             ao = rec.samples[0]['AD'][1]
-            rawvaf = round(ao / (ao + ro),4) #rec.samples[0]['AF'][0]
+
+            if ao + ro == 0:
+                rawvaf = 0
+            else:
+                rawvaf = round(ao / (ao + ro),4) #rec.samples[0]['AF'][0]
+
             for v in rec.filter.keys():
                 nrec.filter.add(v)