def test_parse_hmm_header(self):
gene_name = "gspD"
c_gene = CoreGene(self.model_location, gene_name, self.profile_factory)
report_path = os.path.join(self.cfg.working_dir(),
gene_name + self.cfg.res_search_suffix())
report = GembaseHMMReport(c_gene, report_path, self.cfg)
hmm_hit = [
">> NC_xxxxx_xx_056141 C ATG TAA 6260390 6261757 Valid PA5567 1368 _NP_254254.1_ PA5567 1 6260390 6261757 | tRNA modific"
]
hit_id = report._parse_hmm_header(hmm_hit)
self.assertEqual(hit_id, 'NC_xxxxx_xx_056141')
