0
|
1 <tool id="rnalfold" name="RNALfold" version="2.1.6.0">
|
|
2 <description>calculates locally stable secondary structures of RNA</description>
|
|
3 <expand macro="requirements" />
|
|
4 <expand macro="version_command" />
|
|
5 <expand macro="stdio" />
|
|
6 <macros>
|
|
7 <token name="@EXECUTABLE@">RNALfold</token>
|
|
8 <import>macros.xml</import>
|
|
9 </macros>
|
|
10 <command>
|
|
11 <![CDATA[
|
|
12 RNALfold -L $span -T $temperature -d$dangling < $fasta_input | sed -r 's/(^[.()]*) +(\(.*\)) +(.*)/\1\t\2\t\3/' > $out_file
|
|
13 #if $varExists('$advancedOptions.noconversion')
|
|
14 #if $advancedOptions.noconversion
|
|
15 --noconv
|
|
16 #end if
|
|
17 #if $advancedOptions.gquad
|
|
18 --gquad
|
|
19 #end if
|
|
20 #if $advancedOptions.nolp
|
|
21 --noLP
|
|
22 #end if
|
|
23 #if $advancedOptions.nogu
|
|
24 --noGU
|
|
25 #end if
|
|
26 #if $advancedOptions.noclosinggu
|
|
27 --noClosingGU
|
|
28 #end if
|
|
29 #if $advancedOptions.notetra
|
|
30 --noTetra
|
|
31 #end if
|
|
32 #end if
|
|
33 ]]>
|
|
34 </command>
|
|
35 <inputs>
|
|
36 <param format="fasta" name="fasta_input" type="data" label="FASTA file"/>
|
|
37 <param name="span" type="integer" min="0" max="5000" value="150" label="base pair span" help="Maximal distance between two paired bases. (-L)"/>
|
|
38 <param name="temperature" size="6" type="float" value="37.0" label="Temperature [°C]" help="-T"/>
|
|
39 <param name="dangling" type="select" label="how to treat dangling end energies" help="-d">
|
|
40 <option value="2" selected="true">unpaired bases participate in all dangling ends (2)</option>
|
|
41 <option value="0">ignore dangling ends (0)</option>
|
|
42 <option value="1">unpaired bases participate in one dangling end only (1)</option>
|
|
43 <option value="3">allow coaxial stacking (3)</option>
|
|
44 </param>
|
|
45 <conditional name="advancedOptions">
|
|
46 <param name="advancedSelector" type="select" label="advanced options">
|
|
47 <option value="basic">basic Options</option>
|
|
48 <option value="advanced">advanced Options</option>
|
|
49 </param>
|
|
50 <when value="advanced">
|
|
51 <param name="noconversion" type="boolean" checked="false" label="no conversion" help="--noconv do not convert thymine to uracile (T -> U)."/>
|
|
52 <param name="gquad" type="boolean" checked="false" label="G Quadruplex formation" help="-g take into account G Quadruplex formation"/>
|
|
53 <param name="nolp" type="boolean" checked="false" label="No lonely pairs" help="--noLP don't allow lonely pairs."/>
|
|
54 <param name="nogu" type="boolean" checked="false" label="No GU pairing" help="--noGU don't allow pairing of G and U."/>
|
|
55 <param name="noclosinggu" type="boolean" checked="false" label="No GU pairing at the ends" help="--noClosingGU don't allow pairing of G and U at the ends of helices."/>
|
|
56 <param name="notetra" type="boolean" checked="false" label="No stabilization for loops, hairpins etc." help="--noTetra"/>
|
|
57 </when>
|
|
58 </conditional>
|
|
59 </inputs>
|
|
60 <outputs>
|
|
61 <data format="txt" name="out_file"/>
|
|
62 </outputs>
|
|
63 <tests>
|
|
64 </tests>
|
|
65 <help>
|
|
66 <![CDATA[
|
|
67 **RNALfold**
|
|
68
|
|
69 Compute locally stable RNA secondary structure with a maximal base pair span.
|
|
70 For a sequence of length n and a base pair span of L the algorithm uses only
|
|
71 O(n+L*L) memory and O(n*L*L) CPU time. *Thus it is practical to "scan" very
|
|
72 large genomes for short RNA structures*.
|
|
73 Output consists of a list of secondary structure components of size <= L, one
|
|
74 entry per line. Each output line contains the predicted local structure its
|
|
75 energy in kcal/mol and the starting position of the local structure.
|
|
76
|
|
77 -----
|
|
78
|
|
79 **Input format**
|
|
80
|
|
81 - RNALfold requires one input file in FASTA format
|
|
82
|
|
83 ------
|
|
84
|
|
85 **Outputs**
|
|
86
|
|
87 - text output with dot-bracket notation and free energies of the secondary structures
|
|
88
|
|
89 ]]>
|
|
90 </help>
|
|
91 <expand macro="requirements" />
|
|
92 </tool>
|