<tool id="rnapaln" name="RNApaln" version="2.1.6.0">
    <description>RNA alignment based on sequence base pairing propensities
    </description>
    <expand macro="requirements" />
    <expand macro="version_command" />
    <expand macro="stdio" />
    <macros>
        <token name="@EXECUTABLE@">RNApaln</token>
        <import>macros.xml</import>
    </macros>
    <command>
<![CDATA[
        RNApaln < $input > $output
        --mode=$compare
        -T$temperature -d$dangling
        $backtrack
        --gapo=$gapo
        --gape=$gape
        --seqw=$seqw
        $endgaps

        #if $varExists('$advancedOptions.nolp')
            $advancedOptions.nooconv
            $advancedOptions.nolp
            $advancedOptions.nogu
            $advancedOptions.noclosinggu
            $advancedOptions.notetra
        #end if

        #if $backtrack and str($compare)=="m"
            ; cat backtrack.file >> $outfile
        #end if
]]>
    </command>

    <inputs>
        <param format="fasta" name="input" type="data" label="Fasta file"/>
        <param name="temperature" size="6" type="float" value="37.0" label="temperature [°C]" help="-T"/>
        <param name="dangling" type="select" label="how to treat dangling end energies" help="-d">
            <option value="0">0: ignore dangling ends</option>
            <option value="1">1: unpaired bases participate in one dangling end only</option>
            <option value="2" selected="True">2: unpaired bases participate in all dangling ends</option>
            <option value="3">3: allow coaxial stacking</option>
        </param>
        <param name="compare" type="select" label="Comparison Option" help="-d">
            <option value="p" selected="True">p: pairwise (1st with 2nd, 3rd with 4th, ...)</option>
            <option value="m">m: matrix (each with each, output in matrix form)</option>
            <option value="f" >f: first (1st with 2nd, 1st with 3rd, ...)</option>
            <option value="c">c: continuous (1st with 2nd, 2nd with 3rd, ...)</option>
        </param>

        <!-- In version 2.1.2 and 2.1.7 the alignment score is wrong if -B is not selected, therefore the default for this option is set to true -->
        <param name="backtrack" type="boolean" checked="true" truevalue="-B" falsevalue="" label="Print an alignment" help="-B"/>

        <param name="gapo" size="8" type="float" value="1.5" label="gap open penalty" help="--gapo"/>
        <param name="gape" size="8" type="float" value="0.67" label="gap extension penalty" help="--gape"/>
        <param name="seqw" size="8" type="float" value="0.5" label="weight of sequence compared to structure" help="--seqw"/>
        <param name="endgaps" type="boolean" checked="false" truevalue="--endgaps" falsevalue="" label="Allow free end-gaps" help="--endgaps"/>

        <conditional name="advancedOptions">
            <param name="advancedSelector" type="select" label="advanced options">
                <option value="basic">basic Options</option>
                <option value="advanced">advanced Options</option>
            </param>
            <when value="advanced">
                <param name="noconv" type="boolean" truevalue="--noconv" falsevalue="" checked="false" label="No conversion of T -> U" help="--noconv"/>
                <param name="nolp" type="boolean" truevalue="--noLP" falsevalue="" checked="false" label="No lonely pairs" help="--noLP  don't allow lonely pairs."/>
                <param name="nogu" type="boolean" truevalue="--noGU" falsevalue="" checked="false" label="No GU pairing" help="--noGU  don't allow pairing of G and U."/>
                <param name="noclosinggu" type="boolean" truevalue="--noClosingGU" falsevalue="" checked="false" label="No GU pairing at the ends" help="--noClosingGU  don't allow pairing of G and U at the ends of helices."/>
                <param name="notetra" type="boolean" truevalue="--noTetra" falsevalue="" checked="false" label="No stabilization for loops, hairpins etc." help="--noTetra"/>
            </when>
        </conditional>
    </inputs>
    <outputs>
        <data format="txt" name="output"/>
        <collection name="matrix_outputs" type="list" label="rna_eps outputs">
            <discover_datasets pattern="(?P&lt;designation&gt;.+)_dp\.ps" ext="rna_eps" visible="true"/>
        </collection>
    </outputs>
    <tests>
    </tests>
    <help>
<![CDATA[

**RNApaln**

RNA alignment based on sequence base pairing propensities.

Uses string-alignment techniques to perform fast pairwise structural alignments of RNAs. Similar to RNApdist secondary structure is incorporated in an approximate manner by computing base pair probabilities, which are then reduced to a vector holding the probability that a base is paired upstream, downstream, or remains unpaired. Such pair propsensity vectors can then be compared using standard alignment algorithms. In contrast to RNApdist, RNApaln performs similarity (instead of distance) alignments, considers both sequence and structure information, and uses affine (rather than linear) gap costs. RNApaln can perform semi-local alignments by using free end gaps, a true local alignment mode is planned.

-----

**Input format**

RNApaln requires one input file
- Fasta file

------

**Outputs**

- output of alignment scores and alignments (if -B option is selected)
- dot plot matrices in postscript format bundled together in a tar file

]]>
    </help>
    <expand macro="requirements" />
</tool>