Mercurial > repos > rnateam > viennarna_rnaalifold

changeset 0:c8fff1159bbb 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:31:07 -0500
parents
children 10b56bba8dee
files macros.xml rnaalifold.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, 1021 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:31:07 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/rnaalifold.xml	Tue Dec 06 12:31:07 2016 -0500
@@ -0,0 +1,243 @@
+<tool id="viennarna_rnaalifold" name="@EXECUTABLE@" version="@VERSION@.0">
+    <description>Calculate minimum free energy secondary structures and partition function on a multiple alignment file</description>
+    <macros>
+        <token name="@EXECUTABLE@">RNAalifold</token>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements" />
+    <expand macro="stdio" />
+    <expand macro="version_command" />
+    <command>
+<![CDATA[
+        
+        RNAalifold  < '$input' > '$tabularFile'
+        -T$model_options.temperature -d$model_options.dangling
+        $general_options.verbose
+        $general_options.color
+        $general_options.alignment
+        #if $general_options.layout_type ==0
+            --layout-type=$general_options.layout_type
+        #end if
+        $general_options.noPS
+        #if str($IDs.id_prefix) <> "alignment"
+            --id-prefix='$IDs.id_prefix'
+        #end if
+        $IDs.continuous_ids
+        $IDs.auto_id
+        #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 $constraints.maxBPspan <> -1
+            --maxBPspan=$constraints.maxBPspan
+        #end if
+        #if str($constraints.constraintLocation.constraintSelector) == "fromFile"
+            --constraint='$constraints.constraintLocation.constraintsFile'
+            $constraints.constraintLocation.batch
+            $constraints.ss_cons
+        #end if        
+        #if str($constraints.shapeOption.shapeSelector) == "isUsed"
+            #if len($constraints.shapeOption.shape_files) > 0
+                #set $sf=list()
+                #for $i in $constraints.shapeOption.shape_files
+                    $sf.append(str('$i.shape'))
+                #end for
+                --shape= #echo ','.join($sf)
+                #if $constraints.shapeOption.m <> 1.8 or  $constraints.shapeOption.b <> -0.6
+                    #set $s="Dm"+str($constraints.shapeOption.m)+"b"+str($constraints.shapeOption.b)
+                    --shapeMethod=$s
+                #end if		
+            #end if
+        #end if		
+        #if $algorithm_options.pf <> "-1"
+            --partfunc=$algorithm_options.pf
+        #end if
+        #if $algorithm_options.mea <> 1.0
+            --mea=$algorithm_options.mea
+        #end if
+        $algorithm_options.mis
+        #if $algorithm_options.stochBT_en <> 1
+            --stochBT_en=$algorithm_options.stochBT_en
+        #end if
+        #if $algorithm_options.pfScale <> 1.07
+            --pfScale=$algorithm_options.pfScale
+        #end if
+        $algorithm_options.circular
+        #if $algorithm_options.bppmThreshold <> 1e-6
+            --bppmThreshold=$algorithm_options.bppmThreshold
+        #end if
+        $algorithm_options.gquad
+        $algorithm_options.sci		
+        $model_options.notetra
+        $model_options.nolp
+        $model_options.nogu
+        $model_options.noclosinggu
+        #if $model_options.cfactor <> 1.0
+            --cfactor=$model_options.cfactor
+        #end if
+        #if $model_options.nfactor <> 1.0
+            --nfactor=$model_options.nfactor
+        #end if
+        $model_options.endgaps
+        $model_options.ribosum
+        #if $model_options.nsp
+            --nsp='$model_options.nsp'
+        #end if
+        #if $model_options.betaScale <> 1.0
+            --betaScale=$model_options.betaScale
+        #end if
+        && tar -cf '$imagesFile' *.ps		
+]]>
+    </command>
+    <inputs>
+        <param format="txt" name="input" type="data" label="Clustal file"/>        
+        <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="mea" type="float" value="1.0" label="Gamma Value" help="Calculate an MEA (maximum expected accuracy) structure. The parameter gamma tunes the importance of correctly predicted pairs versus unpaired bases. Thus, for small values of gamma the MEA structure will contain only pairs with very high probability. where the expected accuracy is computed from the pair probabilities: each base pair (i,j) gets a score 2*gamma*p_ij and the score of an unpaired base is given by the probability of not forming a pair." argumen="--MEA"/>
+            <param argument="--mis" type="boolean" checked="false" truevalue="--mis" falsevalue="" label="Most Informative Sequence" help="Output 'most informative sequence' instead of simple consensu"/>							
+            <param name="stochBT_en" type="integer" value="1" label="Number of random structures." help="Compute a certain number of random structures with a probability dependend on the partition function. Print out the energies and probabilities of the backtraced structures. " argument="--stochBT_en"/>	        
+            <param name="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 name="circular" type="boolean" truevalue="--circ" falsevalue="" checked="false" label="Assume circular RNA structure" help="--circ"/>
+            <param name="bppmThreshold" type="float" value="1e-6" 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 name="gquad" type="boolean" truevalue="--gquad" falsevalue="" checked="false" label="G Quadruplex formation" help="Take into account G Quadruplex formation (-g)"/>
+            <param name="sci" type="boolean" truevalue="--sci" falsevalue="" checked="false" label="Structure conservation index" 
+                help="Compute the structure conservation index (SCI) for the MFE consensus structure of the alignment" argument="--sci"/>
+        </section>    
+        <section name="general_options" title="General Options">
+            <param name="verbose" type="boolean" truevalue="--v" falsevalue="" checked="false" label=" Be verbose" argument="--v"/>
+            <param name="color" type="boolean" truevalue="--color" falsevalue="" checked="false" label="Colored secondary structure image" 
+                help="Produce a colored version of the consensus strcture plot" argument="--color"/>
+            <param name="alignment" type="boolean" truevalue="--aln" falsevalue="" checked="false" 
+                label="Colored and annotated alignment image" 
+                help="Produce a colored and structure annotated alignment" 
+                argument="--aln"/>
+            <param name="layout_type" type="select" label="Layout algorithm" argument="--layout-type">
+                <option value="1" selected="true">Default: Naview layout</option>
+                <option value="0">Simple radial layout</option>
+            </param>
+            <param name="noPS" type="boolean" truevalue="" falsevalue="--noPS" checked="true" label="Produce postscript output" help="" argument="--noPS"/>
+        </section>
+        <section name="IDs" title="Naming Conventions">		
+                <param name="continuous_ids" 
+                    type="boolean" truevalue="--continuous-ids" falsevalue="" checked="false" 
+                    label="Continuous alignment ID numbering" 
+                    help="Use continuous alignment ID numbering when no alignment ID can be retrieved from input data." 
+                    argument="--continuous-ids"/>
+                <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="constraints" title="Structure constraints">
+                <param name="maxBPspan" type="integer" value="-1" label="Set the maximum base pair span" help="" argument="--maxBPspan"/>
+                <conditional name="constraintLocation">
+                    <param name="constraintSelector" type="select" label="Constraints" argument="--constraint">
+                        <!--	<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"/>
+                        <param name="batch" type="boolean" checked="false" truevalue="--batch" falsevalue="" 
+                            label="Use constraints for all alignment records" 
+                            help="Usually, constraints provided from input file are only applied to a single sequence alignment. Therefore, RNAalifold will stop its computation and quit after the first input alignment was processed. Using this switch, RNAalifold processes all sequence alignments in the input and applies the same provided constraints to each of them." 
+                            argument="--batch"/>
+                        <param name="ss_cons" type="boolean" checked="false" truevalue="--SS_cons" falsevalue="" 
+                            label="Use consensus structures from Stockholm file (#=GF SS_cons) as constraint" 
+                            help="Stockholm formatted alignment files have the possibility to store a secondary structure string in one of if ('#=GC') column annotation meta tags. The corresponding tag name is usually 'SS_cons', a consensus secondary structure. Activating this flag allows one to use this consensus secondary structure from the input file as structure constraint. Note: Activating this flag implies −−constraint." 
+                            argument="--SS_cons"/>	
+                    </when>
+                </conditional>
+                <conditional name="shapeOption">
+                        <param name="shapeSelector" type="select" label="Shape reactivity data">
+                            <option value="isUsed">Use shape reactivity data</option>
+                            <option value="notUsed" selected="true">Don't use shape reactivity data</option>
+                        </param>
+                        <when value ="isUsed">
+                            <repeat name="shape_files" title="Shape File" default="0" argument="--shape">
+                                <param type="data" name="shape" format="shape,*" label="Shape file"/>
+                            </repeat>
+                                <param name="m" type="float" value="1.8" label="Slope m" />
+                                <param name="b" type="float" value="-0.6" label="Intercept b" help="Currently, the only data conversion method available is that of to Deigan et al 2009. This method is the default. The slope 'm' and the intercept 'b' can be set to a non−default value if necessary. Otherwise m=1.8 and b=−0.6 as stated in the paper mentioned before."/>
+                        </when>
+                        <when value="notUsed">
+                        </when>
+                </conditional>
+        </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="2" selected="True">2: unpaired bases participate in all dangling ends</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="cfactor" type="float" value="1.0" label="Weight of the covariance term" help="--cfactor"/>
+            <param name="nfactor" type="float" value="1.0" label="Penalty for non-compatible sequences in the covariance term" help="--nfactor"/>
+            <param name="endgaps" type="boolean" truevalue="--endgaps" falsevalue="" checked="false" label="Score pairs with endgaps same as gap-gap pairs" help="--endgaps"/>
+            <param name="ribosum" type="boolean" truevalue="--ribosum_scoring" falsevalue="" checked="false" label="Use ribosum scoring matrix" help="--ribosum_scoring"/>
+            <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>
+    </inputs>
+    <outputs>
+        <data format="txt" name="tabularFile"/>
+        <data format="tar" name="imagesFile"/>
+    </outputs>
+    <tests>
+        <test>
+            <param name="input" value="rnaalifold_input1.clustal"/>
+            <output name="output" file="rnaalifold_result1.txt"/>
+        </test>
+    </tests>
+    <help>
+<![CDATA[
+**RNAalifold**
+
+-----
+
+**Input format**
+
+RNAalifold requires one input file
+- Clustal file
+
+------
+
+**Outputs**
+
+- energy of the consensus structures in the clustal file
+- several possible postscript images bundled together in a tar file
+
+
+]]>
+    </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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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:31:07 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