diff --git a/panaroo/generate_output.py b/panaroo/generate_output.py index 72e502a..c65cbb1 100755 --- 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: