def test_blast_lca(taxa_db_simple, simple_m8): test_path = join(util.file.get_test_input_path(), 'TestTaxonomy') expected = textwrap.dedent("""\ C\tM04004:13:000000000-AGV3H:1:1101:12068:2105\t2 C\tM04004:13:000000000-AGV3H:1:1101:13451:2146\t2 C\tM04004:13:000000000-AGV3H:1:1101:13509:2113\t2 C\tM04004:13:000000000-AGV3H:1:1101:14644:2160\t2 C\tM04004:13:000000000-AGV3H:1:1101:18179:2130\t2 C\tM04004:13:000000000-AGV3H:1:1111:10629:2610\t2 C\tM04004:13:000000000-AGV3H:1:1111:10629:26101\t2 """) out = StringIO() with simple_m8 as f: metagenomics.blast_lca(taxa_db_simple, f, out, paired=True) out.seek(0) assert out.read() == expected
def test_blast_lca(taxa_db_simple, simple_m8): test_path = join(util.file.get_test_input_path(), 'TestTaxonomy') expected = textwrap.dedent("""\ M04004:13:000000000-AGV3H:1:1101:12068:2105\t2 M04004:13:000000000-AGV3H:1:1101:13451:2146\t2 M04004:13:000000000-AGV3H:1:1101:13509:2113\t2 M04004:13:000000000-AGV3H:1:1101:14644:2160\t2 M04004:13:000000000-AGV3H:1:1101:18179:2130\t2 M04004:13:000000000-AGV3H:1:1111:10629:2610\t2 M04004:13:000000000-AGV3H:1:1111:10629:26101\t2 """) out = StringIO() with simple_m8 as f: metagenomics.blast_lca(taxa_db_simple, f, out, paired=True) out.seek(0) assert out.read() == expected