<?xml version="1.0"?>
<tool id="hyphy_gard" name="HyPhy-GARD" version="@VERSION@+galaxy0">
  <description>: Recombination Detection</description>
  <macros>
    <import>macros.xml</import>
  </macros>
  <requirements>
    <requirement type="package" version="2.0.4">translate-gard</requirement>
    <requirement type="package" version="2.3.11">hyphy</requirement>
  </requirements>
  <version_command>HYPHYMP --version | tail -n 1 | awk '{print $1}'</version_command>
  <command detect_errors="exit_code"><![CDATA[
    export MPI_NODE_COUNT=\${GALAXY_SLOTS:-2} &&
    ln -s '$input_file' gard_input.fa &&
    echo `pwd`/gard_input.fa > tool_params &&
    echo "010010" >> tool_params &&
    echo $rate_cond.rate >> tool_params &&
    #if str($rate_cond.rate) != "1":
        echo $rate_cond.rate_classes >> tool_params &&
    #end if
    echo `pwd`/"gard_input.GARD" >> tool_params &&
    @HYPHYMPI_INVOCATION@ \$HYPHYMPI_LIB/TemplateBatchFiles/GARD.bf > '$gard_log' &&
    translate-gard -i gard_input.GARD -j gard_input.GARD.json -o gard_input.fa.translated.json
    ]]></command>
  <inputs>
    <param name="input_file" type="data" format="fasta" label="Input FASTA file"/>
    <conditional name="rate_cond">
      <param name="rate" type="select" label="Rate variation">
        <option value="1">None</option>
        <option value="2">General Discrete</option>
        <option value="3">Beta-Gamma</option>
      </param>
      <when value="1"/>
      <when value="2">
        <param name="rate_classes" type="integer" value="2" min="2" max="6" label="Rate classes"/>
      </when>
      <when value="3">
        <param name="rate_classes" type="integer" value="2" min="2" max="6" label="Rate classes"/>
      </when>
    </conditional>
  </inputs>
  <outputs>
    <data name="gard_log" format="txt"/>
    <data name="gard_output" format="nex" from_work_dir="gard_input.GARD_finalout"/>
    <data name="translated" format="hyphy_json_results" from_work_dir="gard_input.fa.translated.json"/>
  </outputs>
  <tests>
    <test>
      <param name="input_file" ftype="fasta" value="gard-in1.fa"/>
      <output name="gard_output" file="gard-out1.nex" compare="sim_size" delta="160000"/>
    </test>
  </tests>
  <help><![CDATA[
GARD (Genetic Algorithm for Recombination Detection) is a method to screen a multiple sequence analysis for the presence of recombination and is extremely useful as a pre-processing step for selection inference. Because recombinant sequences cannot be adequately described with a single phylogenetic history, selection inference on recombinant data often leads to a significant increase in false positives. GARD alleviates this concern by comprehensively screening an alignment for recombination breakpoints and inferring a unique phylogenetic history for each detected recombination block.
        ]]></help>
  <expand macro="citations">
    <citation type="doi">10.1093/molbev/msl051</citation>
  </expand>
</tool>