def METHOD_NAME(run_component):
    """
    Test a simple merge with fully overlapping observations
    """
    run_component([
        "--input", input_sample1_file,
        "--input", input_sample2_file,
        "--output", "merge.h5mu",
        "--output_compression", "gzip"])
    assert Path("merge.h5mu").is_file()
    concatenated_data = md.read("merge.h5mu")
    data_sample1 = md.read(input_sample1_file)
    data_sample2 = md.read(input_sample2_file)
    assert concatenated_data.n_mod == 2
    assert concatenated_data.obsm_keys() == ['prot', 'rna']
    # In this case, the observations overlap perfectly in the two input files
    assert concatenated_data.n_obs == data_sample1.n_obs
    assert concatenated_data.n_obs == data_sample2.n_obs
    pd.testing.assert_index_equal(concatenated_data.obs.index, data_sample1.obs.index)
    pd.testing.assert_index_equal(concatenated_data.obs.index, data_sample2.obs.index)
    assert set(data_sample1.var_keys()) | set(data_sample2.var_keys()) == \
           set(concatenated_data.var_keys())
    assert set(concatenated_data.var_names) ==  (set(data_sample1.var_names) | set(data_sample2.var_names))
    assert concatenated_data.var_keys() == ['gene_id', 'feature_type', 'genome']