Mercurial > repos > rnateam > vienna_rna
diff rnafold.xml @ 0:78d673470d45 draft
Imported from capsule None
author | rnateam |
---|---|
date | Wed, 04 Feb 2015 12:05:27 -0500 |
parents | |
children | 5e58cbf27a05 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/rnafold.xml Wed Feb 04 12:05:27 2015 -0500 @@ -0,0 +1,125 @@ +<tool id="rnafold" name="RNAfold" version="2.1.6.0"> + <description>Calculate minimum free energy secondary structures and partition function of RNAs</description> + <expand macro="requirements" /> + <expand macro="version_command" /> + <expand macro="stdio" /> + <macros> + <token name="@EXECUTABLE@">RNAfold</token> + <import>macros.xml</import> + </macros> + <command> + <![CDATA[ + RNAfold + -T $temperature + --dangles=$dangling + #if $measelect.mea == "yes": + --MEA=$measelect.meavalue + #else + $measelect.pf + #end if + #if $varExists('$advancedOptions.nogu'): + $advancedOptions.noconversion + $advancedOptions.gquad + $advancedOptions.nolp + $advancedOptions.nogu + $advancedOptions.noclosinggu + $advancedOptions.canonicalonly + $advancedOptions.circular + #end if + + < $fasta_input + + > $tabular_file +]]> + </command> + + <inputs> + <param format="fasta" name="fasta_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> + <conditional name="measelect"> + <param name="mea" type="select" label="Calculate Maximum Expected accuracy" help="--MEA"> + <option value="no">No</option> + <option value="yes">Yes</option> + </param> + <when value="yes"> + <param name="meavalue" size="6" type="float" value="1.0" label="Gamma Value"/> + </when> + <when value="no"> + <param name="pf" type="boolean" checked="false" truevalue="--partfunc" falsevalue="" label="Calculate Partition Function" help="--partfunc"/> + </when> + </conditional> + <conditional name="advancedOptions"> + <param name="advancedSelector" type="boolean" checked="false" label=" advanced options"/> + <when value="true"> + <param name="noconversion" type="boolean" truevalue="--noconv" falsevalue="" checked="false" label="no conversion" help="--noconv do not convert thymine to uracile (T -> U)."/> + <param name="gquad" type="boolean" truevalue="--gquad" falsevalue="" checked="false" label="G Quadruplex formation" help="-g take into account G Quadruplex formation"/> + <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"/> + <param name="canonicalonly" type="boolean" truevalue="--canonicalBPonly" falsevalue="" checked="false" label="Canonical basepairing only" help="--canonicalBPonly"/> + <param name="circular" type="boolean" truevalue="--circ" falsevalue="" checked="false" label="Assume circular RNA structure" help="--circ"/> + </when> + </conditional> + </inputs> + <outputs> + <data format="tabular" name="tabular_file"/> + <collection name="sequence_outputs" type="list" label="rna_eps outputs"> + <discover_datasets pattern="(?P<designation>.+)_ss\.ps" ext="eps" /> + </collection> + <collection name="matrix_outputs" type="list" label="rna_eps outputs"> + <filter>measelect['pf'] is True</filter> + <discover_datasets pattern="(?P<designation>.+)_dp\.ps" ext="rna_eps" visible="true"/> + </collection> + + </outputs> + <tests> + </tests> + <help> +<![CDATA[ +**RNAfold** + +The program reads RNA sequences, calculates their minimum free +energy (mfe) structure and the mfe structure in dot-bracket notation. + +If the -p option was given it also computes the +partition function (pf) and base pairing probability matrix. + +The dot plot of the base pairing probability matrix shows a matrix of squares with area proportional to the pairing +probability in the upper right half, and one square for each pair in the +minimum free energy structure in the lower left half. For each pair i-j with +probability p>10E-6 there is a line of the form + +i j sqrt(p) ubox + +in the PostScript file, so that the pair probabilities can be easily extracted. + +The sequences have to be provided in FASTA format. The first word (max. 42 char) of the FASTA header will be used for output file names. PostScript files "name_ss.ps" and "name_dp.ps" are produced for the structure and dot plot, respectively. +The program will read the whole FASTA input file and provide output for each found sequence. + + +----- + +**Input format** + +RNAfold requires one input file +- FASTA file + +------ + +**Outputs** + +- Secondary structures in dot-bracket notation + +- several possible postscript images bundled together in a tar file + - secondary structure for each sequence in the input file + - if partition function is calculated (--MEA or --partfunc is set) then also the pairing probabilty matrix is generated for each sequence + +]]> + </help> + <expand macro="citations" /> +</tool>