Mercurial > repos > rnateam > locarna_pairwise_p
changeset 2:535f2e41d4f7 draft default tip
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:23 +0000 |
| parents | 53583f715a8c |
| children | |
| files | locarna.tar.bz2 locarna_pairwise_p.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, 97 insertions(+), 153 deletions(-) [+] |
line wrap: on
line diff
--- a/locarna_pairwise_p.xml Sat Jan 21 17:38:19 2017 -0500 +++ b/locarna_pairwise_p.xml Mon Dec 05 14:29:23 2022 +0000 @@ -1,4 +1,4 @@ -<tool id="locarna_pairwise_p" name="LocARNA Pairwise Probability Aligner" version="@VERSION@.0"> +<tool id="locarna_pairwise_p" name="LocARNA Pairwise Probability Aligner" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@"> <description> Partition Function and Probabilities from Pairwise Simultaneous Alignment and Folding of RNAs (LocARNA-P) </description> @@ -15,14 +15,12 @@ <command><![CDATA[ locarna_p - + '$inputA' '$inputB' - - --width 60 - + @SCORING_ARGS@ - + @FOLDING_ARGS@ @HEURISTIC_ARGS@ @@ -31,12 +29,7 @@ #end if ## -------------------- other parameters - - #if $Constraints.maxBPspan != -1 - --maxBPspan $Constraints.maxBPspan - #end if - - $Constraints.ignore_constraints + @CONSTRAINT_ARGS_WOLP@ ## -------------------- probability parameters $Probabilities.include_am_in_bm @@ -69,11 +62,11 @@ help="Second sequence in fasta, clustal, dp_ps, or PP2.0 format" /> - <param name="outputs" type="select" display="checkboxes" multiple="True" + <param name="outputs" type="select" display="checkboxes" multiple="true" label="Output options"> - <option value="base_match_probs" selected="True"> + <option value="base_match_probs" selected="true"> Output base match probabilities</option> - <option value="arc_match_probs" selected="True"> + <option value="arc_match_probs" selected="true"> Output arc match probabilities</option> </param> @@ -85,16 +78,16 @@ </param> <section name="Probabilities" title="Probability parameters"> - <param name="include_am_in_bm" argument="include-am-in-bm" + <param argument="--include-am-in-bm" type="boolean" truevalue="--include-am-in-bm" falsevalue="" label="Include arc match cases in computation of base match probabilities" /> - <param name="min_bm_prob" argument="min-bm-prob" + <param argument="--min-bm-prob" type="float" value="0.0005" min="0.0" max="1.0" label="Minimal base match probability" /> - <param name="min_am_prob" argument="min-am-prob" + <param argument="--min-am-prob" type="float" value="0.0005" min="0.0" max="1.0" label="Minimal arc match probability" @@ -107,27 +100,19 @@ <section name="Folding" title="RNA folding parameters"> <expand macro="common_folding_parameters" /> - <expand macro="alifold_consensus_parameter" /> - <param name="consensus_structure" type="select" label="Consensus structure type" - help="Type of consensus structures written to screen and stockholm output" - > - <!-- <option value="mea">mea</option>--> - <option value="none">none</option> - <option value="alifold">alifold</option> - </param> </section> <section name="Heuristics" title="Heuristic parameters"> <expand macro="common_heuristic_parameters" /> <expand macro="max_diff_parameters_aln" /> - <param name="pf_scale" argument="pf-scale" + <param argument="--pf-scale" type="float" value="1.0" label="Scaling of the partition function." help="Use (only) if standard scaling overflows." /> </section> <section name="Constraints" title="Constraint parameters"> - <expand macro="common_constraint_parameters" /> + <expand macro="common_constraint_parameters_woLP" /> </section> </inputs> @@ -143,7 +128,7 @@ </outputs> <tests> - <test> + <test expect_num_outputs="2"> <param name="inputA" value="tRNA_2-1.fa" /> <param name="inputB" value="tRNA_2-2.fa" /> <param name="outputs" value="arc_match_probs,base_match_probs"/>
--- a/macros.xml Sat Jan 21 17:38:19 2017 -0500 +++ b/macros.xml Mon Dec 05 14:29:23 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:19 2017 -0500 +++ b/test-data/archaea-default.stdout Mon Dec 05 14:29:23 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:19 2017 -0500 +++ b/test-data/archaea_relplot.scr Mon Dec 05 14:29:23 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:19 2017 -0500 +++ b/test-data/haca.snoRNA-c.aln Mon Dec 05 14:29:23 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:19 2017 -0500 +++ b/test-data/haca.snoRNA-default.stdout Mon Dec 05 14:29:23 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:19 2017 -0500 +++ b/test-data/haca.snoRNA.aln Mon Dec 05 14:29:23 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:19 2017 -0500 +++ b/test-data/tRNA_2.aln Mon Dec 05 14:29:23 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:19 2017 -0500 +++ b/test-data/tRNA_2.epms Mon Dec 05 14:29:23 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