compound-hets: write summary even if no variants

brentp · Jun 4, 2021 · 5fe122b · 5fe122b
1 parent e6a75f1
commit 5fe122b
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 5 deletions.
diff --git a/CHANGES.md b/CHANGES.md
@@ -1,6 +1,7 @@
-v0.2.3 (dev)
-============
+v0.2.3
+======
 + add `feature_fusion` to impact order list (#92)
++ [compound-hets] write summary even if no variants found
 
 v0.2.2
 ======

diff --git a/src/slivarpkg/comphet.nim b/src/slivarpkg/comphet.nim
@@ -304,11 +304,9 @@ proc main*(dropfirst:bool=false) =
   ovcf.close()
   ivcf.close()
   stderr.write_line &"[slivar compound-hets] wrote {nwritten} variants that were part of a compound het."
-  if ncsqs == 0:
-    stderr.write_line &"[slvar compound-hets] WARNING!!! no variants had any of the requested fields; unable to call compound hets!"
-    quit 0
 
   var summaryPath = getEnv("SLIVAR_SUMMARY_FILE")
+
   if summaryPath == "":
     stderr.write_line counter.tostring(kids)
   else:
@@ -319,6 +317,10 @@ proc main*(dropfirst:bool=false) =
     fh.close()
     stderr.write_line "[slivar] wrote summary table to:" & summaryPath
 
+  if ncsqs == 0:
+    stderr.write_line &"[slvar compound-hets] WARNING!!! no variants had any of the requested fields; unable to call compound hets!"
+    quit 0
+
 when isMainModule:
   import unittest
   suite "compound hets":