comparison protein_blast_grouping.py @ 4:258afe02e34a draft

planemo upload commit 1098c37fcd6cf35d77748a08aa8b3ee9093c5aea
author cpt
date Thu, 08 Aug 2024 03:52:37 +0000
parents c5e0e05ce58a
children 11701cc63084
comparison
equal deleted inserted replaced
3:c5e0e05ce58a 4:258afe02e34a
42 self, num_hits, sort_key="unique_queries", output_file=sys.stdout 42 self, num_hits, sort_key="unique_queries", output_file=sys.stdout
43 ): 43 ):
44 top_hits = self.get_top_hits(num_hits, sort_key) 44 top_hits = self.get_top_hits(num_hits, sort_key)
45 print(f"# Top {num_hits} Hits", file=output_file) 45 print(f"# Top {num_hits} Hits", file=output_file)
46 print( 46 print(
47 "{:>4} {:<20} {:>10} {:>10}".format( 47 "Rank\tPhage Name\tUnique Query Matches\tUnique Subject Hits",
48 "Rank", "Phage Name", "Unique Query Matches", "Unique Subject Hits"
49 ),
50 file=output_file, 48 file=output_file,
51 ) 49 )
52 for rank, (organism, data) in enumerate(top_hits): 50 for rank, (organism, data) in enumerate(top_hits):
53 print( 51 print(
54 "{:>4} {:<50} {:<25} {:<25}".format( 52 f"{rank}\t{organism}\t{len(data['unique_queries'])}\t{len(data['unique_hits'])}",
55 rank,
56 organism,
57 len(data["unique_queries"]),
58 len(data["unique_hits"]),
59 ),
60 file=output_file, 53 file=output_file,
61 ) 54 )
62 55
63 56
64 def main(): 57 def main():