Fix a bug in --strictly-novel

When only one parent has the alternate read, the score was set to -inf
samtools · Sep 4, 2024 · 418e7fc · 418e7fc
1 parent bbe4a1e
commit 418e7fc
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/plugins/trio-dnm2.c b/plugins/trio-dnm2.c
@@ -974,7 +974,9 @@ static double process_trio_ACM(args_t *args, priors_t *priors, int nals, double
         int ial = *al1;
         if ( qs[iMOTHER][ial] + qs[iFATHER][ial] != 0 )
         {
-            double tmp = subtract_log(0,qs[iMOTHER][ial]) + subtract_log(0,qs[iFATHER][ial]);
+            double tmp = 0;
+            if ( qs[iMOTHER][ial] ) tmp += subtract_log(0,qs[iMOTHER][ial]);
+            if ( qs[iFATHER][ial] ) tmp += subtract_log(0,qs[iFATHER][ial]);
             sum = sum_log(sum,tmp);
             max += tmp;
 #if DEBUG