<?xml version="1.0"?>
<tool id="hyphy_meme" name="HyPhy-MEME" version="@VERSION@+galaxy0">
    <description>: Mixed Effects Model of Evolution</description>
    <macros>
        <import>macros.xml</import>
    </macros>
    <expand macro="requirements" />
    <expand macro="version_command" />
    <command detect_errors="exit_code">
    <![CDATA[
    export MPI_NODE_COUNT=\${GALAXY_SLOTS:-1} &&
    ln -s '$input_file' meme_input.fa &&
    ln -s '$input_nhx' meme_input.nhx &&
    echo $gencodeid > tool_params &&
    echo `pwd`/meme_input.fa >> tool_params &&
    echo `pwd`/meme_input.nhx >> tool_params &&
    echo $branches >> tool_params &&
    echo $p_value >> tool_params &&
    @HYPHY_INVOCATION@ \$HYPHY_LIB/TemplateBatchFiles/SelectionAnalyses/MEME.bf > '$meme_log'
    ]]>
    </command>
    <inputs>
        <expand macro="inputs" />
        <expand macro="gencode" />
        <expand macro="branches" />
        <param name="p_value" type="float" value=".1" min="0" max="1" label="P-value" />
    </inputs>
    <outputs>
        <data name="meme_log" format="txt" />
        <data name="meme_output" format="hyphy_json_results" from_work_dir="meme_input.fa.MEME.json" />
    </outputs>
    <tests>
        <test>
            <param name="input_file" ftype="fasta" value="meme-in1.fa" />
            <param name="input_nhx" ftype="nhx" value="meme-in1.nhx" />
            <param name="branches" value="1" />
            <param name="p_value" value="0.1" />
            <output name="meme_output" file="meme-out1.json" compare="sim_size" />
        </test>
    </tests>
    <help>
        <![CDATA[
MEME employs a mixed-effects maximum likelihood approach to test the hypothesis that individual sites have been subject to episodic positive or diversifying selection. In other words,  MEME aims to detect sites evolving under positive selection under a proportion of branches.
        ]]>
    </help>
    <expand macro="citations">
        <citation type="doi">10.1371/journal.pgen.1002764</citation>
    </expand>
</tool>