added check for the number of core genes to prevent a crash if no clu…

…sters passed the core frequency threshold
gtonkinhill · Sep 21, 2023 · 6dc9998 · 6dc9998
1 parent 27044cb
commit 6dc9998
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/panaroo/generate_output.py b/panaroo/generate_output.py
@@ -440,6 +440,10 @@ def generate_core_genome_alignment(
         None
     # Get core nodes
     core_genes = get_core_gene_nodes(G, threshold, num_isolates)
+    if len(core_genes) < 1:
+        print("No gene clusters were present above the core frequency"
+              " threshold! Try adjusting the '--core_threshold' parameter")
+
     core_gene_names = [G.nodes[x]["name"] for x in core_genes]
 
     if codons == True: