Quit rel.abun. calculation if there are no strains

broadinstitute · Apr 13, 2021 · 3435c15 · 3435c15
1 parent c0a6476
commit 3435c15
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/strainge/search_tool.py b/src/strainge/search_tool.py
@@ -304,6 +304,9 @@ def calc_relative_abundance(self, sample, result):
         Applies a regression model to estimate the relative abundances of reported strains.
         """
 
+        if not result.strains:
+            return
+
         # Reset cache because the StrainGST algorithm may have removed k-mers from the reference k-mer sets
         self.pangenome.reset_cache()
         if self.top == 1: