Added print statement in case no genes are sig in gene report.

quadram-institute-bioscience · Dec 13, 2019 · 2ce6166 · 2ce6166
1 parent 9d23aa0
commit 2ce6166
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/albatradis/BlockInsertions.py b/albatradis/BlockInsertions.py
@@ -247,6 +247,7 @@ def gene_statistics(self, forward_plotfile, reverse_plotfile, combined_plotfile,
                 all_genes.append(g)
 
         if len(all_genes) == 0 and len(intergenic_blocks) == 0:
+            print("No significant genes found for chosen parameters.\n")
             return []
 
         self.write_gene_report(all_genes, intergenic_blocks)
@@ -259,6 +260,7 @@ def gene_statistics(self, forward_plotfile, reverse_plotfile, combined_plotfile,
 
     def write_gene_report(self, genes, intergenic_blocks):
         block_filename = os.path.join(self.prefix, "gene_report.csv")
+
         with open(block_filename, 'w') as bf:
             bf.write(str(genes[0].header()) + "\n")
             if not self.use_annotation: