Mercurial > repos > rnateam > locarna_reliability_profile
changeset 3:e84d1fc87829 draft
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/locarna commit 10ccc47885ce71e602d66e157bd475f1facbd042
| author | bgruening |
|---|---|
| date | Mon, 05 Dec 2022 14:29:00 +0000 |
| parents | dd991ca3d8ce |
| children | 7d63908343c3 |
| files | locarna.tar.bz2 locarna_reliability_profile.xml macros.xml test-data/archaea-default.stdout test-data/archaea_relplot.scr test-data/haca.snoRNA-c.aln test-data/haca.snoRNA-default.stdout test-data/haca.snoRNA.aln test-data/tRNA_2.aln test-data/tRNA_2.epms |
| diffstat | 10 files changed, 116 insertions(+), 149 deletions(-) [+] |
line wrap: on
line diff
--- a/locarna_reliability_profile.xml Sat Jan 21 17:38:37 2017 -0500 +++ b/locarna_reliability_profile.xml Mon Dec 05 14:29:00 2022 +0000 @@ -1,4 +1,4 @@ -<tool id="locarna_reliability_profile" name="LocARNA reliability-profile" version="@VERSION@.1"> +<tool id="locarna_reliability_profile" name="LocARNA reliability-profile" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@"> <description> Plot reliability profile for probabilistic mlocarna </description> @@ -76,42 +76,38 @@ or leave blank to plot without projection." /> - <param name="predict" type="boolean" checked="True" + <param argument="--dont-predict" name="predict" type="boolean" checked="true" truevalue="" falsevalue="--dont-predict" - argument="dont-predict" help="Turn on/off predicting." /> - <param name="fit_once_on" type="boolean" checked="False" + <param argument="--fit-once-on" type="boolean" checked="false" truevalue="--fit-once-on" falsevalue="" - argument="fit-once-on" help="Predict exactly one 'on'-region" /> - <param name="show_sw" type="boolean" checked="True" + <param argument="--show-sw" type="boolean" checked="true" truevalue="--show-sw" falsevalue="" - argument="show-sw" help="Show the influence of structure weight." /> - <param name="rev_compl" type="boolean" checked="False" + <param argument="--revcompl" name="rev_compl" type="boolean" checked="false" truevalue="--revcompl" falsevalue="" - argument="revcompl" help="Plot and fit a reverse complement." /> - <param argument="beta" help="Inverse temperature for fitting." + <param argument="--beta" help="Inverse temperature for fitting." type="float" value="12.0" min="1.0" max="25.0" /> - <param name="structure_weight" argument="structure-weight" + <param argument="--structure-weight" type="float" value="1.0" min="0.0" max="8.0" help="Weight of structure against sequence (1.0 means equal)" /> - <param argument="offset" help="Genomic offset." + <param argument="--offset" help="Genomic offset." type="integer" value="1" /> - <param name="signals" type="text" label="List of signals" + <param argument="--signals" type="text" label="List of signals" help="List of '{from to} orientation' signal specifications. Show signals in plot and compared infered signal to them. Ranges 'from to' can be repeated; the orientation must be either +1 or -1. @@ -125,18 +121,18 @@ <param name="output_format" type="select" label="Output options"> - <option value="pdf" selected="True">Plot in pdf + <option value="pdf" selected="true">Plot in pdf format</option> - <option value="png" selected="False">Plot in png format</option> + <option value="png" selected="false">Plot in png format</option> </param> - <param name="write_rscript" type="boolean" argument="write-R-script" - label="Write R script" checked="False" + <param name="write_rscript" type="boolean" argument="--write-R-script" + label="Write R script" checked="false" help="Write R script that can produce the output plot." /> - <param name="write_subsequence" type="boolean" argument="write-subseq" - label="Write Subsequence" checked="False" + <param name="write_subsequence" type="boolean" argument="--write-subseq" + label="Write Subsequence" checked="false" help="Write subsequence in output." /> @@ -145,25 +141,37 @@ <outputs> <data format="pdf" name="pdf_output_plot" - label="${tool.name} reliability plot (pdf) ${on_string}"> + label="${tool.name} on ${on_string}: plot"> <filter>output_format == 'pdf'</filter> </data> <data format="png" name="png_output_plot" - label="${tool.name} reliability plot (png) ${on_string}"> + label="${tool.name} on ${on_string}: plot"> <filter>output_format == 'png'</filter> </data> <data format="txt" name="rscript_output" - label="${tool.name} reliability plot (R script) ${on_string}"> + label="${tool.name} on ${on_string}: R script "> <filter> write_rscript is True </filter> </data> - <data format="txt" name="stdout" label="${tool.name} std out on ${on_string}" /> + <data format="txt" name="stdout" label="${tool.name} on ${on_string}: stdout" /> </outputs> <tests> - <test> + <test expect_num_outputs="3"> <param name="input_archive" value="archaea.tar.gz" /> - <param name="write_rscript" value="True" /> + <param name="write_rscript" value="true" /> + <output name="pdf_output_plot"> + <assert_contents> + <has_text text="%PDF"/> + <has_size value="5401" delta="100"/> + </assert_contents> + </output> <output name="rscript_output" file="archaea_relplot.scr" lines_diff="2" /> + <output name="stdout"> + <assert_contents> + <has_text text="ONOFF"/> + <has_text text="FIT"/> + </assert_contents> + </output> </test> </tests>
--- a/macros.xml Sat Jan 21 17:38:37 2017 -0500 +++ b/macros.xml Mon Dec 05 14:29:00 2022 +0000 @@ -1,9 +1,12 @@ <macros> - <token name="@VERSION@">1.9.0</token> + <token name="@TOOL_VERSION@">1.9.2.3</token> + <token name="@VERSION_SUFFIX@">1</token> + <token name="@PROFILE@">21.09</token> <xml name="requirements"> <requirements> - <requirement type="package" version="@VERSION@">locarna</requirement> + <requirement type="package" version="@TOOL_VERSION@">locarna</requirement> + <requirement type="package" version="4.2.2">r-base</requirement> </requirements> </xml> @@ -46,15 +49,15 @@ </xml> <xml name="common_scoring_parameters"> - <param name="struct_weight" argument="struct-weight" + <param argument="--struct-weight" label="Structure weight" type="integer" value="200" min="0" max="800" /> - <param name="indel_opening" argument="indel-opening" + <param argument="--indel-opening" label="Indel opening score" type="integer" value="-500" max="0" min="-1500" /> - <param argument="indel" label="Indel score" type="integer" + <param argument="--indel" label="Indel score" type="integer" value="-350" min="-1000" max="0" /> - <param argument="tau" type="integer" value="50" + <param argument="--tau" type="integer" value="50" min="0" max="200" label="Sequence contribution at structure match in percent"/> @@ -78,17 +81,17 @@ </xml> <xml name="plfolding_parameters"> - <param name="plfold_span" argument="--plfold-span" + <param argument="--plfold-span" type="integer" value="150" min="-1" max="400" label="Maximum basepair span by RNAplfold (local folding); -1 for global folding" /> - <param name="plfold_winsize" argument="--plfold-winsize" + <param argument="--plfold-winsize" type="integer" value="300" min="-1" max="800" label="Window size for local folding" /> </xml> <xml name="common_folding_parameters"> - <param name="rnafold_temperature" argument="rnafold-temperature" + <param argument="rnafold-temperature" type="float" value="37.0" min="10" max="50" label="Temperature for RNAfold (RNAfold's -T option)" /> </xml> @@ -147,16 +150,18 @@ </token> - <token name="@CONSTRAINT_ARGS@"> - ## -------------------- constraint parameters - $Constraints.lonely_pairs - + <token name="@CONSTRAINT_ARGS_WOLP@"> #if $Constraints.maxBPspan != -1 --maxBPspan $Constraints.maxBPspan #end if $Constraints.ignore_constraints </token> + <token name="@CONSTRAINT_ARGS@"> + ## -------------------- constraint parameters + $Constraints.lonely_pairs + @CONSTRAINT_ARGS_WOLP@ + </token> <token name="@STDOUT_ARGS@"> $stdout_verbosity @@ -171,7 +176,7 @@ label="Restrict alignable positions by maximum difference" help="(max-diff*)"> <option value="off">Off</option> - <option value="max-diff" selected="True">Maximal difference of aligned positions</option> + <option value="max-diff" selected="true">Maximal difference of aligned positions</option> <option value="max-diff-at-am">Maximal difference of aligned positions at arc matches</option> <yield /> </param> @@ -180,12 +185,12 @@ <xml name="max_diff_parameters_standard_cases"> <when value="off" /> <when value="max-diff"> - <param name="max_diff" argument="max-diff" type="integer" + <param argument="--max-diff" type="integer" value="60" min="-1" max="300" label="Maximal difference of aligned positions" /> </when> <when value="max-diff-at-am"> - <param name="max_diff_at_am" argument="max-diff-at-am" type="integer" + <param argument="--max-diff-at-am" type="integer" value="60" min="-1" max="300" label="Maximal difference of aligned positions, only at arc match positions" /> @@ -206,13 +211,13 @@ </expand> <expand macro="max_diff_parameters_standard_cases" /> <when value="max-diff-aln"> - <param name="max_diff_aln" argument="max-diff-aln" type="data" format="clustal" + <param argument="--max-diff-aln" type="data" format="clustal" label="Reference alignment" /> - <param name="max_diff" argument="max-diff" type="integer" + <param argument="--max-diff" type="integer" value="60" min="-1" max="300" label="Maximal difference of aligned positions" /> - <param name="max_diff_relax" argument="max-diff-relax" type="boolean" + <param argument="--max-diff-relax" type="boolean" truevalue="--max-diff-relax" falsevalue="" checked="False" label="Relax deviation constraints in multiple @@ -223,75 +228,79 @@ </xml> <xml name="common_heuristic_parameters"> - <param name="min_prob" argument="min-prob" type="float" value="0.0005" + <param argument="--min-prob" type="float" value="0.0005" min="0.0" max="0.2" label="Minimal / cutoff probability" /> - <param name="max_diff_am" argument="max-diff-am" + <param argument="--max-diff-am" type="integer" value="30" min="-1" max="300" label="Maximal difference for sizes of matched arcs (-1=off)" /> - <param name="max_bps_length_ratio" argument="max-bps-length-ratio" + <param argument="--max-bps-length-ratio" type="float" value="0.0" min="0.0" max="10.0" label="Maximal ratio 'number of considered base pairs' by 'sequence length' (default: 0.0 = no effect)" /> </xml> <xml name="common_heuristic_parameters_sparse"> - <param name="prob_unpaired_in_loop_threshold" argument="prob-unpaired-in-loop-threshold" type="float" value="0.00005" + <param argument="--prob-unpaired-in-loop-threshold" type="float" value="0.00005" min="0.0" max="0.2" label="Probability threshold for unpaired bases in loops" /> - <param name="prob_basepair_in_loop_threshold" argument="prob-basepair-in-loop-threshold" type="float" value="0.0001" + <param argument="--prob-basepair-in-loop-threshold" type="float" value="0.0001" min="0.0" max="0.2" label="Probability threshold for base pairs in loops" /> </xml> <xml name="in_loop_ratio_parameters"> - <param name="max_uil_length_ratio" argument="max-uil-length-ratio" + <param argument="--max-uil-length-ratio" type="float" value="0.0" min="0.0" max="10.0" label="Maximal ratio 'number considered unpaired bases in loops' by 'sequence length' (default: 0.0; no effect)" /> - <param name="max_bpil_length_ratio" argument="max-bpil-length-ratio" + <param argument="--max-bpil-length-ratio" type="float" value="0.0" min="0.0" max="10.0" label="Maximal ratio 'number of considered base pairs in loops' by 'sequence length' (default: 0.0; no effect)" /> </xml> <xml name="exparnap_in_loop_parameters"> - <param name="prob_unpaired_in_loop_threshold" - argument="prob_unpaired_in_loop_threshold" + <param argument="--prob_unpaired_in_loop_threshold" type="float" value="0.01" min="0.0" max="1.0" label="Probability threshold for unpaired bases in loops" /> - <param name="prob_basepair_in_loop_threshold" - argument="prob_basepair_in_loop_threshold" + <param argument="--prob_basepair_in_loop_threshold" type="float" value="0.01" min="0.0" max="1.0" label="Probability threshold for base pairs in loops" /> </xml> <xml name="alifold_consensus_parameter"> - <param name="alifold_consensus_dp" argument="alifold-consensus-dp" + <param argument="--alifold-consensus-dp" type="boolean" checked="False" truevalue="--alifold-consensus-dp" falsevalue="" label="Compute consensus dot plot by alifold" /> </xml> - <xml name="common_constraint_parameters"> - <param name="lonely_pairs" type="boolean" truevalue="--LP" falsevalue="--noLP" - checked="false" label="Allow lonely base-pairs" help="(--LP/--noLP)" /> - <param name="maxBPspan" argument="--maxBPspan" + <xml name="common_constraint_parameters_woLP"> + <yield/> + <param argument="--maxBPspan" type="integer" value="-1" min="-1" max="400" label="Maximum basepair span by RNAfold; -1 for arbitrary span" /> - <param name="ignore_constraints" argument="ignore-constraints" + <param argument="--ignore-constraints" type="boolean" checked="false" truevalue="--ignore-constraints" falsevalue="" help="Ignore all anchor and structure constraints given in the fasta(-ish) input." /> </xml> + <xml name="common_constraint_parameters"> + <expand macro="common_constraint_parameters_woLP"> + <param argument="--noLP" name="lonely_pairs" type="boolean" truevalue="" falsevalue="--noLP" + checked="true" label="Allow lonely base-pairs"/> + </expand> + </xml> + <xml name="standard_outupt"> - <data format="txt" name="stdout" label="${tool.name} std out on ${on_string}"> + <data format="txt" name="stdout" label="${tool.name} on ${on_string}: stdout"> <filter>stdout_verbosity != '--quiet'</filter> </data> </xml> @@ -299,25 +308,25 @@ <xml name="mlocarna_outputs"> <expand macro="standard_outupt" /> <data format="clustal" name="clustal" from_work_dir="mlocarna_results/results/result.aln" - label="${tool.name} alignment (annotated clustal) on ${on_string}"> + label="${tool.name} on ${on_string}: alignment (annotated clustal)"> <filter>'clustal' in outputs</filter> </data> <data format="clustal" name="clustal_strict" from_work_dir="mlocarna_results/results/result.strict-aln" - label="${tool.name} alignment (clustal) on ${on_string}"> - <filter>'clustal_strict' in outputs</filter> + label="${tool.name} on ${on_string}: alignment (clustal)"> + <filter>'strict' in outputs</filter> </data> <data format="stockholm" name="stockholm" from_work_dir="mlocarna_results/results/result.stk" - label="${tool.name} alignment (stockholm) on ${on_string}"> + label="${tool.name} on ${on_string}: alignment (stockholm)"> <filter>'stockholm' in outputs</filter> </data> <data format="txt" name="pp" from_work_dir="mlocarna_results/results/result.pp" - label="${tool.name} alignment (PP 2.0) on ${on_string}"> + label="${tool.name} on ${on_string}: alignment (PP 2.0)"> <filter>'pp' in outputs</filter> </data> <data format="tar.gz" name="mlocarna_results_tgz" - label="${tool.name} results archive on ${on_string}"> + label="${tool.name} on ${on_string}: results archive"> <filter>'mlocarna_results' in outputs</filter> </data> </xml> @@ -330,8 +339,5 @@ <citation type="doi">10.1186/s12859-014-0404-0</citation> </citations> </xml> - - - </macros>
--- a/test-data/archaea-default.stdout Sat Jan 21 17:38:37 2017 -0500 +++ b/test-data/archaea-default.stdout Mon Dec 05 14:29:00 2022 +0000 @@ -6,13 +6,4 @@ Perform progressive alignment ... - -vhuU AG-CUCACAACCGAACCC-AU------------UUGGGAGGUUGUGAGCU- -fwdB AU-GUUGGAGGGGAACCC-GU------------AAGGGACCCUCCAAGAU- -selD UUACGAUGUGCCGAACCCUUU------------AAGGGAGGCACAUCGAAA -hdrA GG--CACCACUCGAAGGC--U------------AAGCCAAAGUGGUG-CU- -vhuD GU--UCUCUCGGGAACCCGUC------------AAGGGACCGAGAGA-AC- -fruA ---CCUCGAGGGGAACCC-GA------------AAGGGACCCGAGAGG--- -fdhA CG-CCACCCUGCGAACCCAAUAUAAAAUAAUACAAGGGAGCAG-GUGGCG- - -alifold ((.(((((((((...(((.................))).))))))))))). (-31.59 = -20.01 + -11.58) +alifold ((.(((((((((...(((.................))).))))))))))). (-31.36 = -20.01 + -11.34)
--- a/test-data/archaea_relplot.scr Sat Jan 21 17:38:37 2017 -0500 +++ b/test-data/archaea_relplot.scr Mon Dec 05 14:29:00 2022 +0000 @@ -48,7 +48,7 @@ par(mar=c(6,2.5,1,1)) # open plot (and draw threshold) -plot(c(0),c(0),type="l", +plot(c(0),c(0),type="l", xlab="",ylab="", xlim=the_xlim,ylim=c(0,maxy), yaxp=c(0,1,2)) @@ -87,15 +87,15 @@ colors<-c(colors,colors); if (0>0) { - + for (i in 1:0) { orientation <- signals[signal_starts[i]+signal_sizes[i]*2]; sig_y <- maxy-i*anno_space; - + for (j in 0:(signal_sizes[i]-1)) { sig_x <- c(signals[signal_starts[i]+j*2],signals[signal_starts[i]+j*2+1]); - + ## draw arrows if (orientation!=0) { the_code <- 1+(orientation+1)/2; @@ -109,9 +109,9 @@ #draw inferred on-signal hit_color <- rgb(0.1,0.6,0.1,0.9) - + if (0!=1) { - + on <- c(0,4,33,44,49); off <- c(3,20,43,48,51);
--- a/test-data/haca.snoRNA-c.aln Sat Jan 21 17:38:37 2017 -0500 +++ b/test-data/haca.snoRNA-c.aln Mon Dec 05 14:29:00 2022 +0000 @@ -1,22 +1,15 @@ -CLUSTAL W --- LocARNA 1.9.0 - -ACA7 ACCUCCUGGGAUCGCAUCUGG-AGAGUGC---CUAGUAUUCUGCCAGCUU-CGGAAAGGG -ACA30 UGGCACUUU--CACAGUUCCU-UCCCCAG---GCAGUGGGGCCAGGAUUUGGUAGCUGGU -ACA5 UGCAGCCGUGUCAAAUUCAGUACCUGUCCUAUGCAUGGUAGGCACUGGC--CCAGAAGGC -ACA59 GCCCAGGGUAUGUUCACGGGGCGAUGCUGCCCUCCCAGCUGG-CCCAUG--GGUGACCCU -#A1 ............................................................ -#A2 ............................................................ +CLUSTAL W --- LocARNA 1.9.2.3 -ACA7 AGGGAAAGCAAGCCUGGCAGAG-GCACCCAUUCCAUUCCCAGCUUGCUCCGUAGCUGGCG -ACA30 GCUGAGAGAAAAC---CC-UUG---AUUGUAUUCUUGCCCUGGG---AUUAUACCAGUGG -ACA5 UGCCACAGAAACAC--UGUGAC-UCAUGG-----GCCCUGUUCCUGUGUCCCAGGCUCAG -ACA59 GGGAACAUUAACUGCCUCACAACGUUUGUGCCUCAGUUACCCGUAGAUGUAGUGAGGGUA -#A1 ....AAAAAA.................................................. -#A2 ....123456.................................................. +ACA7 ACCUCCUGGGAUCGCAUCUGG-AGAGUGC---CUAGUAUUCUGCCAGCUU-CGGAAAGGGAGGGAAAGCAAGCCUGGCAGAGGCACCCAUUC-CAUUCCCAGCUUGCUCCGUAGCUGGCG +ACA30 UGGCACUUU--CACAGUUCCU-UCCCCAG---GCAGUGGGGCCAGGAUUUGGUAGCUGGUGCUGAGAGAAAAC---CC-UUG--AUUGUAUU-CUUGCCCUGGG---AUUAUACCAGUGG +ACA5 UGCAGCCGUGUCAAAUUCAGUACCUGUCCUAUGCAUGGUAGGCACUGGC--CCAGAAGGCUGCCACAGAAACAC--UGUGACUCAUGG------GCCCUGUUCCUGUGUCCCAGGCUCAG +ACA59 GCCCAGGGUAUGUUCACGGGGCGAUGCUGCCCUCCCAGC-UGGCCCAUG--GGUGACCCUGGGAACAUUAACUGCCUCACAACGUUUGUGCCUCAGUUACCCGUAGAUGUAGUGAGGGUA +#A1 ................................................................AAAAAA.................................................. +#A2 ................................................................123456.................................................. -ACA7 AUUGGA--AGA---CACU-CUGCG-----ACA -ACA30 CAACUG--UCA---CUCA-AUGGG-----ACA -ACA5 GGAUAA--AUU---UGGU-UACAG-----ACA +ACA7 AUUGGAA-----GACACUC-UGC-----GACA +ACA30 CAACUGU-----CACUCAA-UGG-----GACA +ACA5 GGAUAAA-----UUUGGUU-ACA-----GACA ACA59 ACAAUACUUACUCUCGUUGGUGAUAAGGAACA #A1 .............................BBB #A2 .............................123
--- a/test-data/haca.snoRNA-default.stdout Sat Jan 21 17:38:37 2017 -0500 +++ b/test-data/haca.snoRNA-default.stdout Mon Dec 05 14:29:00 2022 +0000 @@ -6,28 +6,5 @@ Perform progressive alignment ... - -ACA7 ACCUCCUGGGAUCGCAUCUGG-AGAGUGC---CUAGUAUUCUGCCAGCUUCGGAA-AGGG -ACA30 UGGCACUUUCACAG--UUCCU-UCCCCAG---GCAGUGGGGCCAGGAUUUGGUAGCUGGU -ACA5 UGCAGCCGUGUCAAAUUCAGUACCUGUCCUAUGCAUGGUAGGCACUGGCCCAGAA--GGC -ACA59 GCCCAGGGUAUGUUCACGGGGCGAUGCUGCCCUCCCAGCUGG-CCCAUGGGUGAC--CCU -#A1 ............................................................ -#A2 ............................................................ - -ACA7 AGGGAAAGCAAGCCUGGCAGAG-GCACCCAUUCCAUUCCCAGCUUGCUCCGUAGCUGGCG -ACA30 GCUGAGAGAAAAC---CC-UUG---AUUGUAUUCUUGCCCUGGG---AUUAUACCAGUGG -ACA5 UGCCACAGAAACAC--UGUGAC-UCAUGG-----GCCCUGUUCCUGUGUCCCAGGCUCAG -ACA59 GGGAACAUUAACUGCCUCACAACGUUUGUGCCUCAGUUACCCGUAGAUGUAGUGAGGGUA -#A1 ....AAAAAA.................................................. -#A2 ....123456.................................................. - -ACA7 AUUGGA--AGA---CACU-CUGCG-----ACA -ACA30 CAACUG--UCA---CUCA-AUGGG-----ACA -ACA5 GGAUAA--AUU---UGGU-UACAG-----ACA -ACA59 ACAAUACUUACUCUCGUUGGUGAUAAGGAACA -#A1 .............................BBB -#A2 .............................123 - -alifold .((((((((........((((.(((((((......))))))).))))......))..))) - ))).............((((((.....((.....((((((((((((.....))))))))) - )))......)).....)).))))......... (-61.69 = -33.85 + -27.84) +alifold .((((((((........((((.(((((((......))))))).))))......))..)))))).............(((((((...((......((((((((((((.....))))))))) + ))).........)).))).))))......... (-62.29 = -34.50 + -27.79)
--- a/test-data/haca.snoRNA.aln Sat Jan 21 17:38:37 2017 -0500 +++ b/test-data/haca.snoRNA.aln Mon Dec 05 14:29:00 2022 +0000 @@ -1,16 +1,11 @@ -CLUSTAL W --- LocARNA 1.9.0 - -ACA7 ACCUCCUGGGAUCGCAUCUGG-AGAGUGC---CUAGUAUUCUGCCAGCUU-CGGAAAGGG -ACA30 UGGCACUUU--CACAGUUCCU-UCCCCAG---GCAGUGGGGCCAGGAUUUGGUAGCUGGU -ACA5 UGCAGCCGUGUCAAAUUCAGUACCUGUCCUAUGCAUGGUAGGCACUGGC--CCAGAAGGC -ACA59 GCCCAGGGUAUGUUCACGGGGCGAUGCUGCCCUCCCAGCUGG-CCCAUG--GGUGACCCU +CLUSTAL W --- LocARNA 1.9.2.3 -ACA7 AGGGAAAGCAAGCCUGGCAGAG-GCACCCAUUCCAUUCCCAGCUUGCUCCGUAGCUGGCG -ACA30 GCUGAGAGAAAAC---CC-UUG---AUUGUAUUCUUGCCCUGGG---AUUAUACCAGUGG -ACA5 UGCCACAGAAACAC--UGUGAC-UCAUGG-----GCCCUGUUCCUGUGUCCCAGGCUCAG -ACA59 GGGAACAUUAACUGCCUCACAACGUUUGUGCCUCAGUUACCCGUAGAUGUAGUGAGGGUA +ACA7 ACCUCCUGGGAUCGCAUCUGG-AGAGUGC---CUAGUAUUCUGCCAGCUU-CGGAAAGGGAGGGAAAGCAAGCCUGGCAGAGGCACCCAUUC-CAUUCCCAGCUUGCUCCGUAGCUGGCG +ACA30 UGGCACUUU--CACAGUUCCU-UCCCCAG---GCAGUGGGGCCAGGAUUUGGUAGCUGGUGCUGAGAGAAAAC---CC-UUG--AUUGUAUU-CUUGCCCUGGG---AUUAUACCAGUGG +ACA5 UGCAGCCGUGUCAAAUUCAGUACCUGUCCUAUGCAUGGUAGGCACUGGC--CCAGAAGGCUGCCACAGAAACAC--UGUGACUCAUGG------GCCCUGUUCCUGUGUCCCAGGCUCAG +ACA59 GCCCAGGGUAUGUUCACGGGGCGAUGCUGCCCUCCCAGC-UGGCCCAUG--GGUGACCCUGGGAACAUUAACUGCCUCACAACGUUUGUGCCUCAGUUACCCGUAGAUGUAGUGAGGGUA -ACA7 AUUGGA--AGA---CACU-CUGCG-----ACA -ACA30 CAACUG--UCA---CUCA-AUGGG-----ACA -ACA5 GGAUAA--AUU---UGGU-UACAG-----ACA +ACA7 AUUGGAA-----GACACUC-UGC-----GACA +ACA30 CAACUGU-----CACUCAA-UGG-----GACA +ACA5 GGAUAAA-----UUUGGUU-ACA-----GACA ACA59 ACAAUACUUACUCUCGUUGGUGAUAAGGAACA
--- a/test-data/tRNA_2.aln Sat Jan 21 17:38:37 2017 -0500 +++ b/test-data/tRNA_2.aln Mon Dec 05 14:29:00 2022 +0000 @@ -1,7 +1,4 @@ -CLUSTAL W --- LocARNA 1.9.0 --- Score: 2875 +CLUSTAL W --- LocARNA 1.9.2.3 --- Score: 2587 -D10744 GGAAAAUU-GAUCAUCGGCAAGAUAAGUUAUUUACUAAAUAAUAGGAUUUAAUAACCUGG -AF008220 GGAGGAUUAGCUCAGCUGGGAGAGCAUCUGCCUUACAAGCAGAGGG-----------UCG - -D10744 UGAGUUCGAAUCUCACAUUUUCCG -AF008220 GCGGUUCGAGCCCGUCAUCCUCCA +D10744 GGAAAAUUGAUCAUCGGCAAGAUAAGUUAUUUACUAAAUAAUAGGAUUUAAUAACCUGGUGAGUUCGAAUCUCACAUUUUCCG +AF008220 GGAGGAUUAGC--UCAGCUGGGAGAGCAUCUGCCUUA-------CAAGCAGAGGGUCGGCG-GUUCGAGCCCGUCAUCCUCCA
--- a/test-data/tRNA_2.epms Sat Jan 21 17:38:37 2017 -0500 +++ b/test-data/tRNA_2.epms Mon Dec 05 14:29:00 2022 +0000 @@ -1,4 +1,4 @@ epm_id score structure positions -0 3583 ((())) 1:4 2:5 3:6 80:67 81:68 82:69 +0 3587 ((())) 1:4 2:5 3:6 80:67 81:68 82:69 1 3593 ((())) 1:1 2:2 3:3 80:70 81:71 82:72