# HG changeset patch
# User cpt
# Date 1723088025 0
# Node ID c5e0e05ce58af419e62461b991bef221bf7c9da5
# Parent  8e34cd137d3ec8e9e89041ef544c34d9369359ff
planemo upload commit 120ca2f78836796630f3116ec95ad0326c6e6b6f

diff -r 8e34cd137d3e -r c5e0e05ce58a protein_blast_grouping.py
--- a/protein_blast_grouping.py	Thu Aug 08 03:02:17 2024 +0000
+++ b/protein_blast_grouping.py	Thu Aug 08 03:33:45 2024 +0000
@@ -42,16 +42,20 @@
         self, num_hits, sort_key="unique_queries", output_file=sys.stdout
     ):
         top_hits = self.get_top_hits(num_hits, sort_key)
-        print(f"# Top {num_hits} Hits")
+        print(f"# Top {num_hits} Hits", file=output_file)
         print(
-            "{:<50} {:<25} {:<25}".format(
-                "# Name", "Unique Query Matches", "Unique Subject Hits"
-            )
+            "{:>4} {:<20} {:>10} {:>10}".format(
+                "Rank", "Phage Name", "Unique Query Matches", "Unique Subject Hits"
+            ),
+            file=output_file,
         )
-        for organism, data in top_hits:
+        for rank, (organism, data) in enumerate(top_hits):
             print(
-                "{:<50} {:<25} {:<25}".format(
-                    organism, len(data["unique_queries"]), len(data["unique_hits"])
+                "{:>4} {:<50} {:<25} {:<25}".format(
+                    rank,
+                    organism,
+                    len(data["unique_queries"]),
+                    len(data["unique_hits"]),
                 ),
                 file=output_file,
             )