Mercurial > repos > rnateam > viennarna_rnacofold

changeset 0:f3d4656ed57c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/vienna_rna commit 0065dafe7bbd382bb995b28cc4089c9e4f4eeeb9
author rnateam
date Tue, 06 Dec 2016 12:33:31 -0500
parents
children cb8201484ad8
files macros.xml rnacofold.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
diffstat 51 files changed, 960 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml	Tue Dec 06 12:33:31 2016 -0500
@@ -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>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/rnacofold.xml	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,182 @@
+<tool id="viennarna_rnacofold" name="@EXECUTABLE@" version="@VERSION@.0">
+    <description>Calculate secondary structures of two RNAs with dimerization</description>
+    <macros>
+        <token name="@EXECUTABLE@">RNAcofold</token>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements" />
+    <expand macro="stdio" />
+    <expand macro="version_command" />
+    <command>
+<![CDATA[
+        RNAcofold  < '$input' > '$tabularFile'
+        -T$model_options.temperature -d$model_options.dangling        
+        $general_options.noPS
+        $general_options.noconversion        
+        $IDs.auto_id
+        #if str($IDs.id_prefix) <> "alignment"
+            --id-prefix='$IDs.id_prefix'
+        #end if
+        #if $IDs.id_digits <> 4
+            --id-digits=$IDs.id_digits
+        #end if
+        #if $IDs.id_start <> 1
+            --id-start=$IDs.id_start
+        #end if        
+        #if str($constraints.constraintLocation.constraintSelector) == "fromFile"
+            --constraint='$constraints.constraintLocation.constraintsFile'
+            $constraints.constraintLocation.canonicalBPonly
+            $constraints.constraintLocation.enforceConstraint                
+        #end if
+        $algorithm_options.allpf
+        #if len($algorithm_options.concfiles) == 1
+            #for $i in $algorithm_options.concfiles
+                --concfile='$i.concfile'
+            #end for
+        #end if
+        #if $algorithm_options.pf <> "-1"
+            --partfunc=$algorithm_options.pf
+        #end if
+        #if $algorithm_options.pfScale <> 1.07
+            --pfScale=$algorithm_options.pfScale
+        #end if
+        #if $algorithm_options.bppmThreshold <> 1e-5
+            --bppmThreshold=$algorithm_options.bppmThreshold
+        #end if
+        $algorithm_options.gquad		
+        $model_options.notetra
+        $model_options.nolp
+        $model_options.nogu
+        $model_options.noclosinggu
+        #if $model_options.nsp
+            --nsp='$model_options.nsp'
+        #end if
+        #if $model_options.betaScale <> 1.0
+            --betaScale=$model_options.betaScale
+        #end if
+        
+]]>
+    </command>
+    <inputs>
+        <param format="fasta" name="input" type="data" label="Fasta file"/>
+        <section name="general_options" title="General Options">
+            <param name="noPS" type="boolean" truevalue="" falsevalue="--noPS" checked="true" label="Produce postscript output" help="" argument="--noPS"/>
+            <param name="noconversion" type="boolean" truevalue="" falsevalue="--noconv" checked="true" label="Convert Thymine to Uracil (T -> U)" argument="--noconv" help="Avoids confusion with DNA sequences"/>
+        </section>
+        <section name="IDs" title="Naming Conventions">		
+                <param name="auto_id" 
+                    type="boolean" truevalue="--auto-id" falsevalue="" checked="false" 
+                    label="Automatically generate an ID for each alignment." 
+                    help="If this flag is active, RNAalifold ignores any IDs retrieved from the input and automatically generates an ID for each alignment." 
+                    argument="--auto-id"/>
+                <param name="id_prefix" 
+                    type="text" value="alignment" 
+                    label="Prefix for automatically generated IDs (as used in output file names)" 
+                    help="If this parameter is set, each alignment will be prefixed with the provided string. Hence, the output files will obey the following naming scheme: 'prefix_xxxx_ss.ps' (secondary structure plot), 'prefix_xxxx_dp.ps' (dot−plot), 'prefix_xxxx_aln.ps' (annotated alignment), etc. where xxxx is the alignment number beginning with the second alignment in the input. Use this setting in conjunction with the −−continuous−ids flag to assign IDs beginning with the first input alignment." 
+                    argument="--id-prefix"/>
+                <param name="id_digits" 
+                    type="integer" value="4" min="1" max="18" 
+                    label="The number of digits of the counter in automatically generated alignment IDs" 
+                    help="When alignments IDs are automatically generated, they receive an increasing number, starting with 1. This number will always be left−padded by leading zeros, such that the number takes up a certain width. Using this parameter, the width can be specified to the users need. We allow numbers in the range [1:18]." 
+                    argument="--id-digits"/>
+                <param name="id_start" 
+                    type="integer" value="1" min="0" 
+                    label="First number in automatically generated alignment IDs" 
+                    help="When alignment IDs are automatically generated, they receive an increasing number, usually starting with 1. Using this parameter, the first number can be specified to the users requirements. Note: negative numbers are not allowed. Note: Setting this parameter implies continuous alignment IDs, i.e. it activates the −−continuous−ids flag.." 
+                    argument="--id-start"/>
+        </section>
+        <section name="model_options" title="Model Options" expanded="true">
+            <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: only unpaired bases can participate in at most one dangling end, this is the default for mfe folding but unsupported for the partition function folding</option>
+                <option value="2" selected="True">2: unpaired bases participate in all dangling ends</option>
+                <option value="3">3: mfe folding will allow coaxial stacking of adjacent helices in multi−loops.</option>
+            </param>
+            <param name="nolp" type="boolean" truevalue="" falsevalue="--noLP" checked="true" label="Allow lonely base-pairs" help="(--noLP)"/>
+            <param name="nogu" type="boolean" truevalue="" falsevalue="--noGU" checked="true" label="Allow GU pairing" help="--noGU"/>
+            <param name="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. --noClosingGU"/>
+            <param name="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. (--noTetra)"/>
+            <param name="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." 
+               argument="--nsp"/>
+            <param name="betaScale" type="float" value="1.0" label="Scaling of Boltzman factors" 
+                help=" The argument provided with this option enables to scale the thermodynamic temperature used in the Boltzmann factors independently from the temperature used to scale the individual energy contributions of the loop types." 
+                argument="--betaScale"/>
+        </section>
+        <section name="constraints" title="Structure constraints">
+                <conditional name="constraintLocation">
+                    <param name="constraintSelector" type="select" label="Constraints">
+                    <!--	<option value="fromInput">The constraints are included in the input file</option> -->
+                        <option value="fromFile">The constraints are in a seperate file</option>
+                        <option value="none" selected="true">Don't use constraints</option>
+                    </param>
+                    <!-- <when value="fromInput"></when> -->
+                    <when value="none"></when>
+                    <when value="fromFile">
+                        <param name="constraintsFile" type="data" format="*" label="Constraints file" argument="--constraint"/>
+                        <param argument="--canonicalBPonly" type="boolean" truevalue="--canonicalBPonly" falsevalue="" checked="false" label="Remove non-canonical base pairs from he structure constraint" />
+                        <param argument="--enforceConstraint" type="boolean" truevalue="--enforceConstraint" falsevalue="" checked="false" label="Enforce base pair given by round brackets () in structure constraint" />
+                    </when>
+                </conditional>				
+        </section>
+        <section name="algorithm_options" title="Algorithm Options">
+            <param name="pf" type="select" label="Calculate partition function" help="Calculate the partition function and base pairing probability matrix in addition to the mfe structure. Default is calculation of mfe structure only." argument="-p">
+                <option value="-1" selected="true">None</option>
+                <option value="0">0: Deactivates the calculation of the pair probabilities, saving about 50% in runtime</option>
+                <option value="1" >1: Calculate the partition function and base pairing probability matrix</option>				
+            </param>
+            <param name="allpf" type="boolean" checked="false" truevalue="--all_pf" falsevalue="" label="Calculate homo-dimers as well as A and B monomers" help="--all_pf"/>
+            <param argument="-c" type="boolean" checked="false" truevalue="-c" falsevalue="" label="Concentrations" help="In addition to everything listed under the -a option, read in initial monomer concentrations and compute the expected equilibrium concentrations of the 5 possible species (AB, AA, BB, A, B)."/>
+            <repeat name="concfiles" title="A file with initial concentrations for the two sequences" default="0" max="1" help="The table consits of arbitrary many lines with just two numbers (the concentration of sequence A and B). This option will automatically toggle the −c (and thus -a and --all_pf) options">
+                <param type="data" name="concfile" format="*" label="Concentrations file"/>
+            </repeat>			
+            <param argument="--pfScale" type="float" value="1.07" label="Scaling factor" help="In the calculation of the pf use scale*mfe as an estimate for the ensemble free energy (used to avoid overflows). The default is 1.07, useful values are 1.0 to 1.2. Occasionally needed for long sequences." argumen="--pfScale"/>
+            <param argument="--bppmThreshold" type="float" value="1e-5" label="Threshold for base pair probabilities" help="By setting the threshold the base pair probabilities that are included in the output can be varied. By default only those exceeding 1e−5 in probability will be shown as squares in the dot plot. Changing the threshold to any other value allows for increase or decrease of data." argumen="--bppmThreshold"/>
+            <param argument="--gquad" type="boolean" truevalue="--gquad" falsevalue="" checked="false" label="G Quadruplex formation" help="take into account G Quadruplex formation"/>
+        </section>      
+    </inputs>
+    <outputs>
+        <data format="txt" name="tabularFile"/>
+        <collection name="sequence_outputs" type="list" label="rna_eps outputs">
+            <discover_datasets pattern="(?P&lt;designation&gt;.+)_ss\.ps" ext="eps" />
+        </collection>
+        <collection name="matrix_outputs" type="list" label="rna_eps outputs">
+            <filter>pf is True</filter>
+            <discover_datasets pattern="(?P&lt;designation&gt;.+)_dp\.ps" ext="rna_eps" visible="true"/>
+        </collection>
+    </outputs>
+    <tests>
+        <test>
+            <param name="input" value="rnacofold_input1.fas"/>
+            <output name="out_file" file="rnacofold_result1.txt"/>
+        </test>
+    </tests>
+    <help>
+<![CDATA[
+**RNAcofold**
+
+
+-----
+
+**Input format**
+
+RNAcofold requires one input file:
+
+- Fasta file
+
+Two different RNA sequences can be specified concatenated with the '&' character. The folding of the two sequences with each other will be calculated.
+
+
+------
+
+**Outputs**
+
+- Energies of the RNA sequences and the dimers
+- Structure images
+- Dot Plot Matrix images (if --partfunc is used)
+
+]]>
+    </help>
+    <expand macro="citations" />
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/kinfold_input.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,1 @@
+ACUGAUCGUAGUCAC
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rna2dfold_input1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,4 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872)  Ala (AGC) 73 bp  Sc: 49.55
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+.((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))..
+.((((((..(((..........))).(((((.......))))).....(((((.......)))))))))))..
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rna2dfold_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaaliduplex_input1.clustal	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,7 @@
+CLUSTAL 2.1 multiple sequence alignment
+
+
+Anolis_carolinensis_chrUn_GL34      TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA
+Anolis_carolinensis_chrUn_GL35      GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA
+                                     ************************************** ************** ******* ********* 
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaaliduplex_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,1 @@
+(((((((..(.......((((..(..((((((..((((((.......((((((..(..(((((..(((((((.&)))))))..).......))))..)..))))))..)))))).......))))))..)..)))))..))))))).   1,73  :   1,73  (-40.30)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaalifold_input1.clustal	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,7 @@
+CLUSTAL 2.1 multiple sequence alignment
+
+
+Anolis_carolinensis_chrUn_GL34      TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA
+Anolis_carolinensis_chrUn_GL35      GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA
+                                     ************************************** ************** ******* ********* 
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaalifold_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,2 @@
+GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUCGAUGCCCACAUUCUCCA
+(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). (-26.45 = -26.20 +  -0.25)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnacofold_input1.fas	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,2 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872)  Ala (AGC) 73 bp  Sc: 49.55
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA&UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnacofold_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,3 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-A
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA&UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+.((((((..((((........)))).(((((.......))))).....(((((((((((.(((((((((((..&.))))))..((((........)))).(((((.......))))).....))))))))))).))))))))))).. (-55.90)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnadistance_input1.dbn	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,2 @@
+.((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))..
+(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))).
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnadistance_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,1 @@
+f: 4  
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaduplex_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaduplex_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,3 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC
+>Anolis_carolinensis_chrUn_GL343207.trna3-AlaAGC
+.((((((..((......((((...((((..(((.((((((.......((((((..(..((((...(((((((.&)))))))..........))))..)..))))))..))))))..)))..)))).......)))).)))))))).   1,73  :   1,72  (-40.70)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaeval_input1.dbn	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,3 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872)  Ala (AGC) 73 bp  Sc: 49.55
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+.((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))..
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaeval_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,3 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872)  Ala (AGC) 73 bp  Sc: 49.55
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnafold_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnafold_input2.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnafold_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,6 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-A
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00)
+>Anolis_carolinensis_chrUn_GL343207.trna3-A
+GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA
+(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))). (-29.60)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnafold_result2.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,6 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-A
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+..(((.....((((....((..(((....)))..))...)))).....(((((.......)))))....))). ( -3.37)
+>Anolis_carolinensis_chrUn_GL343207.trna3-A
+GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUUGAUGCCCGCAUUCUCCA
+(((((..((.(((.....((..(((....)))..))......))).))(((((.......)))))..))))). ( -5.02)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnafold_result3.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,3 @@
+>Sequence
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-22.00)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaheat_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,2 @@
+> comment 1
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaheat_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnainverse_input1.clu	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,2 @@
+(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))).
+gggccccnnaauunnnnnnnnaauunGGGGcnnnnnnnAAAAAnnnnnuuuuannnnnnnuaaaaggggcccn
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnalalifold_input1.clustal	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,7 @@
+CLUSTAL 2.1 multiple sequence alignment
+
+
+Anolis_carolinensis_chrUn_GL34      TGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGTGAGAGGTAGTGGGATCGATGCCCACATTCTCCA
+Anolis_carolinensis_chrUn_GL35      GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTAGCATGCGAGAGGTAGCGGGATTGATGCCCGCATTCTCCA
+                                     ************************************** ************** ******* ********* 
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnalalifold_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,4 @@
+(((.(((((.(((((.......)))))..)).(((((.......)))))..)))))) (-19.05)   17 -   73
+((((........)))). ( -5.10)   10 -   26
+GGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGCGAGAGGUAGCGGGAUCGAUGCCCACAUUCUCCA
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnalfold_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnalfold_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnapaln_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnapaln_input1.fas	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,1 @@
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA&UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnapaln_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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
+(((((((..((((........)))).(((((.......))))).....(((((.......)))))))))))).
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnapdist_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnapdist_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnapkplex_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnapkplex_result1.txt	Tue Dec 06 12:33:31 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)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaplex_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaplex_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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>
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaplfold_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaplfold_result1.ps	Tue Dec 06 12:33:31 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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaplfold_result2.ps	Tue Dec 06 12:33:31 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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaplot_input1.dbn	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,3 @@
+>Anolis_carolinensis_chrUn_GL343590.trna2-AlaAGC (218800-218872)  Ala (AGC) 73 bp  Sc: 49.55
+UGGGAAUUAGCUCAAAUGGUAGAGCGCUCGCUUAGCAUGUGAGAGGUAGUGGGAUCGAUGCCCACAUUCUCCA
+.((((((..((((........)))).(((((.......))))).....(((((.......))))))))))).. (-21.90)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaplot_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnasnoop_input1a.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,2 @@
+>homo
+CGCGACCUCAGAUCAGACGUGGCGACCCGCUGAAUUUAAGCAUAUUAGUCAGCGGAGGAAAAGAAACUAACCAGGAUUCCCUCAGUAACGGCGAGUGAACAGGGAAGAGCCCAGCGCCGAAUCCCCGCCCCGCGGGGCGCGGGACAUGUGGCGUACGGAAGACCCGCUCCCCGGCGCCGCUCGUGGGGGGCCCAAGUCCUUCUGAUCGAGGCCCAGCCCGUGGACGGUGUGAGGCCGGUAGCGGCCGGCGCGCGCCCGGGUCUUCCCGGAGUCGGGUUGCUUGGGAAUGCAGCCCAAAGCGGGUGGUAAACUCCAUCUAAGGCUAAAUACCGGCACGAGACCGAUAGUCAACAAGUACCGUAAGGGAAAGUUGAAAAGAACUUUGAAGAGAGAGUUCAAGAGGGCGUGAAACCGUUAAGAGGUAAACGGGUGGGGUCCGCGCAGUCCGCCCGGAGGAUUCAACCCGGCGGCGGGUCCGGCCGUGUCGGCGGCCCGGCGGAUCUUUCCCGCCCCCCGUUCCUCCCGACCCCUCCACCCGCCCUCCCUUCCCCCGCCGCCCCUCCUCCUCCUCCCCGGAGGGGGCGGGCUCCGGCGGGUGCGGGGGUGGGCGGGCGGGGCCGGGGGUGGGGUCGGCGGGGGACCGUCCCCCGACCGGCGACCGGCCGCCGCCGGGCGCAUUUCCACCGCGGCGGUGCGCCGCGACCGGCUCCGGGACGGCUGGGAAGGCCCGGCGGGGAAGGUGGCUCGGGGGGCCCCGUCCGUCCGUCCGUCCUCCUCCUCCCCCGUCUCCGCCCCCCGGCCCCGCGUCCUCCCUCGGGAGGGCGCGCGGGUCGGGGCGGCGGCGGCGGCGGCGGUGGCGGCGGCGGCGGGGGCGGCGGGACCGAAACCCCCCCCGAGUGUUACAGCCCCCCCGGCAGCAGCACUCGCCGAAUCCCGGGGCCGAGGGAGCGAGACCCGUCGCCGCGCUCUCCCCCCUCCCGGCGCCCACCCCCGCGGGGAAUCCCCCGCGAGGGGGGUCUCCCCCGCGGGGGCGCGCCGGCGUCUCCUCGUGGGGGGGCCGGGCCACCCCUCCCACGGCGCGACCGCUCUCCCACCCCUCCUCCCCGCGCCCCCGCCCCGGCGACGGGGGGGGUGCCGCGCGCGGGUCGGGGGGCGGGGCGGACUGUCCCCAGUGCGCCCCGGGCGGGUCGCGCCGUCGGGCCCGGGGGAGGUUCUCUCGGGGCCACGCGCGCGUCCCCCGAAGAGGGGGACGGCGGAGCGAGCGCACGGGGUCGGCGGCGACGUCGGCUACCCACCCGACCCGUCUUGAAACACGGACCAAGGAGUCUAACACGUGCGCGAGUCGGGGGCUCGCACGAAAGCCGCCGUGGCGCAAUGAAGGUGAAGGCCGGCGCGCUCGCCGGCCGAGGUGGGAUCCCGAGGCCUCUCCAGUCCGCCGAGGGCGCACCACCGGCCCGUCUCGCCCGCCGCGCCGGGGAGGUGGAGCACGAGCGCACGUGUUAGGACCCGAAAGAUGGUGAACUAUGCCUGGGCAGGGCGAAGCCAGAGGAAACUCUGGUGGAGGUCCGUAGCGGUCCUGACGUGCAAAUCGGUCGUCCGACCUGGGUAUAGGGGCGAAAGACUAAUCGAACCAUCUAGUAGCUGGUUCCCUCCGAAGUUUCCCUCAGGAUAGCUGGCGCUCUCGCAGACCCGACGCACCCCCGCCACGCAGUUUUAUCCGGUAAAGCGAAUGAUUAGAGGUCUUGGGGCCGAAACGAUCUCAACCUAUUCUCAAACUUUAAAUGGGUAAGAAGCCCGGCUCGCUGGCGUGGAGCCGGGCGUGGAAUGCGAGUGCCUAGUGGGCCACUUUUGGUAAGCAGAACUGGCGCUGCGGGAUGAACCGAACGCCGGGUUAAGGCGCCCGAUGCCGACGCUCAUCAGACCCCAGAAAAGGUGUUGGUUGAUAUAGACAGCAGGACGGUGGCCAUGGAAGUCGGAAUCCGCUAAGGAGUGUGUAACAACUCACCUGCCGAAUCAACUAGCCCUGAAAAUGGAUGGCGCUGGAGCGUCGGGCCCAUACCCGGCCGUCGCCGGCAGUCGAGAGUGGACGGGAGCGGCGGGGGCGGCGCGCGCGCGCGCGCGUGUGGUGUGCGUCGGAGGGCGGCGGCGGCGGCGGCGGCGGGGGUGUGGGGUCCUUCCCCCGCCCCCCCCCCCACGCCUCCUCCCCUCCUCCCGCCCACGCCCCGCUCCCCGCCCCCGGAGCCCCGCGGACGCUACGCCGCGACGAGUAGGAGGGCCGCUGCGGUGAGCCUUGAAGCCUAGGGCGCGGGCCCGGGUGGAGCCGCCGCAGGUGCAGAUCUUGGUGGUAGUAGCAAAUAUUCAAACGAGAACUUUGAAGGCCGAAGUGGAGAAGGGUUCCAUGUGAACAGCAGUUGAACAUGGGUCAGUCGGUCCUGAGAGAUGGGCGAGCGCCGUUCCGAAGGGACGGGCGAUGGCCUCCGUUGCCCUCGGCCGAUCGAAAGGGAGUCGGGUUCAGAUCCCCGAAUCCGGAGUGGCGGAGAUGGGCGCCGCGAGGCGUCCAGUGCGGUAACGCGACCGAUCCCGGAGAAGCCGGCGGGAGCCCCGGGGAGAGUUCUCUUUUCUUUGUGAAGGGCAGGGCGCCCUGGAAUGGGUUCGCCCCGAGAGAGGGGCCCGUGCCUUGGAAAGCGUCGCGGUUCCGGCGGCGUCCGGUGAGCUCUCGCUGGCCCUUGAAAAUCCGGGGGAGAGGGUGUAAAUCUCGCGCCGGGCCGUACCCAUAUCCGCAGCAGGUCUCCAAGGUGAACAGCCUCUGGCAUGUUGGAACAAUGUAGGUAAGGGAAGUCGGCAAGCCGGAUCCGUAACUUCGGGAUAAGGAUUGGCUCUAAGGGCUGGGUCGGUCGGGCUGGGGCGCGAAGCGGGGCUGGGCGCGCGCCGCGGCUGGACGAGGCGCGCGCCCCCCCCACGCCCGGGGCACCCCCCUCGCGGCCCUCCCCCGCCCCACCCGCGCGCGCCGCUCGCUCCCUCCCCACCCCGCGCCCUCUCUCUCUCUCUCUCCCCCGCUCCCCGUCCUCCCCCCUCCCCGGGGGAGCGCCGCGUGGGGGCGCGGCGGGGGGAGAAGGGUCGGGGCGGCAGGGGCCGCGCGGCGGCCGCCGGGGCGGCCGGCGGGGGCAGGUCCCCGCGAGGGGGGCCCCGGGGACCCGGGGGGCCGGCGGCGGCGCGGACUCUGGACGCGAGCCGGGCCCUUCCCGUGGAUCGCCCCAGCUGCGGCGGGCGUCGCGGCCGCCCCCGGGGAGCCCGGCGGCGGCGCGGCGCGCCCCCCACCCCCACCCCACGUCUCGGUCGCGCGCGCGUCCGCUGGGGGCGGGAGCGGUCGGGCGGCGGCGGUCGGCGGGCGGCGGGGCGGGGCGGUUCGUCCCCCCGCCCUACCCCCCCGGCCCCGUCCGCCCCCCGUUCCCCCCUCCUCCUCGGCGCGCGGCGGCGGCGGCGGCAGGCGGCGGAGGGGCCGCGGGCCGGUCCCCCCCGCCGGGUCCGCCCCCGGGGCCGCGGUUCCGCGCGCGCCUCGCCUCGGCCGGCGCCUAGCAGCCGACUUAGAACUGGUGCGGACCAGGGGAAUCCGACUGUUUAAUUAAAACAAAGCAUCGCGAAGGCCCGCGGCGGGUGUUGACGCGAUGUGAUUUCUGCCCAGUGCUCUGAAUGUCAAAGUGAAGAAAUUCAAUGAAGCGCGGGUAAACGGCGGGAGUAACUAUGACUCUCUUAAGGUAGCCAAAUGCCUCGUCAUCUAAUUAGUGACGCGCAUGAAUGGAUGAACGAGAUUCCCACUGUCCCUACCUACUAUCCAGCGAAACCACAGCCAAGGGAACGGGCUUGGCGGAAUCAGCGGGGAAAGAAGACCCUGUUGAGCUUGACUCUAGUCUGGCACGGUGAAGAGACAUGAGAGGUGUAGAAUAAGUGGGAGGCCCCCGGCGCCCCCCCGGUGUCCCCGCGAGGGGCCCGGGGCGGGGUCCGCGGCCCUGCGGGCCGCCGGUGAAAUACCACUACUCUGAUCGUUUUUUCACUGACCCGGUGAGGCGGGGGGGCGAGCCCGAGGGGCUCUCGCUUCUGGCGCCAAGCGCCCGCCCGGCCGGGCGCGACCCGCUCCGGGGACAGUGCCAGGUGGGGAGUUUGACUGGGGCGGUACACCUGUCAAACGGUAACGCAGGUGUCCUAAGGCGAGCUCAGGGAGGACAGAAACCUCCCGUGGAGCAGAAGGGCAAAAGCUCGCUUGAUCUUGAUUUUCAGUACGAAUACAGACCGUGAAAGCGGGGCCUCACGAUCCUUCUGACCUUUUGGGUUUUAAGCAGGAGGUGUCAGAAAAGUUACCACAGGGAUAACUGGCUUGUGGCGGCCAAGCGUUCAUAGCGACGUCGCUUUUUGAUCCUUCGAUGUCGGCUCUUCCUAUCAUUGUGAAGCAGAAUUCGCCAAGCGUUGGAUUGUUCACCCACUAAUAGGGAACGUGAGCUGGGUUUAGACCGUCGUGAGACAGGUUAGUUUUACCCUACUGAUGAUGUGUUGUUGCCAUGGUAAUCCUGCUCAGUACGAGAGGAACCGCAGGUUCAGACAUUUGGUGUAUGUGCUUGGCUGAGGAGCCAAUGGGGCGAAGCUACCAUCUGUGGGAUUAUGACUGAACGCCUCUAAGUCAGAAUCCCGCCCAGGCGAACGAUACGGCAGCGCCGCGGAGCCUCGGUUGGCCUCGGAUAGCCGGUCCCCCGCCUGUCCCCGCCGGCGGGCCGCCCCCCCCUCCACGCGCCCCGCCGCGGGAGGGCGCGUGCCCCGCCGCGCGCCGGGACCGGGGUCCGGUGCGGAGUGCCCUUCGUCCUGGGAAACGGGGCGCGGCCGGAAAGGCGGCCGCCCCCUCGCCCGUCACGCACCGCACGUUCGUGGGGAACCUGGCGCUAAACCAUUCGUAGACGACCUGCUUCUGGGUCGGGGUUUCGUACGUAGCAGAGCAGCUCCCUCGCUGCGAUCUAUUGAAAGUCAGCCCUCGACACAAGGGUUUGUC
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnasnoop_input1b.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -0,0 +1,2 @@
+>ACA51
+AACCTACCCCATATACACCTCAGCTCAGGCCCTGTGCCTGGTCTGTATTGTGAATGGGGGAACATAG
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnasnoop_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnasubopt_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnasubopt_result1.txt	Tue Dec 06 12:33:31 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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaup_input1.fa	Tue Dec 06 12:33:31 2016 -0500
@@ -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
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaup_result1.txt	Tue Dec 06 12:33:31 2016 -0500
@@ -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
Binary file vienna_rna.tar.gz has changed