def test_basic(): params = {'-g': join(INPATH, 'map_red.tsv'), '-t': join(INPATH, 'go-basic.obo.1'), '-o': join(OUTPATH, 'output.tsv')} response = runner.invoke(_propagate_GO, f"{dict2str(params)}") assert response.exit_code == 0 genes_out = load_df(join(OUTPATH, "output.tsv")) genes_exp = load_df(join(EXPPATH, "new_tree_propagated_GO.tsv")) pdt.assert_frame_equal(genes_out, genes_exp)
def test_basic(): params = { '-i': join(INPATH, 'deepfri_reduced.tsv'), '-g': join(INPATH, 'GO_informative.txt'), '-o': join(OUTPATH, 'output.tsv') } response = runner.invoke(_process_deepfri, f"{dict2str(params)}") assert response.exit_code == 0 genes_out = load_df(join(OUTPATH, "output.tsv")) genes_exp = load_df(join(EXPPATH, "deepfri_out.tsv")) pdt.assert_frame_equal(genes_out, genes_exp)
def test_missing_checkm_and_gtdb(): params = { '-r': join(INPATH, 'cluster_genes/nr.reduced.clstr'), '-g': join(INPATH, 'cluster_genes/sample_genes.fa'), '-c': join(INPATH, 'assemble/merged.min500.contigs.fa'), '-b': join(INPATH, 'metabat2/'), '-t': join(INPATH, 'empty/'), # difference wrt basic '-m': join(INPATH, 'empty/'), # difference wrt basic '-e': join(INPATH, 'eggnog-mapper/eggNOG_reduced.tsv'), '-p': join(OUTPATH), '-o': 'example_split' } response = runner.invoke(_perform_mapping, f"{dict2str(params)}") assert response.exit_code == 0 out = load_df(join(OUTPATH, f"example_split.tsv")) exp = load_df(join(EXPPATH, f"table_missing_checkm_gtdb.tsv")) pdt.assert_frame_equal(out, exp)
def test_missing_checkm_and_gtdb(): params = { '-r': join(INPATH, 'cluster_genes/nr.reduced.clstr'), '-g': join(INPATH, 'cluster_genes/sample_genes.fa'), '-c': join(INPATH, 'assemble/merged.min500.contigs.fa'), '-b': join(INPATH, 'metabat2/'), '-t': join(INPATH, 'empty/'), # difference wrt basic '-m': join(INPATH, 'empty/'), # difference wrt basic '-e': join(INPATH, 'eggnog-mapper/eggNOG_reduced.tsv'), '-o': join(OUTPATH) } response = runner.invoke(_perform_mapping, f"{dict2str(params)}") assert response.exit_code == 0 for name in ['Mapped_genes_cluster', 'Individual_mapped_genes', 'MAGS']: out = load_df(os.path.join(OUTPATH, f"{name}.tsv")) exp = load_df(os.path.join(EXPPATH, f"{name}_missing_checkm_gtdb.tsv")) pdt.assert_frame_equal(out, exp)
def test_basic_split(): params = { '-r': join(INPATH, 'cluster_genes/nr.reduced.clstr'), '-g': join(INPATH, 'cluster_genes/sample_genes.fa'), '-c': join(INPATH, 'assemble/merged.min500.contigs.fa'), '-b': join(INPATH, 'metabat2/'), '-t': join(INPATH, 'gtdbtk/'), '-m': join(INPATH, 'checkm/'), '-e': join(INPATH, 'eggnog-mapper/eggNOG_reduced.tsv'), '-s': '', '-p': join(OUTPATH), '-o': 'example_split' } response = runner.invoke(_perform_mapping, f"{dict2str(params)}") assert response.exit_code == 0 for name in ['mapped_genes_cluster', 'individual_mapped_genes', 'MAGS']: out = load_df(join(OUTPATH, f"example_split_{name}.tsv")) exp = load_df(join(EXPPATH, f"{name}.tsv")) pdt.assert_frame_equal(out, exp)