def test_best_hit(self):
gene_name = 'gspD'
c_gene = CoreGene(self.model_location, gene_name, self.profile_factory)
shutil.copy(
self.find_data("hmm", gene_name + self.cfg.res_search_suffix()),
self.cfg.working_dir())
report_path = os.path.join(self.cfg.working_dir(),
gene_name + self.cfg.res_search_suffix())
report = GembaseHMMReport(c_gene, report_path, self.cfg)
self.assertIsNone(report.best_hit())
report.extract()
best_hit = report.best_hit()
hit_expected = CoreHit(c_gene, "NC_xxxxx_xx_056141", 803,
"NC_xxxxx_xx", 141, float(2e-236), float(779.2),
float(1.000000), (741.0 - 104.0 + 1) / 803, 104,
741)
self.assertEqual(hit_expected, best_hit)
