| Next changeset 1:671b25f012e9 (2016-12-21) |
|
Commit message:
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/vienna_rna commit 0065dafe7bbd382bb995b28cc4089c9e4f4eeeb9 |
|
added:
macros.xml rnaduplex.xml test-data/kinfold_input.txt test-data/rna2dfold_input1.txt test-data/rna2dfold_result1.txt test-data/rnaaliduplex_input1.clustal test-data/rnaaliduplex_result1.txt test-data/rnaalifold_input1.clustal test-data/rnaalifold_result1.txt test-data/rnacofold_input1.fas test-data/rnacofold_result1.txt test-data/rnadistance_input1.dbn test-data/rnadistance_result1.txt test-data/rnaduplex_input1.fa test-data/rnaduplex_result1.txt test-data/rnaeval_input1.dbn test-data/rnaeval_result1.txt test-data/rnafold_input1.fa test-data/rnafold_input2.fa test-data/rnafold_result1.txt test-data/rnafold_result2.txt test-data/rnafold_result3.txt test-data/rnaheat_input1.fa test-data/rnaheat_result1.txt test-data/rnainverse_input1.clu test-data/rnalalifold_input1.clustal test-data/rnalalifold_result1.txt test-data/rnalfold_input1.fa test-data/rnalfold_result1.txt test-data/rnapaln_input1.fa test-data/rnapaln_input1.fas test-data/rnapaln_result1.txt test-data/rnapdist_input1.fa test-data/rnapdist_result1.txt test-data/rnapkplex_input1.fa test-data/rnapkplex_result1.txt test-data/rnaplex_input1.fa test-data/rnaplex_result1.txt test-data/rnaplfold_input1.fa test-data/rnaplfold_result1.ps test-data/rnaplfold_result2.ps test-data/rnaplot_input1.dbn test-data/rnaplot_input1.fa test-data/rnasnoop_input1a.fa test-data/rnasnoop_input1b.fa test-data/rnasnoop_result1.txt test-data/rnasubopt_input1.fa test-data/rnasubopt_result1.txt test-data/rnaup_input1.fa test-data/rnaup_result1.txt vienna_rna.tar.gz |
| b |
| diff -r 000000000000 -r 204093ddb00a macros.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/macros.xml Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,26 @@ +<macros> + <xml name="requirements"> + <requirements> + <requirement version="2.2.10">viennarna</requirement> + </requirements> + </xml> + <token name="@VERSION@">2.2.10</token> + <xml name="version_command"> + <version_command>@EXECUTABLE@ --version</version_command> + </xml> + + <xml name="stdio"> + <stdio> + <exit_code range="1:" /> + <exit_code range=":-1" /> + <regex match="Error:" /> + <regex match="Exception:" /> + </stdio> + </xml> + + <xml name="citations"> + <citations> + <citation type="doi">10.1186/1748-7188-6-26</citation> + </citations> + </xml> +</macros> |
| b |
| diff -r 000000000000 -r 204093ddb00a rnaduplex.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/rnaduplex.xml Tue Dec 06 12:33:48 2016 -0500 |
| [ |
| @@ -0,0 +1,85 @@ +<tool id="viennarna_rnaduplex" name="@EXECUTABLE@" version="@VERSION@.0"> + <description>Compute the structure upon hybridization of two RNA strands</description> + <macros> + <token name="@EXECUTABLE@">RNAduplex</token> + <import>macros.xml</import> + </macros> + <expand macro="requirements" /> + <expand macro="stdio" /> + <expand macro="version_command" /> + <command> +<![CDATA[ + RNAduplex < '$input' > '$tabularFile' + -T$temperature -d$dangling --deltaEnergy=$denergy + #if $varExists('$advancedOptions.noconversion') + $advancedOptions.noconversion + $advancedOptions.noLP + $advancedOptions.noGU + $advancedOptions.noClosingGU + $advancedOptions.noTetra + #if $advancedOptions.nsp + --nsp='$advancedOptions.nsp' + #end if + #end if +]]> + </command> + <inputs> + <param format="fasta" name="input" type="data" label="Fasta file"/> + <param name="temperature" 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 argument="--denergy" type="float" value="0.0" label="delta energy range [kcal/mol] for suboptimal structures" help=""/> + <param argument="--sorted" type="boolean" truevalue="--sorted" falsevalue="" checked="false" label="Sort by energy" help=""/> + <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 argument="--noconversion" type="boolean" truevalue="" falsevalue="--noconv" checked="true" label="Convert Thymine to Uracil (T -> U)" help="Avoids confusion with DNA sequences"/> + <param argument="--noLP" type="boolean" truevalue="" falsevalue="--noLP" checked="true" label="Allow lonely base-pairs" help=""/> + <param argument="--noGU" type="boolean" truevalue="" falsevalue="--noGU" checked="true" label="Allow GU pairing" help="Allow pairing of G and U at the ends of helices."/> + <param argument="--noClosingGU" type="boolean" truevalue="" falsevalue="--noClosingGU" checked="true" label="Allow GU pairing at the ends" help="Allow pairing of G and U at the ends of helices."/> + <param argument="--noTetra" type="boolean" truevalue="" falsevalue="--noTetra" checked="true" label="Allow stabilization for loops, hairpins etc." help=" Include special tabulated stabilizing energies for tri-, tetra- and hexaloop hairpins. Mostly for testing."/> + <param argument="--nsp" type="text" value="" label="Allow other pairs in addition to the usual AU,GC,and GU pairs." + help="Its argument is a comma separated list of additionally allowed pairs. If the first character is '-' then AB will imply that AB and BA are allowed pairs. e.g. RNAfold -nsp -GA will allow GA and AG pairs. Nonstandard pairs are given 0 stacking energy."/> + </when> + <when value="basic"> + </when> + </conditional> + </inputs> + <outputs> + <data format="txt" name="tabularFile"/> + </outputs> + <tests> + <test> + <param name="input" value="rnaduplex_input1.fa"/> + <output name="out_file" file="rnaduplex_result1.txt"/> + </test> + </tests> + <help> +<![CDATA[ +**RNAduplex** +RNA duplex reads two RNA sequences and computes optimal and suboptimal secondary structures for their hybridization. The calculation is simplified by allowing only inter−molecular base pairs, for the general case use RNAcofold. The computed optimal and suboptimal structure are written to stdout, one structure per line. Suboptimal structures are computed with energy in a certain range of the optimum (kcal/mol). Default is calculation of mfe structure only. + +----- + +**Input format** + +RNAduplex requires one input file in Fasta format. The secondary strucures will be calculated pairwise: 1st sequence with 2nd sequence, 3rd sequence with 4th sequence, etc. + +------ + +**Outputs** + +Each line consist of: The structure in dot bracket format with a "&" separating the two strands. The range of the structure in the two sequences in the format "from,to : from,to"; the energy of duplex structure in kcal/mol. The format is especially useful for computing the hybrid structure between a small probe sequence and a long target sequence. + + +]]> + </help> + <expand macro="citations" /> +</tool> |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/kinfold_input.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/kinfold_input.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,1 @@ +ACUGAUCGUAGUCAC |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rna2dfold_input1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rna2dfold_input1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. +.((((((..(((..........))).(((((.......))))).....(((((.......))))))))))).. |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rna2dfold_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rna2dfold_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,6 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00) +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00) <ref 1> +.((((((..(((..........))).(((((.......))))).....(((((.......))))))))))).. (-17.30) <ref 2> +k l en structure |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaaliduplex_input1.clustal --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaaliduplex_input1.clustal Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,7 @@ +CLUSTAL 2.1 multiple sequence alignment + + +Anolis_carolinensis_chrUn_GL34 TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +Anolis_carolinensis_chrUn_GL35 GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA + ************************************** ************** ******* ********* + |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaaliduplex_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaaliduplex_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,1 @@ +(((((((..(.......((((..(..((((((..((((((.......((((((..(..(((((..(((((((.&)))))))..).......))))..)..))))))..)))))).......))))))..)..)))))..))))))). 1,73 : 1,73 (-40.30) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaalifold_input1.clustal --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaalifold_input1.clustal Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,7 @@ +CLUSTAL 2.1 multiple sequence alignment + + +Anolis_carolinensis_chrUn_GL34 TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +Anolis_carolinensis_chrUn_GL35 GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA + ************************************** ************** ******* ********* + |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaalifold_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaalifold_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUCGAUGCCCACAUUCUCCA +(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). (-26.45 = -26.20 + -0.25) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnacofold_input1.fas --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnacofold_input1.fas Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA&UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnacofold_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnacofold_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,3 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-A +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA&UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((((((((.(((((((((((..&.))))))..((((........)))).(((((.......))))).....))))))))))).))))))))))).. (-55.90) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnadistance_input1.dbn --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnadistance_input1.dbn Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. +(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnadistance_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnadistance_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,1 @@ +f: 4 |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaduplex_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaduplex_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaduplex_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaduplex_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,3 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC +.((((((..((......((((...((((..(((.((((((.......((((((..(..((((...(((((((.&)))))))..........))))..)..))))))..))))))..)))..)))).......)))).)))))))). 1,73 : 1,72 (-40.70) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaeval_input1.dbn --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaeval_input1.dbn Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,3 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaeval_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaeval_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,3 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnafold_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnafold_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnafold_input2.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnafold_input2.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnafold_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnafold_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,6 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-A +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00) +>Anolis_carolinensis_chrUn_GL343207.trna3-A +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA +(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). (-29.60) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnafold_result2.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnafold_result2.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,6 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-A +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +..(((.....((((....((..(((....)))..))...)))).....(((((.......)))))....))). ( -3.37) +>Anolis_carolinensis_chrUn_GL343207.trna3-A +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA +(((((..((.(((.....((..(((....)))..))......))).))(((((.......)))))..))))). ( -5.02) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnafold_result3.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnafold_result3.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,3 @@ +>Sequence +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaheat_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaheat_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +> comment 1 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaheat_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaheat_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,102 @@ +> comment 1 +0 0.0407267 +1 0.0435052 +2 0.0470227 +3 0.0505646 +4 0.0553391 +5 0.0603894 +6 0.0662047 +7 0.0718158 +8 0.0782007 +9 0.0863514 +10 0.0950516 +11 0.10505 +12 0.115613 +13 0.128245 +14 0.141712 +15 0.157281 +16 0.17497 +17 0.195056 +18 0.216799 +19 0.241743 +20 0.270432 +21 0.302132 +22 0.337644 +23 0.377263 +24 0.421546 +25 0.470799 +26 0.525589 +27 0.587016 +28 0.653829 +29 0.726606 +30 0.806458 +31 0.892656 +32 0.984857 +33 1.08272 +34 1.18601 +35 1.29332 +36 1.40305 +37 1.51467 +38 1.62618 +39 1.73524 +40 1.84032 +41 1.9407 +42 2.03385 +43 2.11859 +44 2.19443 +45 2.26182 +46 2.31999 +47 2.36969 +48 2.41182 +49 2.44826 +50 2.47967 +51 2.50769 +52 2.534 +53 2.55942 +54 2.5851 +55 2.61217 +56 2.64151 +57 2.67329 +58 2.70759 +59 2.74465 +60 2.78458 +61 2.8279 +62 2.8742 +63 2.92425 +64 2.97698 +65 3.03283 +66 3.09163 +67 3.15591 +68 3.23033 +69 3.3067 +70 3.38122 +71 3.45791 +72 3.54606 +73 3.63843 +74 3.73392 +75 3.83847 +76 3.94646 +77 4.05462 +78 4.16571 +79 4.28366 +80 4.39944 +81 4.50863 +82 4.61035 +83 4.70186 +84 4.78345 +85 4.85204 +86 4.90509 +87 4.9392 +88 4.95425 +89 4.95214 +90 4.93688 +91 4.91166 +92 4.87691 +93 4.8314 +94 4.78561 +95 4.75271 +96 4.71866 +97 4.66056 +98 4.57671 +99 4.48395 +100 4.39006 |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnainverse_input1.clu --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnainverse_input1.clu Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). +gggccccnnaauunnnnnnnnaauunGGGGcnnnnnnnAAAAAnnnnnuuuuannnnnnnuaaaaggggcccn |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnalalifold_input1.clustal --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnalalifold_input1.clustal Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,7 @@ +CLUSTAL 2.1 multiple sequence alignment + + +Anolis_carolinensis_chrUn_GL34 TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +Anolis_carolinensis_chrUn_GL35 GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA + ************************************** ************** ******* ********* + |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnalalifold_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnalalifold_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +(((.(((((.(((((.......)))))..)).(((((.......)))))..)))))) (-19.05) 17 - 73 +((((........)))). ( -5.10) 10 - 26 +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUCGAUGCCCACAUUCUCCA + |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnalfold_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnalfold_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnalfold_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnalfold_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,28 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +.((((....)))). ( -0.10) 56 +.(((..((((....)))).))). ( -3.40) 51 +.(((((.......))))). ( -7.70) 48 +.((((..(((((.......)))))..)))). (-10.30) 42 +.(((((...(((((.......))))).))))). (-11.50) 40 +.((.(((((...(((((.......))))).))))))) (-12.10) 37 +.(((((.......))))). ( -5.80) 26 +.((.(((((.......)))))..)). ( -6.70) 23 +.((((....)))). ( -3.20) 21 +.((.((((....)))).)). ( -4.70) 18 +.((((........)))). ( -5.30) 9 +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))). (-22.00) 1 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA + (-22.00) +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +.(((..((((....)))).))). ( -3.40) 51 +.(((((.......))))). ( -9.20) 48 +.((((..(((((.......)))))..)))). (-11.80) 42 +.(((((...(((((.......))))).))))). (-13.30) 40 +.((.(((((...(((((.......))))).))))))) (-13.60) 37 +.(((((.......))))). ( -7.70) 26 +.((.(((((.......)))))..)). ( -8.60) 23 +.(((.(((((.(((((.......)))))..)).(((((.......)))))..)))))) (-20.00) 16 +.((((........)))). ( -5.30) 9 +(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). (-29.60) 1 +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA + (-29.60) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnapaln_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnapaln_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnapaln_input1.fas --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnapaln_input1.fas Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,1 @@ +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA&UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnapaln_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnapaln_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,7 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +68.8844 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +,((((((..((((........)))).(((((.......))))).....(((((.......))))))))),)), +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA +(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnapdist_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnapdist_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnapdist_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnapdist_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,3 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +8.66253 |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnapkplex_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnapkplex_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnapkplex_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnapkplex_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| [ |
| @@ -0,0 +1,6 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).[[[.(((((]]]....))))))))))).. (-31.90) +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA +(((((((..((((..[[[.[[)))).(((((.]].]]]))))).....(((((.......)))))))))))). (-38.54) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaplex_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaplex_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaplex_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaplex_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +.((((((..((......((((...((((..(((.((((((.......((((((..(..((((...(((((((.&)))))))..........))))..)..))))))..))))))..)))..)))).......)))).)))))))). 1,73 : 1,72 (-40.70) i:72,j:1 <-41.26> + |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaplfold_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaplfold_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaplfold_result1.ps --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaplfold_result1.ps Tue Dec 06 12:33:48 2016 -0500 |
| [ |
| @@ -0,0 +1,242 @@ +%!PS-Adobe-3.0 EPSF-3.0 +%%Title: RNA Dot Plot +%%Creator: ViennaRNA-2.2.10 +%%CreationDate: Tue Oct 4 14:45:55 2016 +%%BoundingBox: 66 530 520 650 +%%DocumentFonts: Helvetica +%%Pages: 1 +%%EndComments + +%Options: +%This file contains the square roots of the base pair probabilities in the form +% i j sqrt(p(i,j)) ubox + +%%BeginProlog +/DPdict 100 dict def +DPdict begin +/logscale false def +/lpmin 1e-05 log def + +/box { %size x y box - draws box centered on x,y + 2 index 0.5 mul sub % x -= 0.5 + exch 2 index 0.5 mul sub exch % y -= 0.5 + 3 -1 roll dup rectfill +} bind def + +/ubox { + logscale { + log dup add lpmin div 1 exch sub dup 0 lt { pop 0 } if + } if + 3 1 roll + exch len exch sub 1 add box +} bind def + +/lbox { + 3 1 roll + len exch sub 1 add box +} bind def + +/drawseq { +% print sequence along all 4 sides +[ [0.7 -0.3 0 ] + [0.7 0.7 len add 0] + [-0.3 len sub -0.4 -90] + [-0.3 len sub 0.7 len add -90] +] { + gsave + aload pop rotate translate + 0 1 len 1 sub { + dup 0 moveto + sequence exch 1 getinterval + show + } for + grestore + } forall +} bind def + +/drawgrid{ + 0.01 setlinewidth + len log 0.9 sub cvi 10 exch exp % grid spacing + dup 1 gt { + dup dup 20 div dup 2 array astore exch 40 div setdash + } { [0.3 0.7] 0.1 setdash } ifelse + 0 exch len { + dup dup + 0 moveto + len lineto + dup + len exch sub 0 exch moveto + len exch len exch sub lineto + stroke + } for + [] 0 setdash + 0.04 setlinewidth + currentdict /cutpoint known { + cutpoint 1 sub + dup dup -1 moveto len 1 add lineto + len exch sub dup + -1 exch moveto len 1 add exch lineto + stroke + } if + 0.5 neg dup translate +} bind def + +end +%%EndProlog +DPdict begin +%delete next line to get rid of title +270 665 moveto /Helvetica findfont 14 scalefont setfont (Anolis_carolinensis_chrUn_GL343590.trna2-A) show + +/sequence { (\ +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA\ +) } def +/winSize 70 def +/len { sequence length } bind def + +292 416 translate +72 6 mul len 1 add winSize add 2 sqrt mul div dup scale +/Helvetica findfont 0.95 scalefont setfont + +/drawseq_turn {% print sequence at bottom + gsave + len 2 sqrt div dup neg 0.28 add exch 0.78 sub translate + 0 1 len 1 sub { + dup dup 2 sqrt mul 0 moveto + sequence exch 1 getinterval + show + } for + grestore +} bind def +/drawgrid_turn{ + 0.01 setlinewidth + len log 0.9 sub cvi 10 exch exp % grid spacing + dup 1 gt { + dup dup 20 div dup 2 array astore exch 40 div setdash + } { [0.3 0.7] 0.1 setdash } ifelse + 0 exch len { %for (0, gridspacing, len) + dup dup %duplicate what - gridspacing?? + dup len exch sub moveto %moveto diagonal? + dup winSize gt + {dup dup len exch sub winSize add lineto} + {dup len lineto}ifelse + dup len exch sub moveto %moveto diagonal? + dup len winSize sub le + {dup dup len exch sub dup winSize exch sub len add exch lineto} + {dup dup len exch sub len exch lineto}ifelse stroke pop pop + } for + len log 0.9 sub cvi 10 exch exp % grid spacing + dup 1 gt { + dup dup 20 div dup 2 array astore exch 40 div setdash + } { [0.3 0.7] 0.1 setdash } ifelse + 0 exch len { %for (0, gridspacing, len) + dup dup %duplicate what - gridspacing?? + dup len exch sub moveto %moveto diagonal? + len exch sub 0.7 sub exch 0.7 sub exch lineto + stroke + }for + winSize len moveto len winSize lineto stroke + [] 0 setdash + 0.04 setlinewidth + currentdict /cutpoint known { + cutpoint 1 sub + dup dup -1 moveto len 1 add lineto + len exch sub dup + -1 exch moveto len 1 add exch lineto + stroke + } if + 0.5 neg dup translate +} bind def + +0.5 dup translate +drawseq_turn +45 rotate + + +%draw the grid +drawgrid_turn + +%start of base pair probability data +2 70 0.1568 ubox +2 71 0.9619 ubox +3 69 0.1395 ubox +3 70 0.7414 ubox +3 72 0.8060 ubox +4 68 0.1157 ubox +4 69 0.6748 ubox +4 71 0.4682 ubox +5 67 0.1065 ubox +5 68 0.6724 ubox +5 70 0.3765 ubox +6 47 0.1250 ubox +6 67 0.6497 ubox +7 46 0.1273 ubox +7 66 0.6008 ubox +8 45 0.1294 ubox +8 48 0.2252 ubox +9 47 0.2335 ubox +10 25 0.7863 ubox +11 24 0.7884 ubox +11 43 0.5215 ubox +11 45 0.2330 ubox +12 23 0.7883 ubox +12 42 0.5493 ubox +12 44 0.2317 ubox +13 22 0.7882 ubox +13 41 0.5528 ubox +13 43 0.2306 ubox +14 40 0.5338 ubox +15 20 0.1027 ubox +16 38 0.5325 ubox +16 40 0.1646 ubox +17 37 0.5639 ubox +17 39 0.1633 ubox +18 36 0.5591 ubox +18 38 0.1125 ubox +19 36 0.2406 ubox +20 34 0.5341 ubox +20 35 0.2514 ubox +21 33 0.4064 ubox +22 32 0.2491 ubox +22 33 0.4413 ubox +23 32 0.5541 ubox +24 31 0.6100 ubox +25 30 0.6092 ubox +26 36 0.2144 ubox +27 35 0.2146 ubox +27 43 0.7269 ubox +28 34 0.2043 ubox +28 42 0.7445 ubox +29 41 0.7467 ubox +30 40 0.7468 ubox +31 39 0.7470 ubox +38 73 0.3242 ubox +39 72 0.3055 ubox +40 73 0.1211 ubox +41 71 0.2953 ubox +41 72 0.1146 ubox +42 70 0.2588 ubox +43 69 0.2613 ubox +43 71 0.2848 ubox +44 68 0.2587 ubox +44 70 0.3418 ubox +45 67 0.2339 ubox +45 68 0.1772 ubox +45 69 0.3617 ubox +45 70 0.1014 ubox +45 72 0.3111 ubox +46 67 0.2550 ubox +46 68 0.3039 ubox +46 69 0.1150 ubox +46 71 0.2540 ubox +47 66 0.2925 ubox +48 67 0.1378 ubox +49 65 0.9938 ubox +50 64 0.9971 ubox +51 63 0.9980 ubox +52 62 0.9980 ubox +53 61 0.9963 ubox +54 59 0.1628 ubox +55 60 0.1625 ubox +showpage +end +%%EOF |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaplfold_result2.ps --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaplfold_result2.ps Tue Dec 06 12:33:48 2016 -0500 |
| [ |
| @@ -0,0 +1,212 @@ +%!PS-Adobe-3.0 EPSF-3.0 +%%Title: RNA Dot Plot +%%Creator: ViennaRNA-2.2.10 +%%CreationDate: Tue Oct 4 14:45:55 2016 +%%BoundingBox: 66 530 520 650 +%%DocumentFonts: Helvetica +%%Pages: 1 +%%EndComments + +%Options: +%This file contains the square roots of the base pair probabilities in the form +% i j sqrt(p(i,j)) ubox + +%%BeginProlog +/DPdict 100 dict def +DPdict begin +/logscale false def +/lpmin 1e-05 log def + +/box { %size x y box - draws box centered on x,y + 2 index 0.5 mul sub % x -= 0.5 + exch 2 index 0.5 mul sub exch % y -= 0.5 + 3 -1 roll dup rectfill +} bind def + +/ubox { + logscale { + log dup add lpmin div 1 exch sub dup 0 lt { pop 0 } if + } if + 3 1 roll + exch len exch sub 1 add box +} bind def + +/lbox { + 3 1 roll + len exch sub 1 add box +} bind def + +/drawseq { +% print sequence along all 4 sides +[ [0.7 -0.3 0 ] + [0.7 0.7 len add 0] + [-0.3 len sub -0.4 -90] + [-0.3 len sub 0.7 len add -90] +] { + gsave + aload pop rotate translate + 0 1 len 1 sub { + dup 0 moveto + sequence exch 1 getinterval + show + } for + grestore + } forall +} bind def + +/drawgrid{ + 0.01 setlinewidth + len log 0.9 sub cvi 10 exch exp % grid spacing + dup 1 gt { + dup dup 20 div dup 2 array astore exch 40 div setdash + } { [0.3 0.7] 0.1 setdash } ifelse + 0 exch len { + dup dup + 0 moveto + len lineto + dup + len exch sub 0 exch moveto + len exch len exch sub lineto + stroke + } for + [] 0 setdash + 0.04 setlinewidth + currentdict /cutpoint known { + cutpoint 1 sub + dup dup -1 moveto len 1 add lineto + len exch sub dup + -1 exch moveto len 1 add exch lineto + stroke + } if + 0.5 neg dup translate +} bind def + +end +%%EndProlog +DPdict begin +%delete next line to get rid of title +270 665 moveto /Helvetica findfont 14 scalefont setfont (Anolis_carolinensis_chrUn_GL343207.trna3-A) show + +/sequence { (\ +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA\ +) } def +/winSize 70 def +/len { sequence length } bind def + +292 416 translate +72 6 mul len 1 add winSize add 2 sqrt mul div dup scale +/Helvetica findfont 0.95 scalefont setfont + +/drawseq_turn {% print sequence at bottom + gsave + len 2 sqrt div dup neg 0.28 add exch 0.78 sub translate + 0 1 len 1 sub { + dup dup 2 sqrt mul 0 moveto + sequence exch 1 getinterval + show + } for + grestore +} bind def +/drawgrid_turn{ + 0.01 setlinewidth + len log 0.9 sub cvi 10 exch exp % grid spacing + dup 1 gt { + dup dup 20 div dup 2 array astore exch 40 div setdash + } { [0.3 0.7] 0.1 setdash } ifelse + 0 exch len { %for (0, gridspacing, len) + dup dup %duplicate what - gridspacing?? + dup len exch sub moveto %moveto diagonal? + dup winSize gt + {dup dup len exch sub winSize add lineto} + {dup len lineto}ifelse + dup len exch sub moveto %moveto diagonal? + dup len winSize sub le + {dup dup len exch sub dup winSize exch sub len add exch lineto} + {dup dup len exch sub len exch lineto}ifelse stroke pop pop + } for + len log 0.9 sub cvi 10 exch exp % grid spacing + dup 1 gt { + dup dup 20 div dup 2 array astore exch 40 div setdash + } { [0.3 0.7] 0.1 setdash } ifelse + 0 exch len { %for (0, gridspacing, len) + dup dup %duplicate what - gridspacing?? + dup len exch sub moveto %moveto diagonal? + len exch sub 0.7 sub exch 0.7 sub exch lineto + stroke + }for + winSize len moveto len winSize lineto stroke + [] 0 setdash + 0.04 setlinewidth + currentdict /cutpoint known { + cutpoint 1 sub + dup dup -1 moveto len 1 add lineto + len exch sub dup + -1 exch moveto len 1 add exch lineto + stroke + } if + 0.5 neg dup translate +} bind def + +0.5 dup translate +drawseq_turn +45 rotate + + +%draw the grid +drawgrid_turn + +%start of base pair probability data +2 70 0.1914 ubox +2 71 0.9787 ubox +3 69 0.1660 ubox +3 70 0.7676 ubox +3 72 0.8449 ubox +4 68 0.1377 ubox +4 69 0.7113 ubox +4 71 0.4928 ubox +5 67 0.1266 ubox +5 68 0.7090 ubox +5 70 0.3955 ubox +6 67 0.6860 ubox +7 66 0.6345 ubox +10 25 0.9888 ubox +11 24 0.9915 ubox +12 23 0.9914 ubox +13 22 0.9913 ubox +15 20 0.1291 ubox +17 73 0.1217 ubox +18 72 0.1091 ubox +27 43 0.9522 ubox +28 42 0.9836 ubox +29 41 0.9874 ubox +30 40 0.9886 ubox +31 39 0.9883 ubox +33 37 0.1014 ubox +38 73 0.1170 ubox +39 72 0.1080 ubox +41 71 0.1523 ubox +42 70 0.1341 ubox +43 69 0.1387 ubox +43 71 0.2800 ubox +44 68 0.1392 ubox +44 70 0.3364 ubox +45 67 0.1266 ubox +45 68 0.1838 ubox +45 69 0.3586 ubox +45 72 0.3252 ubox +46 67 0.2524 ubox +46 68 0.3007 ubox +46 69 0.1142 ubox +46 71 0.2655 ubox +47 66 0.2807 ubox +48 67 0.1394 ubox +49 65 0.9981 ubox +50 64 0.9997 ubox +51 63 0.9997 ubox +52 62 0.9997 ubox +53 61 0.9981 ubox +54 59 0.1565 ubox +55 60 0.3242 ubox +showpage +end +%%EOF |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaplot_input1.dbn --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaplot_input1.dbn Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,3 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-21.90) |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaplot_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaplot_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnasnoop_input1a.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnasnoop_input1a.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +>homo +CGCGACCUCAGAUCAGACGUGGCGACCCGCUGAAUUUAAGCAUAUUAGUCAGCGGAGGAAAAGAAACUAACCAGGAUUCCCUCAGUAACGGCGAGUGAACAGGGAAGAGCCCAGCGCCGAAUCCCCGCCCCGCGGGGCGCGGGACAUGUGGCGUACGGAAGACCCGCUCCCCGGCGCCGCUCGUGGGGGGCCCAAGUCCUUCUGAUCGAGGCCCAGCCCGUGGACGGUGUGAGGCCGGUAGCGGCCGGCGCGCGCCCGGGUCUUCCCGGAGUCGGGUUGCUUGGGAAUGCAGCCCAAAGCGGGUGGUAAACUCCAUCUAAGGCUAAAUACCGGCACGAGACCGAUAGUCAACAAGUACCGUAAGGGAAAGUUGAAAAGAACUUUGAAGAGAGAGUUCAAGAGGGCGUGAAACCGUUAAGAGGUAAACGGGUGGGGUCCGCGCAGUCCGCCCGGAGGAUUCAACCCGGCGGCGGGUCCGGCCGUGUCGGCGGCCCGGCGGAUCUUUCCCGCCCCCCGUUCCUCCCGACCCCUCCACCCGCCCUCCCUUCCCCCGCCGCCCCUCCUCCUCCUCCCCGGAGGGGGCGGGCUCCGGCGGGUGCGGGGGUGGGCGGGCGGGGCCGGGGGUGGGGUCGGCGGGGGACCGUCCCCCGACCGGCGACCGGCCGCCGCCGGGCGCAUUUCCACCGCGGCGGUGCGCCGCGACCGGCUCCGGGACGGCUGGGAAGGCCCGGCGGGGAAGGUGGCUCGGGGGGCCCCGUCCGUCCGUCCGUCCUCCUCCUCCCCCGUCUCCGCCCCCCGGCCCCGCGUCCUCCCUCGGGAGGGCGCGCGGGUCGGGGCGGCGGCGGCGGCGGCGGUGGCGGCGGCGGCGGGGGCGGCGGGACCGAAACCCCCCCCGAGUGUUACAGCCCCCCCGGCAGCAGCACUCGCCGAAUCCCGGGGCCGAGGGAGCGAGACCCGUCGCCGCGCUCUCCCCCCUCCCGGCGCCCACCCCCGCGGGGAAUCCCCCGCGAGGGGGGUCUCCCCCGCGGGGGCGCGCCGGCGUCUCCUCGUGGGGGGGCCGGGCCACCCCUCCCACGGCGCGACCGCUCUCCCACCCCUCCUCCCCGCGCCCCCGCCCCGGCGACGGGGGGGGUGCCGCGCGCGGGUCGGGGGGCGGGGCGGACUGUCCCCAGUGCGCCCCGGGCGGGUCGCGCCGUCGGGCCCGGGGGAGGUUCUCUCGGGGCCACGCGCGCGUCCCCCGAAGAGGGGGACGGCGGAGCGAGCGCACGGGGUCGGCGGCGACGUCGGCUACCCACCCGACCCGUCUUGAAACACGGACCAAGGAGUCUAACACGUGCGCGAGUCGGGGGCUCGCACGAAAGCCGCCGUGGCGCAAUGAAGGUGAAGGCCGGCGCGCUCGCCGGCCGAGGUGGGAUCCCGAGGCCUCUCCAGUCCGCCGAGGGCGCACCACCGGCCCGUCUCGCCCGCCGCGCCGGGGAGGUGGAGCACGAGCGCACGUGUUAGGACCCGAAAGAUGGUGAACUAUGCCUGGGCAGGGCGAAGCCAGAGGAAACUCUGGUGGAGGUCCGUAGCGGUCCUGACGUGCAAAUCGGUCGUCCGACCUGGGUAUAGGGGCGAAAGACUAAUCGAACCAUCUAGUAGCUGGUUCCCUCCGAAGUUUCCCUCAGGAUAGCUGGCGCUCUCGCAGACCCGACGCACCCCCGCCACGCAGUUUUAUCCGGUAAAGCGAAUGAUUAGAGGUCUUGGGGCCGAAACGAUCUCAACCUAUUCUCAAACUUUAAAUGGGUAAGAAGCCCGGCUCGCUGGCGUGGAGCCGGGCGUGGAAUGCGAGUGCCUAGUGGGCCACUUUUGGUAAGCAGAACUGGCGCUGCGGGAUGAACCGAACGCCGGGUUAAGGCGCCCGAUGCCGACGCUCAUCAGACCCCAGAAAAGGUGUUGGUUGAUAUAGACAGCAGGACGGUGGCCAUGGAAGUCGGAAUCCGCUAAGGAGUGUGUAACAACUCACCUGCCGAAUCAACUAGCCCUGAAAAUGGAUGGCGCUGGAGCGUCGGGCCCAUACCCGGCCGUCGCCGGCAGUCGAGAGUGGACGGGAGCGGCGGGGGCGGCGCGCGCGCGCGCGCGUGUGGUGUGCGUCGGAGGGCGGCGGCGGCGGCGGCGGCGGGGGUGUGGGGUCCUUCCCCCGCCCCCCCCCCCACGCCUCCUCCCCUCCUCCCGCCCACGCCCCGCUCCCCGCCCCCGGAGCCCCGCGGACGCUACGCCGCGACGAGUAGGAGGGCCGCUGCGGUGAGCCUUGAAGCCUAGGGCGCGGGCCCGGGUGGAGCCGCCGCAGGUGCAGAUCUUGGUGGUAGUAGCAAAUAUUCAAACGAGAACUUUGAAGGCCGAAGUGGAGAAGGGUUCCAUGUGAACAGCAGUUGAACAUGGGUCAGUCGGUCCUGAGAGAUGGGCGAGCGCCGUUCCGAAGGGACGGGCGAUGGCCUCCGUUGCCCUCGGCCGAUCGAAAGGGAGUCGGGUUCAGAUCCCCGAAUCCGGAGUGGCGGAGAUGGGCGCCGCGAGGCGUCCAGUGCGGUAACGCGACCGAUCCCGGAGAAGCCGGCGGGAGCCCCGGGGAGAGUUCUCUUUUCUUUGUGAAGGGCAGGGCGCCCUGGAAUGGGUUCGCCCCGAGAGAGGGGCCCGUGCCUUGGAAAGCGUCGCGGUUCCGGCGGCGUCCGGUGAGCUCUCGCUGGCCCUUGAAAAUCCGGGGGAGAGGGUGUAAAUCUCGCGCCGGGCCGUACCCAUAUCCGCAGCAGGUCUCCAAGGUGAACAGCCUCUGGCAUGUUGGAACAAUGUAGGUAAGGGAAGUCGGCAAGCCGGAUCCGUAACUUCGGGAUAAGGAUUGGCUCUAAGGGCUGGGUCGGUCGGGCUGGGGCGCGAAGCGGGGCUGGGCGCGCGCCGCGGCUGGACGAGGCGCGCGCCCCCCCCACGCCCGGGGCACCCCCCUCGCGGCCCUCCCCCGCCCCACCCGCGCGCGCCGCUCGCUCCCUCCCCACCCCGCGCCCUCUCUCUCUCUCUCUCCCCCGCUCCCCGUCCUCCCCCCUCCCCGGGGGAGCGCCGCGUGGGGGCGCGGCGGGGGGAGAAGGGUCGGGGCGGCAGGGGCCGCGCGGCGGCCGCCGGGGCGGCCGGCGGGGGCAGGUCCCCGCGAGGGGGGCCCCGGGGACCCGGGGGGCCGGCGGCGGCGCGGACUCUGGACGCGAGCCGGGCCCUUCCCGUGGAUCGCCCCAGCUGCGGCGGGCGUCGCGGCCGCCCCCGGGGAGCCCGGCGGCGGCGCGGCGCGCCCCCCACCCCCACCCCACGUCUCGGUCGCGCGCGCGUCCGCUGGGGGCGGGAGCGGUCGGGCGGCGGCGGUCGGCGGGCGGCGGGGCGGGGCGGUUCGUCCCCCCGCCCUACCCCCCCGGCCCCGUCCGCCCCCCGUUCCCCCCUCCUCCUCGGCGCGCGGCGGCGGCGGCGGCAGGCGGCGGAGGGGCCGCGGGCCGGUCCCCCCCGCCGGGUCCGCCCCCGGGGCCGCGGUUCCGCGCGCGCCUCGCCUCGGCCGGCGCCUAGCAGCCGACUUAGAACUGGUGCGGACCAGGGGAAUCCGACUGUUUAAUUAAAACAAAGCAUCGCGAAGGCCCGCGGCGGGUGUUGACGCGAUGUGAUUUCUGCCCAGUGCUCUGAAUGUCAAAGUGAAGAAAUUCAAUGAAGCGCGGGUAAACGGCGGGAGUAACUAUGACUCUCUUAAGGUAGCCAAAUGCCUCGUCAUCUAAUUAGUGACGCGCAUGAAUGGAUGAACGAGAUUCCCACUGUCCCUACCUACUAUCCAGCGAAACCACAGCCAAGGGAACGGGCUUGGCGGAAUCAGCGGGGAAAGAAGACCCUGUUGAGCUUGACUCUAGUCUGGCACGGUGAAGAGACAUGAGAGGUGUAGAAUAAGUGGGAGGCCCCCGGCGCCCCCCCGGUGUCCCCGCGAGGGGCCCGGGGCGGGGUCCGCGGCCCUGCGGGCCGCCGGUGAAAUACCACUACUCUGAUCGUUUUUUCACUGACCCGGUGAGGCGGGGGGGCGAGCCCGAGGGGCUCUCGCUUCUGGCGCCAAGCGCCCGCCCGGCCGGGCGCGACCCGCUCCGGGGACAGUGCCAGGUGGGGAGUUUGACUGGGGCGGUACACCUGUCAAACGGUAACGCAGGUGUCCUAAGGCGAGCUCAGGGAGGACAGAAACCUCCCGUGGAGCAGAAGGGCAAAAGCUCGCUUGAUCUUGAUUUUCAGUACGAAUACAGACCGUGAAAGCGGGGCCUCACGAUCCUUCUGACCUUUUGGGUUUUAAGCAGGAGGUGUCAGAAAAGUUACCACAGGGAUAACUGGCUUGUGGCGGCCAAGCGUUCAUAGCGACGUCGCUUUUUGAUCCUUCGAUGUCGGCUCUUCCUAUCAUUGUGAAGCAGAAUUCGCCAAGCGUUGGAUUGUUCACCCACUAAUAGGGAACGUGAGCUGGGUUUAGACCGUCGUGAGACAGGUUAGUUUUACCCUACUGAUGAUGUGUUGUUGCCAUGGUAAUCCUGCUCAGUACGAGAGGAACCGCAGGUUCAGACAUUUGGUGUAUGUGCUUGGCUGAGGAGCCAAUGGGGCGAAGCUACCAUCUGUGGGAUUAUGACUGAACGCCUCUAAGUCAGAAUCCCGCCCAGGCGAACGAUACGGCAGCGCCGCGGAGCCUCGGUUGGCCUCGGAUAGCCGGUCCCCCGCCUGUCCCCGCCGGCGGGCCGCCCCCCCCUCCACGCGCCCCGCCGCGGGAGGGCGCGUGCCCCGCCGCGCGCCGGGACCGGGGUCCGGUGCGGAGUGCCCUUCGUCCUGGGAAACGGGGCGCGGCCGGAAAGGCGGCCGCCCCCUCGCCCGUCACGCACCGCACGUUCGUGGGGAACCUGGCGCUAAACCAUUCGUAGACGACCUGCUUCUGGGUCGGGGUUUCGUACGUAGCAGAGCAGCUCCCUCGCUGCGAUCUAUUGAAAGUCAGCCCUCGACACAAGGGUUUGUC |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnasnoop_input1b.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnasnoop_input1b.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +>ACA51 +AACCTACCCCATATACACCTCAGCTCAGGCCCTGTGCCTGGTCTGTATTGTGAATGGGGGAACATAG |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnasnoop_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnasnoop_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>ACA51 +>homo +<<<<<.<<<<.|.<<.<<<&...(((>>>.>>.((((...............................))))..>>>>.>>>>>))) 4468,4486;4479 : 8,65 (-35.60 = -16.10 + -7.60 + -12.40 + -3.60 + 4.1 ) (-23.60) +UGUUCACCCACUAAUAGGG&AACCUACCCCAUAUACACCUCAGCUCAGGCCCUGUGCCUGGUCUGUAUUGUGAAUGGGGGAACAUAG |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnasubopt_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnasubopt_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnasubopt_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnasubopt_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| [ |
| @@ -0,0 +1,7 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-A [0] +UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA -22.00 0.00 +(((((((..((((........)))).(((((.......))))).....(((((.......))))))))).))) -22.00 +.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. -22.00 +>Anolis_carolinensis_chrUn_GL343207.trna3-A [0] +GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA -29.60 0.00 +(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). -29.60 |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaup_input1.fa --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaup_input1.fa Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,4 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872) Ala (AGC) 73 bp Sc: 49.55 +TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC (1513626-1513698) Ala (AGC) 73 bp Sc: 56.15 +GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA |
| b |
| diff -r 000000000000 -r 204093ddb00a test-data/rnaup_result1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rnaup_result1.txt Tue Dec 06 12:33:48 2016 -0500 |
| b |
| @@ -0,0 +1,6 @@ +>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC + 55, 58 (0.019) for u= 4 +RNAup output in file: Anolis_carolinensis_chrUn_GL343590.trna2-A_u1.out +>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC + 16, 19 (0.004) for u= 4 +RNAup output in file: Anolis_carolinensis_chrUn_GL343207.trna3-A_u1.out |
| b |
| diff -r 000000000000 -r 204093ddb00a vienna_rna.tar.gz |
| b |
| Binary file vienna_rna.tar.gz has changed |