Repository 'hyphy_busted'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/hyphy_busted

Changeset 26:70fab57f3c05 (2021-04-27)
Previous changeset 25:b6509de48a62 (2021-04-20) Next changeset 27:364eda3cc633 (2021-06-09)
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/hyphy/ commit 3908306ffa6f6cbc7b83303b2d4c581406d2fad9"
modified:
hyphy_busted.xml
macros.xml
test-data/busted-out1.json
test-data/gard-out1.json
test-data/gard-out1.nex
test-data/relax-out1.json
added:
test-data/busted-in2.nex
test-data/busted-out2-1.json
test-data/busted-out2-2.json
test-data/relax-in2.nex
test-data/relax-out2-1.json
test-data/relax-out2-2.json
b
diff -r b6509de48a62 -r 70fab57f3c05 hyphy_busted.xml
--- a/hyphy_busted.xml Tue Apr 20 19:54:45 2021 +0000
+++ b/hyphy_busted.xml Tue Apr 27 18:06:02 2021 +0000
[
@@ -1,33 +1,30 @@
 <?xml version="1.0"?>
-<tool id="hyphy_busted" name="HyPhy-BUSTED" version="@VERSION@+galaxy0" profile="19.09">
+<tool id="hyphy_busted" name="HyPhy-BUSTED" version="@VERSION@+galaxy1" profile="19.09">
     <description>Branch-site Unrestricted Statistical Test for Episodic Diversification</description>
     <macros>
         <import>macros.xml</import>
-        <token name="@OPERATION@">BUSTED</token>
-        <token name="@operation@">busted</token>
     </macros>
     <expand macro="requirements"/>
     <command detect_errors="exit_code"><![CDATA[
         @SYMLINK_FILES@
         hyphy busted
             --alignment ./$input_file
-            #if $input_nhx:
-                --tree ./input.nhx
-            #end if
+            @INPUT_TREE@
             --code $gencodeid
-            --branches '$branches'
+            @branch_options@
             --output '$busted_output'
-            > '$busted_log' ;
-        @CATCH_ERROR@
+            @advanced_options@
+        @ERRORS@
         ]]></command>
     <inputs>
         <expand macro="inputs"/>
         <expand macro="gencode"/>
         <expand macro="branches"/>
+        <expand macro="srv"/>
     </inputs>
     <outputs>
-        <data name="busted_log" format="txt"/>
         <data name="busted_output" format="hyphy_results.json" />
+        <expand macro="alternative_model_output" />
     </outputs>
     <tests>
         <test>
@@ -35,6 +32,15 @@
             <param name="input_nhx" ftype="nhx" value="absrel-in1.nhx"/>
             <output name="busted_output" file="busted-out1.json" compare="sim_size"/>
         </test>
+        <test>
+            <param name="input_file" ftype="nex" value="busted-in2.nex"/>
+            <param name="branches" value="specify" />
+            <param name="branch_label" value="B_US_90_WEAU160" />
+            <param name="srv_options" value="specify" />
+            <param name="save_alternative_model" value="true" />
+            <output name="busted_output" file="busted-out2-1.json" compare="sim_size"/>
+            <output name="alternative_model" file="busted-out2-2.json" compare="sim_size"/>
+        </test>
     </tests>
     <help><![CDATA[
 
b
diff -r b6509de48a62 -r 70fab57f3c05 macros.xml
--- a/macros.xml Tue Apr 20 19:54:45 2021 +0000
+++ b/macros.xml Tue Apr 27 18:06:02 2021 +0000
[
@@ -93,13 +93,45 @@
         </param>
     </xml>
 
+    <xml name="srv">
+        <conditional name="advanced" label="Advanced options">
+            <param name="srv_options" type="select" label="Set advanced parameters" help="If these parameters are not set, RELAX will calculate appropriate values at runtime">
+                <option value="defaults">Calculate</option>
+                <option value="specify">Specify values</option>
+            </param>
+            <when value="defaults" />
+            <when value="specify">
+                <param argument="--grid-size" type="integer" value="250" label="Points in the initial distributional guess for likelihood fitting" />
+                <param argument="--starting-points" type="integer" value="1" label="Initial random guesses to seed rate values optimization" />
+                <param argument="--syn-rates" type="integer" min="1" max="10" value="3" label="Alpha rate classes to include in the model" />
+                <param argument="--rates" type="integer" min="2" max="10" value="3" label="Omega rate classes to include in the model" />
+                <param argument="--srv" type="boolean" truevalue="Yes" falsevalue="No" label="Include synonymous rate variation" />
+                <param name="save_alternative_model" type="boolean" label="Save alternative model fit" />
+            </when>
+        </conditional>
+    </xml>
+
     <xml name="branches">
-        <param name="branches" type="select" label="Set of branches to test">
-            <option value="All">All branches</option>
-            <option value="Internal">Internal branches</option>
-            <option value="Leaves">Leaf branches</option>
-            <option value="'Unlabeled-branches'">Unlabeled branches</option>
-        </param>
+        <conditional name="branch_cond">
+            <param name="branch_sel" type="select" label="Set of branches to test">
+                <option value="All">All branches</option>
+                <option value="Internal">Internal branches</option>
+                <option value="Leaves">Leaf branches</option>
+                <option value="'Unlabeled-branches'">Unlabeled branches</option>
+                <option value="specify">Enter a branch label</option>
+            </param>
+            <when value="specify">
+                <param name="branch_label" type="text" value="Test" optional="false" label="Branch label">
+                    <sanitizer invalid_char="">
+                        <valid initial="default" />
+                    </sanitizer>
+                </param>
+            </when>
+            <when value="All" />
+            <when value="Internal" />
+            <when value="Leaves" />
+            <when value="'Unlabeled-branches'" />
+        </conditional>
     </xml>
     <xml name="citations">
         <citations>
@@ -114,17 +146,40 @@
             <yield/>
         </requirements>
     </xml>
+    <xml name="alternative_model_output">
+        <data name="alternative_model" format="hyphy_results.json" from_work_dir="alternative_model.json">
+            <filter>advanced['save_alternative_model']</filter>
+        </data>
+    </xml>
+    <token name="@INPUT_TREE@"><![CDATA[
+            #if $input_nhx:
+                --tree input.nhx
+            #end if
+        ]]>
+    </token>
+    <token name="@branch_options@"><![CDATA[
+        #if $branch_cond.branch_sel == 'specify':
+            --branches '$branch_cond.branch_label'
+        #else:
+            --branches $branch_cond.branch_sel
+        #end if
+        ]]>
+    </token>
+    <token name="@advanced_options@"><![CDATA[
+        #if $advanced.srv_options == 'specify':
+            --grid-size $advanced.grid_size
+            --starting-points $advanced.starting_points
+            --syn-rates $advanced.syn_rates
+            --rates $advanced.rates
+            --srv $advanced.srv
+            #if $advanced.save_alternative_model:
+                --save-fit alternative_model.json
+            #end if
+        #end if
+    ]]></token>
     <token name="@HYPHYMPI@">\${GALAXY_MPIRUN:-mpirun -mca orte_tmpdir_base "\${TMPDIR:-.}" -np \${GALAXY_SLOTS:-1}} HYPHYMPI</token>
-    <token name="@CATCH_ERROR@"><![CDATA[
-        EC=\$? ;
-        if [ \$EC -ne 0 ] ; then
-            if [ -f errors.log.mpinode0 ] ; then
-                cat errors.log.mpinode0 >&2 ;
-            else 
-                cat errors.log >&2 ;
-            fi ;
-        fi ;
-        exit \$EC
+    <token name="@ERRORS@"><![CDATA[
+        ; EC=\$? ; [[ -f errors.* ]] && cat errors.* >&2 ; exit \$EC
     ]]></token>
     <token name="@HYPHY_ENVIRONMENT@"><![CDATA[
         export HYPHY=`which hyphy` &&
@@ -134,11 +189,10 @@
         @HYPHY_ENVIRONMENT@ hyphy LIBPATH=\$HYPHY_LIB
     ]]></token>
     <token name="@SYMLINK_FILES@"><![CDATA[
-        ln -s '$input_file' input.$input_file.extension &&
-        ln -s '$@operation@_output' input.${input_file.extension}.@OPERATION@.json &&
-        #set $input_file = 'input.%s' % $input_file.extension
-        #if $input_nhx:
+        #if $input_nhx is not None:
             ln -s '$input_nhx' input.nhx &&
         #end if
+        ln -s '$input_file' input.$input_file.extension &&
+        #set $input_file = 'input.%s' % $input_file.extension
     ]]></token>
 </macros>
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/busted-in2.nex
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/busted-in2.nex Tue Apr 27 18:06:02 2021 +0000
b
b"@@ -0,0 +1,30 @@\n+#NEXUS\n+\n+BEGIN TAXA;\n+\tDIMENSIONS NTAX = 8;\n+\tTAXLABELS\n+\t\t'B_FR_83_HXB2' 'B_US_83_RF' 'B_US_86_JRFL' 'B_US_90_WEAU160' 'D_CD_83_ELI' 'D_CD_83_NDK' 'D_CD_84_84ZR085' 'D_UG_94_94UG114' ;\n+END;\n+\n+BEGIN CHARACTERS;\n+\tDIMENSIONS NCHAR = 1320;\n+\tFORMAT\n+\t\tDATATYPE = DNA\n+\t\tGAP=-\n+\t\tMISSING=?\n+\t;\n+\n+MATRIX\n+\t'B_FR_83_HXB2'     CCCATTAGCCCTATTGAGACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTTAAACAATGGCCATTGACAGAAGAAAAAATAAAAGCATTAGTAGAAATTTGTACAGAGATGGAAAAGGAAGGGAAAATTTCAAAAATTGGGCCTGAAAATCCATACAATACTCCAGTATTTGCCATAAAGAAAAAAGACAGTACTAAATGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGACTTCTGGGAAGTTCAATTAGGAATACCACATCCCGCAGGGTTAAAAAAGAAAAAATCAGTAACAGTACTGGATGTGGGTGATGCATATTTTTCAGTTCCCTTAGATGAAGACTTCAGGAAGTATACTGCATTTACCATACCTAGTATAAACAATGAGACACCAGGGATTAGATATCAGTACAATGTGCTTCCACAGGGATGGAAAGGATCACCAGCAATATTCCAAAGTAGCATGACAAAAATCTTAGAGCCTTTTAGAAAACAAAATCCAGACATAGTTATCTATCAATACATGGATGATTTGTATGTAGGATCTGACTTAGAAATAGGGCAGCATAGAACAAAAATAGAGGAGCTGAGACAACATCTGTTGAGGTGGGGACTTACCACACCAGACAAAAAACATCAGAAAGAACCTCCATTCCTTTGGATGGGTTATGAACTCCATCCTGATAAATGGACAGTACAGCCTATAGTGCTGCCAGAAAAAGACAGCTGGACTGTCAATGACATACAGAAGTTAGTGGGGAAATTGAATTGGGCAAGTCAGATTTACCCAGGGATTAAAGTAAGGCAATTATGTAAACTCCTTAGAGGAACCAAAGCACTAACAGAAGTAATACCACTAACAGAAGAAGCAGAGCTAGAACTGGCAGAAAACAGAGAGATTCTAAAAGAACCAGTACATGGAGTGTATTATGACCCATCAAAAGACTTAATAGCAGAAATACAGAAGCAGGGGCAAGGCCAATGGACATATCAAATTTATCAAGAGCCATTTAAAAATCTGAAAACAGGAAAATATGCAAGAATGAGGGGTGCCCACACTAATGATGTAAAACAATTAACAGAGGCAGTGCAAAAAATAACCACAGAAAGCATAGTAATATGGGGAAAGACTCCTAAATTTAAACTGCCCATACAAAAGGAAACATGGGAAACATGGTGGACAGAGTATTGGCAAGCCACCTGGATTCCTGAGTGGGAGTTTGTTAATACCCCTCCCTTAGTGAAATTATGGTACCAGTTAGAGAAAGAACCCATAGTAGGAGCAGAAACCTTC\n+\t'B_US_83_RF'       CCCATTAGTCCTATTGAAACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTTAAACAATGGCCATTGACAGAGGAAAAAATAAAAGCATTGGTAGAAATTTGTACAGAAATGGAAAAGGAAGGAAAAATTTCCAAAATTGGGCCTGAAAATCCATACAATACTCCAGTATTTGCCATAAAGAAAAAAGACAGTACTAAATGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGACTTCTGGGAAGTTCAGTTAGGAATACCACATCCTGCAGGGTTAAAAAAGAAGAAATCAGTAACAGTATTGGATGTGGGTGATGCATATTTTTCAGTTCCCTTAGATAAAGAGTTCAGGAAGTATACTGCATTTACCATACCTAGTATAAACAATGAAACACCACGGATTAGATATCAGTACAATGTGCTTCCACAAGGGTGGAAAGGATCACCAGCAATATTCCAAAGTAGTATGACAAAAATCTTAGAGCCTTTTAAAAAACAAAATCCAGAAATAGTTATCTATCAATACATGGATGATTTGTATGTAGGATCTGATTTAGAAATAGGGCAGCATAGAATAAAAATAGAGGAACTGAGAGAACATCTGTTAAAGTGGGGGTTTACCACACCGGACAAGAAACATCAGAAAGAACCTCCATTTCTTTGGATGGGTTATGAACTCCATCCTGATAAATGGACAGTACAGCCTATAGTGCTGCCAGAAAAAGACAGCTGGACTGTCAATGACATACAGAAGTTAGTGGGAAAATTGAATTGGGCAAGTCAGATTTATGCAGGGATTAAAGTAAAGCAATTATGTAAACTCCTTAGGGGAACCAAAGCACTAACAGAAGTAGTACAACTAACAAAAGAAGCAGAGCTAGAACTGGCAGAAAATAGGGAGATTCTAAAAGAACCAGTACATGGAGTGTATTATGACCCATCAAAAGACTTAATAGCAGAAATACAGAAGCAGGGGCAAGGCCAATGGACATACCAAATTTATCAAGAGCCATTTAAAAACCTGAAAACAGGAAAGTATGCAAGAATGAGGGGTGCCCACACTAATGATGTAAAACAATTAACAGAGGCAGTACAAAAAGTAGCCACAGAAAGCATAGTAATATGGGGAAAGACTCCTAAATTTAAACTACCCATACAAAAAGAAACATGGGAGGCATGGTGGACAGAGTATTGGCAAGCCACCTGGATTCCTGAGTGGGAGTTTGTCAATACCCCTCCCTTAGTAAAATTGTGGTACCAGTTAGAAAAAGAACCCATAATAGGAGCAGAAACTTTC\n+\t'B_US_86_JRFL'     CCCATTAGTCCTATTGAAACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTCAAACAATGGCCATTGACAGAAGAAAAAATAAAAGCATTAGTAGAAATTTGTACAGAAATGGAAAAGGAAGGGAAAATTTCAAAAATTGGGCCTGAAAATCCATACAATACTCCAGTATTTGCCATAAAGAAAAAGGACAGTACTAAATGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAAAACTCAAGACTTCTGGGAAGTTCAATTAGGAATACCACATCCCGCAGGGTTAAAAAAGAGAAAATCAGTAACAGTACTGGATGTGGGTGATGCATATTTTTCAGTTCCCTTAGATAAAGACTTCAGGAAATATACTGCATTTACCATACCTAGTATAAACAATGAGACACCAGGGATTAGGTATCAGTACAATGTGCTTCCGCAGGGATGGAAAGGATCACCAGCAATATTCCAAAGTAGCATGACAAAAATCTTAGAGCCTTTTAGAAAACAAAATCCAGACATAATTATCTATCAATACATGGATGATTTGTATGTAGGATCTGACTTAGAGATAGGGCAGCATAGAGCAAAAATAGAGGAATTGAGACAACATCTGTTGAGGTGGGGGTTTACCACACCAGACAAAAAACATCAGAAAGAACCTCCATTCCTTTGGATGGGTTATGAACTCCATCCTGACAAATGGACAGTACAGCCTATAGTGCTGCCAGAAAAAGACAGCTGGACTGTCAATGACATACAGAAGTTAGTGGGAAAATTAAATTGGGCAAGTCAGATTTACGCAGGGATTAAAGTAAAGCAATTATGTAAACTCCTTAGGGGAACCAAAGCACTAACAGAAGTAATACCACTAACAGAAGAAGCAGAGCTAGAACTGGCAGAAAACAGGGAGATTCTAAAAGAGCCAGTACATGGAGTGTATTATGACCCATCAAAAGAC"..b"ACTCCAATATTTGCCATAAAGAAAAAAGACAGTACCAAGTGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGATTTCTGGGAGGTTCAATTAGGAATACCGCATCCTGCAGGGCTGAAAAAGAAAAAATCAGTAACAGTACTGGATGTGGGTGATGCATATTTCTCAGTTCCCTTAGATGAAGATTTTAGGAAATATACCGCATTTACCATACCTAGTATAAACAATGAGACACCAGGGATTAGATATCAGTACAATGTGCTCCCACAGGGATGGAAAGGATCACCGGCAATATTCCAAAGTAGCATGACAAAAATCTTAGAGCCCTTTAGAAAACAAAATCCAGAAATAGTTATCTATCAATACATGGATGATTTGTATGTAGGATCTGACTTAGAAATAGGGCAGCATAGAACAAAAATAGAGGAATTAAGAGAACATCTATTGAGGTGGGGATTTACCACACCAGATAAAAAACATCAGAAAGAACCTCCATTTCTTTGGATGGGTTATGAACTCCATCCTGATAAATGGACAGTACAGCCTATAAACCTGCCAGAAAAAGAAAGCTGGACTGTCAATGATATACAGAAGTTAGTGGGGAAATTAAACTGGGCAAGCCAGATTTATGCAGGAATTAAAGTAAAGCAATTATGTAAACTCCTTAGGGGAACCAAAGCACTAACAGAAGTAGTACCACTAACAGAAGAAGCAGAATTAGAACTGGCAGAAAACAGGGAAATTCTAAAAGAACCAGTACATGGAGTGTATTATGACCCATCAAAAGACTTAATAGCAGAACTACAGAAACAAGGGGACGGCCAATGGACATACCAAATTTATCAAGAACCATTTAAAAATCTAAAAACAGGAAAGTATGCAAGAACGAGGGGTGCCCACACTAATGATGTAAAACAATTAACAGAGGCAGTGCAAAAAATAGCCACAGAAAGCATAGTGATATGGGGAAAGACTCCTAAATTTAAACTACCCATACAAAAGGAAACATGGGAAACATGGTGGATAGAGTATTGGCAAGCCACCTGGATTCCTGAGTGGGAATTTGTCAATACCCCTCCTTTAGTAAAATTATGGTACCAGTTAGAGAAGGAACCCATAATAGGAGCAGAAACTTTC\n+\t'D_CD_84_84ZR085'  CCAATTAGTCCTATTGAAACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTTAAACAATGGCCGTTGACAGAAGAAAAAATAAAAGCATTAACAGAAATTTGTACAGATATGGAAAAGGAAGGAAAAATTTCAAGAATTGGGCCTGAAAATCCATACAATACTCCAATATTTGCCATAAAGAAAAAAGACAGTACTAAGTGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGACTTCTGGGAAGTTCAATTAGGGATACCACATCCTGCAGGATTAAAGAAGAAAAAGTCAATAACAGTACTGGATGTGGGCGATGCATATTTTTCAATTCCCTTATGTGAAGACTTTAGGAAGTACACTGCATTTACCATACCTAGTATAAACAATGAGACACCAGGGATTAGATATCAGTACAATGTACTTCCACAGGGATGGAAAGGATCACCAGCAATATTCCAAAGTAGCATGATAAAAATCTTAGAGCCCTTTAGAAAACAAAATCCAGAAGTAGTTATCTATCAATACATGGATGATTTGTATGTAGGATCTGATTTAGAAATAGGACAGCATAGAGCAAAAATAGAGAAATTAAGAGAACATCTGTTGAGGTGGGGGCTTACCACACCAGACAAAAAACATCAGAAAGAACCTCCATTTCTTTGGATGGGTTATGAACTCCATCCTGATAAGTGGACAGTACAGTCTATAACACTGCCAGAGAAAGAAAGCTGGACTGTCAATGATATACAGAAGTTAGTGGGAAAATTAAATTGGGCAAGCCAGATTTATCCAGGAATTAAAGTAAAGCAATTATGTAAACTCCTTAGGGGAACCAAGGCACTAACAGAGGTAATACCACTAACAGAAGAAGCAGAATTAGAACTGGCAGAAAACAGGGAGATTCTAAAGGAACCAATGCATGGAGTGTATTATGACCCATCAAAAGACTTAATAGCAGAATTACAGAAACAAGGGCAAGGTCAATGGACATATCAAATTTATCAAGAACCATTTAAAAATCTGAAAACAGGAAAGTATGCAAGAATGAGGGGTGCCCACACTAATGATGTAAAACAGTTAACAGAGGCAGTGCAAAAAATAGCCATAGAAAGCATAGTGATATGGGGAAAGACTCCTAAATTTAGACTACCCATACAAAAGGAAACATGGGAAACATGGTGGATAGACTATTGGCAAGCCACCTGGATTCCTGAGTGGGAATTTGTCAATACCCCTCCTTTAGTAAAATTATGGTACCAGTTAGAGAAGGAACCCATAATAGGAGCAGAAACTTTC\n+\t'D_UG_94_94UG114'  CCAATTAGTCCTATTGAAACTGTACCAGTAAAATTAAAGCCAGGGATGGATGGCCCAAAAGTTAAACAATGGCCGTTGACAGAAGAAAAAATAAAAGCACTAATAGAAATTTGTTCAGAACTAGAAAAGGAAGGAAAAATTTCAAAAATTGGGCCTGAAAACCCATACAATACTCCAATATTTGCCATAAAGAAAAAAGACAGTACTAAGTGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGACTTTTGGGAAGTTCAACTAGGAATACCACATCCTGCAGGGCTAAAAAAGAAAAAATCAGTAACAGTACTGGATGTGGGTGACGCATATTTTTCAGTTCCCTTACATGAAGACTTTAGAAAATATACCGCATTCACCATACCTAGTACAAACAATGAGACACCAGGAATTAGATATCAGTACAATGTGCTTCCACAAGGATGGAAAGGATCACCAGCAATATTCCAAAGTAGCATGACAAAAATCTTAGAACCTTTTAGAAAACAAAATCCAGAAATGATTATCTATCAATACATGGATGATTTGTATGTAGGATCTGACTTAGAAATAGGGCAGCATAGAATAAAAATAGAGGAATTAAGGGGACACCTCTTGAAGTGGGGATTTACCACACCAGACAAAAAGTATCAGAAAGAACCCCCATTTCTTTGGATGGGTTATGAACTCCATCCTGATAAGTGGACAGTACAGCCTATACATCTGCCAGAAAAGGAAAGCTGGACTGTCAATGATATACAGAAGTTAGTGGGAAAATTAAATTGGGCAAGCCAGATTTATCCAGGAATTAAAGTAAGACAATTATGCAAATGCCTTAGGGGAGCCAAAGCACTGACAGAAGTAATACCACTGACAGCAGAAGCAGAATTAGAACTGGCAGAAAACAGGGAAATACTAAAAGAACCAGTACATGGAGCGTATTATGACCCATCAAAAGACTTAATAGCAGAAATACAGAAACAAGGGCAAGATCAATGGACATATCAAATATATCAAGAACAATATAAAAATCTGAAAACAGGAAAGTATGCGAAAATGAGGGGTACCCACACTAATGATGTAAAACAATTAACAGAGGCAGTGCAGAAAATAGCCCAAGAATGTATAGTAATATGGGGAAAGACTCCTAAATTTAGACTACCCATACAAAAGGAAACATGGGAAACATGGTGGACAGAGTATTGGCAGGCCACCTGGATTCCTGAGTGGGAGTATGTCAACACCCCTCCTTTAGTTAAATTATGGTATCAGTTAGAGAAGGAACCCATAGTAGGAGCAGAAACTTTC;\n+END;\n+\n+BEGIN TREES;\n+\tTREE tree = ((((D_CD_83_ELI,D_CD_83_NDK),D_UG_94_94UG114),D_CD_84_84ZR085),B_US_83_RF,((B_FR_83_HXB2,B_US_86_JRFL),B_US_90_WEAU160));\n+END;\n"
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/busted-out1.json
--- a/test-data/busted-out1.json Tue Apr 20 19:54:45 2021 +0000
+++ b/test-data/busted-out1.json Tue Apr 27 18:06:02 2021 +0000
[
b'@@ -1,73 +1,73 @@\n {\n  "Evidence Ratios":{\n    "constrained":    [\n-[1.766356514902524, 0.9484374805763239, 0.9571091930814299, 1.845676838710576, 1.513819401189324, 2.11256226458989, 7.366103227625151, 3.937756516229049, 2.586067836907899, 0.722379841928909, 0.3151239046214458, 1.502917468627881, 0.6499879033472316, 1.047489715978757, 1.060206108705189, 0.5038296392390414, 0.8532658611088545, 1.114374823999701, 0.5611409841208096, 0.6702870396165613, 0.4604115167079821, 0.3988433224622866, 2.47019807349414, 0.9858342091732886, 1.106653406611749, 1.067348070135177, 1.77412165088655, 1.321424873505688, 1.90457149631063, 0.7888758683589224, 0.7229203540914775, 0.8813401886687667, 0.6954297746421804, 4.613654710605887, 0.4975309466317643, 1.84565076502724, 0.6214894235132998, 0.9480915704528677, 0.9669462933548731, 1.488081789562878, 1.368539685217942, 1.726522735566884, 1.222915324989829, 0.7085094082639855, 0.3168181005885327, 2.290014391709923, 0.6343005734674315, 0.8841633809961478, 1.508078414734822, 2.666352901083321, 1.732923983135172, 3.102011384229053, 1.585109653073893, 3.345181824497555, 7.193115253213268, 3.850321954428519, 3.233128500312913, 1.977820966308441, 2.047653340750355, 0.6643430465826362, 3.624272911205855, 1.67814994913239, 1.752817561416471, 3.277248182160172, 0.3406457937754578, 0.4841828665456315, 1.399775481724447, 0.5510187898414104, 0.6276769404480749, 0.4535822272441153, 0.4226923979745105, 2.950051361271033, 0.9999999999999982, 1.027733519332009, 2.963278030618241, 2.8935730521221, 1.881916301985038, 0.4490192221083247, 3.102588216900087, 0.7028469232610687, 2.792137897266834, 0.3889955246330017, 0.9360790804326219, 0.4572233386404743, 2.908919410553763, 0.8995968442944373, 0.3261344314757346, 1.080257865601265, 1.399763176607792, 1.534646955858424, 0.4277164444783946, 2.361993404905034, 2.175584822495568, 3.122529901767528, 0.4453847841440882, 1.711825080925545, 0.8458710509917922, 2.741896592382772, 0.576439632416141, 1.142294357680513, 0.8086330800808955, 2.943414796294152, 0.833291758377592, 1.615921141634484, 0.7252129604805458, 2.833655116179036, 0.9138779988368123, 0.6174720822006171, 1.834532902825704, 0.2316812685220093, 2.479235579982967, 0.4643862206171616, 6.530720225707355, 0.4975309466317643, 1.158125132567293, 0.2987944638927862, 3.78696693228073, 0.7343025200073082, 1.93339742618843, 0.3642707474760721, 1.028034824972818, 1.511057897004292, 0.2969890376651506, 0.987605711796378, 0.3709886418939366, 4.185741574214752, 3.002372601666337, 0.424383367760914, 0.4658081554322532, 0.4024556436060775, 2.540164510486854, 0.5199349857684811, 0.4626485320487528, 0.9349816613814522, 0.5525607476128181, 0.5855612460554309, 1.85872747574379, 1.064964125964651, 3.177687923317396, 2.593304657270028, 5.234187858164534, 2.966427299030309, 1.920561682721974, 2.245869881969924, 3.811870935099101, 2.099654636695842, 3.369504334172069, 0.6379165155453101, 3.119867019567397, 0.9211859677702382, 2.195241379859611, 1.767535349740496, 1.225677095999665, 1.347676192394215, 1.671173168407676, 1.179430472478473, 1.873629958128968, 1.972321068578892, 1.863249013159319, 1.501711930393218, 2.85458165498371, 0.5332817404709185, 0.7514817421990073, 0.2887801591578705, 1.93925520496415, 0.2163710715496819, 2.457406602553581, 0.4895227461773333, 2.731573653164397, 0.8923840424261161, 0.6288736610053396, 3.562511382204229, 0.8043427141817341, 2.97831687935463, 3.021636055180988, 1.919432168129229, 1.231674183402057, 3.382911261347008, 1.455440822286478, 1.714894279904311, 0.4763853457361669, 0.907985271671091, 0.9191891327337078, 0.8704988461364657, 0.4660185191281315, 0.6176233729242167, 0.9987767701341157] \n+[1.751697313438667, 0.9490151622254345, 0.9575938099872893, 1.827482077500528, 1.505094606485015, 2.088346790071941, 7.070051504373732, 3.867967952992799, 2.545750919582881, 0.7247660693488289, 0.3196737937441124, 1.489449650902744, 0.6530151351740718, 1.04337285373284, 1.058333334239147, 0.5088801718500566, 0.8'..b'0.01189287130232129] \n       ],\n-     "Log Likelihood":-3467.004808147542,\n+     "Log Likelihood":-3467.012706765863,\n      "Rate Distributions":{\n        "non-synonymous/synonymous rate ratio for *test*":        [\n-[0.9950576192409499, 1] \n+[0.9959842858928788, 1] \n         ]\n       },\n      "display order":1,\n      "estimated parameters":31\n     },\n    "Nucleotide GTR":{\n-     "AIC-c":7112.85745601076,\n+     "AIC-c":7112.857310169033,\n      "Equilibrium frequencies":      [\n [0.3592490842490842],\n       [0.181959706959707],\n       [0.240018315018315],\n       [0.2187728937728937] \n       ],\n-     "Log Likelihood":-3532.321297387654,\n+     "Log Likelihood":-3532.32122446679,\n      "Rate Distributions":{\n-       "Substitution rate from nucleotide A to nucleotide C":0.5497129763145497,\n+       "Substitution rate from nucleotide A to nucleotide C":0.5501644800943087,\n        "Substitution rate from nucleotide A to nucleotide G":1,\n-       "Substitution rate from nucleotide A to nucleotide T":0.264666902964203,\n-       "Substitution rate from nucleotide C to nucleotide G":0.4916303643305098,\n-       "Substitution rate from nucleotide C to nucleotide T":1.027924764411573,\n-       "Substitution rate from nucleotide G to nucleotide T":0.304615404570875\n+       "Substitution rate from nucleotide A to nucleotide T":0.2647746107611847,\n+       "Substitution rate from nucleotide C to nucleotide G":0.4917423221832032,\n+       "Substitution rate from nucleotide C to nucleotide T":1.02837595903795,\n+       "Substitution rate from nucleotide G to nucleotide T":0.304509117524575\n       },\n      "display order":0,\n      "estimated parameters":24\n     },\n    "Unconstrained model":{\n-     "AIC-c":6910.766610587704,\n-     "Log Likelihood":-3414.486640449675,\n+     "AIC-c":6910.776988890891,\n+     "Log Likelihood":-3414.491829601268,\n      "Rate Distributions":{\n        "Synonymous site-to-site rates":{\n          "0":{\n-           "proportion":0.2645541257391169,\n-           "rate":0.2747941124675511\n+           "proportion":0.2655037162058198,\n+           "rate":0.2749303122566129\n           },\n          "1":{\n-           "proportion":0.6325830460054784,\n-           "rate":0.9734675888351726\n+           "proportion":0.6321791898264187,\n+           "rate":0.9740711294478059\n           },\n          "2":{\n-           "proportion":0.1028628282554047,\n-           "rate":3.028333913894314\n+           "proportion":0.1023170939677615,\n+           "rate":3.041695877803997\n           }\n         },\n        "Test":{\n          "0":{\n            "omega":0,\n-           "proportion":0.04809457608200497\n+           "proportion":0.386646389823596\n           },\n          "1":{\n-           "omega":0,\n-           "proportion":0.3410012743721951\n+           "omega":0.9964438830494696,\n+           "proportion":0\n           },\n          "2":{\n-           "omega":1.89363684081788,\n-           "proportion":0.6109041495458\n+           "omega":1.876622752045829,\n+           "proportion":0.613353610176404\n           }\n         }\n       },\n@@ -341,7 +341,7 @@\n     }\n   },\n  "input":{\n-   "file name":"/tmp/tmp7qHfo2/job_working_directory/000/3/working/./busted_input.fa",\n+   "file name":"/tmp/tmpfqvtss8b/job_working_directory/000/3/working/./input.fasta",\n    "number of sequences":10,\n    "number of sites":187,\n    "partition count":1,\n@@ -350,8 +350,8 @@\n     }\n   },\n  "test results":{\n-   "LRT":10.17362366082671,\n-   "p-value":0.003088842014430382\n+   "LRT":10.17667489929045,\n+   "p-value":0.003084133210483175\n   },\n  "tested":{\n    "0":{\n@@ -376,19 +376,19 @@\n  "timers":{\n    "Constrained BUSTED model fitting":{\n      "order":3,\n-     "timer":62\n+     "timer":70\n     },\n    "Overall":{\n      "order":0,\n-     "timer":228\n+     "timer":451\n     },\n    "Preliminary model fitting":{\n      "order":1,\n-     "timer":3\n+     "timer":10\n     },\n    "Unconstrained BUSTED model fitting":{\n      "order":2,\n-     "timer":161\n+     "timer":370\n     }\n   }\n }\n\\ No newline at end of file\n'
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/busted-out2-1.json
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/busted-out2-1.json Tue Apr 27 18:06:02 2021 +0000
[
b'@@ -0,0 +1,338 @@\n+{\n+ "Evidence Ratios":{\n+   "constrained":    [\n+[0.935557169656361, 0.9478485707814517, 0.9365713272643597, 0.9348516866501475, 0.9478485707814517, 0.9470782825238451, 0.9404255507056813, 0.9445269614625446, 0.9357171171944089, 0.9445269614625446, 0.951038343352286, 0.9642191069363555, 0.9509681609268663, 0.9357171171944089, 0.9365508162861411, 0.9330337048596494, 0.9453769311406178, 0.9354507653189286, 0.9357171171944089, 0.951038343352286, 0.944959794342383, 0.951038343352286, 0.9554450473318812, 0.957741198595587, 0.9357223141306646, 0.9635859922205572, 0.9400812190552175, 0.9471227276423632, 0.9471155585144405, 0.951038343352286, 0.9425085292914046, 0.951038343352286, 0.9382075921692818, 0.9681439336221066, 9.055411306947674, 0.9467899324264936, 0.9478485707814517, 0.9398061397692539, 1.097085569751409, 1.371682580705311, 1.704754646360477, 0.9467899324264936, 0.9509681609268663, 0.9467899324264936, 0.9364415899930733, 0.951038343352286, 0.9478485707814517, 0.9525827524658199, 1.146813959542357, 0.9478485707814517, 0.9362276370626634, 0.9348516866501475, 0.9467899324264936, 0.9486832156560699, 0.9357171171944089, 0.9541686971206805, 0.9484511531422709, 0.9404255507056813, 0.9357171171944089, 1.072863229391129, 0.947891315956519, 0.938167627826613, 0.9425085292914046, 0.9509681609268663, 0.951038343352286, 0.9514408733975561, 0.9395660064639854, 0.9365128450433999, 0.9404014893688138, 0.9514344449206511, 0.957741198595587, 0.9388792918800136, 0.951038343352286, 0.9642191069363555, 0.9445269614625446, 0.9453769311406178, 0.9418928080136506, 0.9388792918800136, 0.9467899324264936, 0.9399880829471382, 0.9484511531422709, 0.9509681609268663, 1.153990916133998, 0.9404255507056813, 0.9554450473318812, 0.9412338215355154, 0.942333955072914, 0.957741198595587, 0.9470815762343416, 0.9449667231570663, 0.9557119925844111, 0.9676532717080474, 0.9365050172103501, 0.9425085292914046, 0.9357012899378949, 0.9407207493748144, 0.9348483936808639, 1.953541836729934, 0.9362670029301522, 0.9644232801771531, 0.9515108647313498, 0.9509681609268663, 1.168382764503998, 0.9515108647313498, 0.949860271872106, 1.047791522454514, 0.9400812190552175, 0.9445269614625446, 0.9566068656126535, 0.9453769311406178, 0.9434809129833114, 0.936045980720387, 0.9455636023716772, 0.9382075921692818, 0.9549292932384192, 0.9479532753114727, 0.949860271872106, 1.071924199060109, 0.9347728559684337, 0.9642191069363555, 61.30480048708166, 1.392589724214399, 1.068829280572215, 0.9452448851782124, 0.9383322213205718, 0.9520356968106081, 0.9548771157697881, 0.9403109573428395, 0.9383425168225155, 0.9481200597562297, 0.9401597754529221, 0.9425085292914046, 1.081655864332756, 0.9365128450433999, 0.958884569235499, 0.9416389648192759, 0.9484511531422709, 0.9474364959906919, 0.9400812190552175, 0.9357171171944089, 1.013189964747481, 0.9478485707814517, 0.939061940389555, 0.9549292932384192, 0.9544426850176303, 0.9541087300514087, 0.9484511531422709, 0.9456450378925056, 0.9399863915898927, 0.9357147779844784, 0.9552428339477506, 0.9365032950653777, 0.957741198595587, 0.951038343352286, 0.9364880871021735, 0.949860271872106, 0.9357012899378949, 0.9382075921692818, 0.9425085292914046, 0.9418928080136506, 0.9554450473318812, 0.9365128450433999, 0.9333687495168144, 0.9330337048596494, 1.030183566875361, 0.951038343352286, 0.9530190340820841, 0.9642191069363555, 0.9471668565289959, 0.9348308160631009, 0.947891315956519, 1.060517303033454, 0.951038343352286, 0.9554450473318812, 0.9484511531422709, 0.9357171171944089, 1.428578774814407, 1.256482256153546, 1.21301128797193, 0.9447052888159178, 0.9549292932384192, 0.9554450473318812, 0.9541087300514087, 0.9330337048596494, 0.9453769311406178, 0.9453769311406178, 0.9635859922205572, 0.9549292932384192, 0.9445269614625446, 0.9364880871021735, 0.9422504512769631, 0.9408248422372665, 0.9642191069363555, 0.9471081975229871, 0.9425085292914046, 0.9362670029301522, 0.9544426850176303, 0.9407207493748144, 0.9390630876730207,'..b'6807635995],\n+      [0.01505600521072376],\n+      [0.002882972930821469],\n+      [0.005825499564956559],\n+      [0.005903516300053625],\n+      [0.02274367484777162],\n+      [0.004355032959661065],\n+      [0.008800028033786655],\n+      [0.008917880494044083],\n+      [0.06143704799998004],\n+      [0.01176416611541591],\n+      [0.0237713451467977],\n+      [0.0240896977132239],\n+      [0.0262463496652457],\n+      [0.005025736545566795],\n+      [0.01015528996016046],\n+      [0.01029129247081581],\n+      [0.02256567087217255],\n+      [0.004320948178469224],\n+      [0.008731154380523415],\n+      [0.00884808446536876],\n+      [0.034087812394816],\n+      [0.006527245376826684],\n+      [0.01318932435908579],\n+      [0.01336595951511979],\n+      [0.01178447449308153],\n+      [0.02413128354890585],\n+      [0.0262916585104385],\n+      [0.005034412430862358],\n+      [0.01017282094890981],\n+      [0.01030905823951665],\n+      [0.004328407393735399],\n+      [0.008746226896405704],\n+      [0.008863358836639712],\n+      [0.03414665788893068],\n+      [0.006538513303760755],\n+      [0.0132120929750345],\n+      [0.01338903305480236] \n+      ],\n+     "Log Likelihood":-3176.217093030587,\n+     "Rate Distributions":{\n+       "non-synonymous/synonymous rate ratio for *test*":        [\n+[0.231720664913378, 1] \n+        ]\n+      },\n+     "display order":1,\n+     "estimated parameters":28\n+    },\n+   "Nucleotide GTR":{\n+     "AIC-c":6683.094244259512,\n+     "Equilibrium frequencies":      [\n+[0.4044507575757576],\n+      [0.1662878787878788],\n+      [0.2095643939393939],\n+      [0.2196969696969697] \n+      ],\n+     "Log Likelihood":-3320.503280793639,\n+     "Rate Distributions":{\n+       "Substitution rate from nucleotide A to nucleotide C":0.1815218119292989,\n+       "Substitution rate from nucleotide A to nucleotide G":1,\n+       "Substitution rate from nucleotide A to nucleotide T":0.08433479296165451,\n+       "Substitution rate from nucleotide C to nucleotide G":0.1775869374255926,\n+       "Substitution rate from nucleotide C to nucleotide T":1.184827038118431,\n+       "Substitution rate from nucleotide G to nucleotide T":0.04288318031815568\n+      },\n+     "display order":0,\n+     "estimated parameters":21\n+    },\n+   "Unconstrained model":{\n+     "AIC-c":6391.938375474656,\n+     "Log Likelihood":-3163.666348620609,\n+     "Rate Distributions":{\n+       "Test":{\n+         "0":{\n+           "omega":0.1490485393860332,\n+           "proportion":0.02658974082564446\n+          },\n+         "1":{\n+           "omega":0.169565784922581,\n+           "proportion":0.9641001770967245\n+          },\n+         "2":{\n+           "omega":15.09387961636913,\n+           "proportion":0.009310082077631076\n+          }\n+        }\n+      },\n+     "display order":2,\n+     "estimated parameters":32\n+    }\n+  },\n+ "input":{\n+   "file name":"/tmp/tmpfqvtss8b/job_working_directory/000/5/working/./input.nex",\n+   "number of sequences":8,\n+   "number of sites":440,\n+   "partition count":1,\n+   "trees":{\n+     "0":"((((D_CD_83_ELI,D_CD_83_NDK)Node3,D_UG_94_94UG114)Node2,D_CD_84_84ZR085)Node1,B_US_83_RF,((B_FR_83_HXB2,B_US_86_JRFL)Node10,B_US_90_WEAU160)Node9)"\n+    }\n+  },\n+ "test results":{\n+   "LRT":8.711160142067456,\n+   "p-value":0.006417496108925114\n+  },\n+ "tested":{\n+   "0":{\n+     "B_FR_83_HXB2":"test",\n+     "B_US_83_RF":"test",\n+     "B_US_86_JRFL":"test",\n+     "B_US_90_WEAU160":"test",\n+     "D_CD_83_ELI":"test",\n+     "D_CD_83_NDK":"test",\n+     "D_CD_84_84ZR085":"test",\n+     "D_UG_94_94UG114":"test",\n+     "Node1":"test",\n+     "Node10":"test",\n+     "Node2":"test",\n+     "Node3":"test",\n+     "Node9":"test"\n+    }\n+  },\n+ "timers":{\n+   "Constrained BUSTED model fitting":{\n+     "order":3,\n+     "timer":23\n+    },\n+   "Overall":{\n+     "order":0,\n+     "timer":133\n+    },\n+   "Preliminary model fitting":{\n+     "order":1,\n+     "timer":4\n+    },\n+   "Unconstrained BUSTED model fitting":{\n+     "order":2,\n+     "timer":105\n+    }\n+  }\n+}\n\\ No newline at end of file\n'
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/busted-out2-2.json
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/busted-out2-2.json Tue Apr 27 18:06:02 2021 +0000
[
b"@@ -0,0 +1,1730 @@\n+#NEXUS\n+\n+BEGIN TAXA;\n+\tDIMENSIONS NTAX = 8;\n+\tTAXLABELS\n+\t\t'B_FR_83_HXB2' 'B_US_83_RF' 'B_US_86_JRFL' 'B_US_90_WEAU160' 'D_CD_83_ELI' 'D_CD_83_NDK' 'D_CD_84_84ZR085' 'D_UG_94_94UG114' ;\n+END;\n+\n+BEGIN CHARACTERS;\n+\tDIMENSIONS NCHAR = 1320;\n+\tFORMAT\n+\t\tDATATYPE = DNA\n+\t\tGAP=-\n+\t\tMISSING=?\n+\t;\n+\n+MATRIX\n+\t'B_FR_83_HXB2'    CCCATTAGCCCTATTGAGACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTTAAACAATGGCCATTGACAGAAGAAAAAATAAAAGCATTAGTAGAAATTTGTACAGAGATGGAAAAGGAAGGGAAAATTTCAAAAATTGGGCCTGAAAATCCATACAATACTCCAGTATTTGCCATAAAGAAAAAAGACAGTACTAAATGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGACTTCTGGGAAGTTCAATTAGGAATACCACATCCCGCAGGGTTAAAAAAGAAAAAATCAGTAACAGTACTGGATGTGGGTGATGCATATTTTTCAGTTCCCTTAGATGAAGACTTCAGGAAGTATACTGCATTTACCATACCTAGTATAAACAATGAGACACCAGGGATTAGATATCAGTACAATGTGCTTCCACAGGGATGGAAAGGATCACCAGCAATATTCCAAAGTAGCATGACAAAAATCTTAGAGCCTTTTAGAAAACAAAATCCAGACATAGTTATCTATCAATACATGGATGATTTGTATGTAGGATCTGACTTAGAAATAGGGCAGCATAGAACAAAAATAGAGGAGCTGAGACAACATCTGTTGAGGTGGGGACTTACCACACCAGACAAAAAACATCAGAAAGAACCTCCATTCCTTTGGATGGGTTATGAACTCCATCCTGATAAATGGACAGTACAGCCTATAGTGCTGCCAGAAAAAGACAGCTGGACTGTCAATGACATACAGAAGTTAGTGGGGAAATTGAATTGGGCAAGTCAGATTTACCCAGGGATTAAAGTAAGGCAATTATGTAAACTCCTTAGAGGAACCAAAGCACTAACAGAAGTAATACCACTAACAGAAGAAGCAGAGCTAGAACTGGCAGAAAACAGAGAGATTCTAAAAGAACCAGTACATGGAGTGTATTATGACCCATCAAAAGACTTAATAGCAGAAATACAGAAGCAGGGGCAAGGCCAATGGACATATCAAATTTATCAAGAGCCATTTAAAAATCTGAAAACAGGAAAATATGCAAGAATGAGGGGTGCCCACACTAATGATGTAAAACAATTAACAGAGGCAGTGCAAAAAATAACCACAGAAAGCATAGTAATATGGGGAAAGACTCCTAAATTTAAACTGCCCATACAAAAGGAAACATGGGAAACATGGTGGACAGAGTATTGGCAAGCCACCTGGATTCCTGAGTGGGAGTTTGTTAATACCCCTCCCTTAGTGAAATTATGGTACCAGTTAGAGAAAGAACCCATAGTAGGAGCAGAAACCTTC\n+\t'B_US_83_RF'      CCCATTAGTCCTATTGAAACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTTAAACAATGGCCATTGACAGAGGAAAAAATAAAAGCATTGGTAGAAATTTGTACAGAAATGGAAAAGGAAGGAAAAATTTCCAAAATTGGGCCTGAAAATCCATACAATACTCCAGTATTTGCCATAAAGAAAAAAGACAGTACTAAATGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAGAACTCAAGACTTCTGGGAAGTTCAGTTAGGAATACCACATCCTGCAGGGTTAAAAAAGAAGAAATCAGTAACAGTATTGGATGTGGGTGATGCATATTTTTCAGTTCCCTTAGATAAAGAGTTCAGGAAGTATACTGCATTTACCATACCTAGTATAAACAATGAAACACCACGGATTAGATATCAGTACAATGTGCTTCCACAAGGGTGGAAAGGATCACCAGCAATATTCCAAAGTAGTATGACAAAAATCTTAGAGCCTTTTAAAAAACAAAATCCAGAAATAGTTATCTATCAATACATGGATGATTTGTATGTAGGATCTGATTTAGAAATAGGGCAGCATAGAATAAAAATAGAGGAACTGAGAGAACATCTGTTAAAGTGGGGGTTTACCACACCGGACAAGAAACATCAGAAAGAACCTCCATTTCTTTGGATGGGTTATGAACTCCATCCTGATAAATGGACAGTACAGCCTATAGTGCTGCCAGAAAAAGACAGCTGGACTGTCAATGACATACAGAAGTTAGTGGGAAAATTGAATTGGGCAAGTCAGATTTATGCAGGGATTAAAGTAAAGCAATTATGTAAACTCCTTAGGGGAACCAAAGCACTAACAGAAGTAGTACAACTAACAAAAGAAGCAGAGCTAGAACTGGCAGAAAATAGGGAGATTCTAAAAGAACCAGTACATGGAGTGTATTATGACCCATCAAAAGACTTAATAGCAGAAATACAGAAGCAGGGGCAAGGCCAATGGACATACCAAATTTATCAAGAGCCATTTAAAAACCTGAAAACAGGAAAGTATGCAAGAATGAGGGGTGCCCACACTAATGATGTAAAACAATTAACAGAGGCAGTACAAAAAGTAGCCACAGAAAGCATAGTAATATGGGGAAAGACTCCTAAATTTAAACTACCCATACAAAAAGAAACATGGGAGGCATGGTGGACAGAGTATTGGCAAGCCACCTGGATTCCTGAGTGGGAGTTTGTCAATACCCCTCCCTTAGTAAAATTGTGGTACCAGTTAGAAAAAGAACCCATAATAGGAGCAGAAACTTTC\n+\t'B_US_86_JRFL'    CCCATTAGTCCTATTGAAACTGTACCAGTAAAATTAAAGCCAGGAATGGATGGCCCAAAAGTCAAACAATGGCCATTGACAGAAGAAAAAATAAAAGCATTAGTAGAAATTTGTACAGAAATGGAAAAGGAAGGGAAAATTTCAAAAATTGGGCCTGAAAATCCATACAATACTCCAGTATTTGCCATAAAGAAAAAGGACAGTACTAAATGGAGAAAATTAGTAGATTTCAGAGAACTTAATAAGAAAACTCAAGACTTCTGGGAAGTTCAATTAGGAATACCACATCCCGCAGGGTTAAAAAAGAGAAAATCAGTAACAGTACTGGATGTGGGTGATGCATATTTTTCAGTTCCCTTAGATAAAGACTTCAGGAAATATACTGCATTTACCATACCTAGTATAAACAATGAGACACCAGGGATTAGGTATCAGTACAATGTGCTTCCGCAGGGATGGAAAGGATCACCAGCAATATTCCAAAGTAGCATGACAAAAATCTTAGAGCCTTTTAGAAAACAAAATCCAGACATAATTATCTATCAATACATGGATGATTTGTATGTAGGATCTGACTTAGAGATAGGGCAGCATAGAGCAAAAATAGAGGAATTGAGACAACATCTGTTGAGGTGGGGGTTTACCACACCAGACAAAAAACATCAGAAAGAACCTCCATTCCTTTGGATGGGTTATGAACTCCATCCTGACAAATGGACAGTACAGCCTATAGTGCTGCCAGAAAAAGACAGCTGGACTGTCAATGACATACAGAAGTTAGTGGGAAAATTAAATTGGGCAAGTCAGATTTACGCAGGGATTAAAGTAAAGCAATTATGTAAACTCCTTAGGGGAACCAAAGCACTAACAGAAGTAATACCACTAACAGAAGAAGCAGAGCTAGAACTGGCAGAAAACAGGGAGATTCTAAAAGAGCCAGTACATGGAGTGTATTATGACCCATCAAAAGACT"..b'nt_3[60][31]:=0.1712746635706107*t*busted.test.theta_CT*busted.test.omega3;\n+busted.test_Q_component_3[60][47]:=0.2567033972679326*t*busted.test.theta_GT*busted.test.omega3;\n+busted.test_Q_component_3[60][49]:=0.4256503784977386*t*busted.test.theta_AT*busted.test.omega3;\n+busted.test_Q_component_3[60][53]:=0.1818409092376876*t*busted.test.theta_CT*busted.test.omega3;\n+busted.test_Q_component_3[60][56]:=0.1563402972714242*t*busted.test.theta_GT*busted.test.omega3;\n+busted.test_Q_component_3[60][57]:=0.5074830879161225*t*busted.test.theta_AT*busted.test.omega3;\n+busted.test_Q_component_3[60][58]:=0.09717451507454297*t*busted.test.theta_CT;\n+busted.test_Q_component_3[60][59]:=0.1963563685387479*t*busted.test.theta_GT*busted.test.omega3;\n+\n+\n+busted.test_pi={\n+{0.07535940345366139} \n+{0.01443006409728228} \n+{0.02915821068672351} \n+{0.02954870567752823} \n+{0.03219408024958758} \n+{0.006164627375802654} \n+{0.01245659774045292} \n+{0.01262342001471309} \n+{0.0276793164843978} \n+{0.005300125700761717} \n+{0.01070973633984303} \n+{0.01085316415297191} \n+{0.0418125103783387} \n+{0.008006395714089126} \n+{0.01617817990958604} \n+{0.0163948426631036} \n+{0.0409913146415633} \n+{0.007849150478922548} \n+{0.01586044121725868} \n+{0.01607284872448998} \n+{0.01751178502783737} \n+{0.003353213651232053} \n+{0.006775694789785172} \n+{0.006866436807635995} \n+{0.01505600521072376} \n+{0.002882972930821469} \n+{0.005825499564956559} \n+{0.005903516300053625} \n+{0.02274367484777162} \n+{0.004355032959661065} \n+{0.008800028033786655} \n+{0.008917880494044083} \n+{0.06143704799998004} \n+{0.01176416611541591} \n+{0.0237713451467977} \n+{0.0240896977132239} \n+{0.0262463496652457} \n+{0.005025736545566795} \n+{0.01015528996016046} \n+{0.01029129247081581} \n+{0.02256567087217255} \n+{0.004320948178469224} \n+{0.008731154380523415} \n+{0.00884808446536876} \n+{0.034087812394816} \n+{0.006527245376826684} \n+{0.01318932435908579} \n+{0.01336595951511979} \n+{0.01178447449308153} \n+{0.02413128354890585} \n+{0.0262916585104385} \n+{0.005034412430862358} \n+{0.01017282094890981} \n+{0.01030905823951665} \n+{0.004328407393735399} \n+{0.008746226896405704} \n+{0.008863358836639712} \n+{0.03414665788893068} \n+{0.006538513303760755} \n+{0.0132120929750345} \n+{0.01338903305480236} \n+};\n+Model busted.test=(\'Exp(busted.test_Q_component_1)*1*busted.test.bsrel_mixture_aux_0+Exp(busted.test_Q_component_2)*1*(1-busted.test.bsrel_mixture_aux_0)*busted.test.bsrel_mixture_aux_1+Exp(busted.test_Q_component_3)*1*(1-busted.test.bsrel_mixture_aux_0)*(1-busted.test.bsrel_mixture_aux_1)\',busted.test_pi,EXPLICIT_FORM_MATRIX_EXPONENTIAL);\n+TRY_NUMERIC_SEQUENCE_MATCH=0;\n+ACCEPT_ROOTED_TREES=0;\n+\n+UseModel (busted.test);\n+Tree cuWdXMfP.tree_id_0=((((D_CD_83_ELI,D_CD_83_NDK)Node3,D_UG_94_94UG114)Node2,D_CD_84_84ZR085)Node1,B_US_83_RF,((B_FR_83_HXB2,B_US_86_JRFL)Node10,B_US_90_WEAU160)Node9);\n+\n+cuWdXMfP.tree_id_0.D_CD_83_ELI.t=0.1377903766097694;\n+cuWdXMfP.tree_id_0.D_CD_83_NDK.t=0.06607109645594529;\n+cuWdXMfP.tree_id_0.Node3.t=0.0890081009446355;\n+cuWdXMfP.tree_id_0.D_UG_94_94UG114.t=0.4322474154485661;\n+cuWdXMfP.tree_id_0.Node2.t=0.01436667144151429;\n+cuWdXMfP.tree_id_0.D_CD_84_84ZR085.t=0.2152052309993384;\n+cuWdXMfP.tree_id_0.Node1.t=0.1867555551194544;\n+cuWdXMfP.tree_id_0.B_US_83_RF.t=0.2093873697545922;\n+cuWdXMfP.tree_id_0.B_FR_83_HXB2.t=0.08552069057797769;\n+cuWdXMfP.tree_id_0.B_US_86_JRFL.t=0.1313117411090917;\n+cuWdXMfP.tree_id_0.Node10.t=0.01181372524396516;\n+cuWdXMfP.tree_id_0.B_US_90_WEAU160.t=0.1600546391700751;\n+cuWdXMfP.tree_id_0.Node9.t=0.03506352054252541;SetParameter (DEFER_CONSTRAINT_APPLICATION, 1, 0);\n+SetParameter (DEFER_CONSTRAINT_APPLICATION, 0, 0);\n+\n+DataSet busted.codon_data = ReadDataFile(USE_NEXUS_FILE_DATA);\n+DataSetFilter busted.filter.default = CreateFilter(busted.codon_data,3,"0-1319","4,5,7,6,1,0,2,3","TAA,TAG,TGA");\n+ASSUME_REVERSIBLE_MODELS=1;\n+USE_LAST_RESULTS=0;\n+LikelihoodFunction cuWdXMfP.likelihoodFunction = (busted.filter.default,cuWdXMfP.tree_id_0);\n+\n+END;\n\\ No newline at end of file\n'
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/gard-out1.json
--- a/test-data/gard-out1.json Tue Apr 20 19:54:45 2021 +0000
+++ b/test-data/gard-out1.json Tue Apr 27 18:06:02 2021 +0000
[
b'@@ -3,316 +3,53 @@\n    "authors":"Sergei L Kosakovsky Pond",\n    "citation":"**Automated Phylogenetic Detection of Recombination Using a Genetic Algorithm**, _Mol Biol Evol 23(10), 1891\xe2\x80\x931901",\n    "contact":"spond@temple.edu",\n-   "info":"GARD : Genetic Algorithms for Recombination Detection. Implements a heuristic\\n    approach to screening alignments of sequences for recombination, by using the CHC genetic algorithm to search for phylogenetic\\n    incongruence among different partitions of the data. The number of partitions is determined using a step-up procedure, while the\\n    placement of breakpoints is searched for with the GA. The best fitting model (based on c-AIC) is returned; and additional post-hoc\\n    tests run to distinguish topological incongruence from rate-variation.",\n+   "info":"GARD : Genetic Algorithms for Recombination Detection. Implements a heuristic\\n    approach to screening alignments of sequences for recombination, by using the CHC genetic algorithm to search for phylogenetic\\n    incongruence among different partitions of the data. The number of partitions is determined using a step-up procedure, while the\\n    placement of breakpoints is searched for with the GA. The best fitting model (based on c-AIC) is returned; and additional post-hoc\\n    tests run to distinguish topological incongruence from rate-variation. v0.2 adds  and spooling results to JSON after each breakpoint search conclusion",\n    "requirements":"A sequence alignment.",\n-   "version":"0.1"\n+   "version":"0.2"\n   },\n- "baselineScore":8884.998566352424,\n+ "baselineScore":1450.968340148184,\n+ "bestModelAICc":1450.968340148184,\n  "breakpointData":{\n-   "0":{\n-     "bps":      [\n-[1, 461] \n-      ],\n-     "tree":"(((((((TREESPARROW_HENAN_1_2004:0.01973996359641645,TREESPARROW_HENAN_2_2004:0)Node6:0.006601594605919903,CHICKEN_HEBEI_326_2005:0.02877904543700591)Node5:0.002960697265622527,SWINE_ANHUI_2004:0.01220336915594569)Node4:0.004301817483024902,((DUCK_GUANGZHOU_20_2005:0.006496379151101181,(GOOSE_SHANTOU_2216_2005:0.008666763019733412,PEREGRINEFALCON_HK_D0028_2004:0.006494907750172948)Node14:0)Node12:0.006529386248252737,CK_HK_WF157_2003:0.008749819170042668)Node11:0.006517868346508392)Node3:0.004369843971142481,((TREESPARROW_HENAN_3_2004:0.004310725480928182,((CHICKEN_HONGKONG_915_97:0,(HONGKONG_97_98:0.007234718327388172,HONGKONG_538_97:0)Node24:0)Node22:0,(GOOSE_HONGKONG_W355_97:0,DUCK_HONGKONG_Y283_97:0)Node27:0.002160209380619512)Node21:0.02651640051975179)Node19:0.002209468499615626,TREESPARROW_HENAN_4_2004:0.008707376458939551)Node18:0)Node2:0.01527162191938894,MALLARD_VIETNAM_16_2003:0.00216300836609909)Node1:0,HUMAN_VIETNAM_CL105_2005:0.01089634139979296,(((DUCK_VIETNAM_376_2005:0.002165747961920766,CHICKEN_THAILAND_KANCHANABURI_CK_160_2005:0.01092964139649115)Node35:0.002165571630784742,VIETNAM_3062_2004:0)Node34:0.002160659783816526,DUCK_VIETNAM_272_2005:0.006501625005083654)Node33:0)"\n-    },\n-   "1":{\n-     "bps":      [\n-[462, 650] \n-      ],\n-     "tree":"((((((TREESPARROW_HENAN_1_2004:0,TREESPARROW_HENAN_3_2004:0.01064277225673711)Node5:0.00531987705161906,TREESPARROW_HENAN_2_2004:0)Node4:0.005231591103466107,CHICKEN_HEBEI_326_2005:0.05043188754547531)Node3:0.005435427857192481,CK_HK_WF157_2003:0)Node2:0.005367407450692448,((TREESPARROW_HENAN_4_2004:0.005403333058615208,(((CHICKEN_HONGKONG_915_97:0,(GOOSE_HONGKONG_W355_97:0,DUCK_HONGKONG_Y283_97:0)Node18:0.005303893633648769)Node16:0,HONGKONG_97_98:0)Node15:0,HONGKONG_538_97:0)Node14:0.0160450339521957)Node12:0.01065583332597595,((DUCK_GUANGZHOU_20_2005:0.02156617952480121,GOOSE_SHANTOU_2216_2005:0.005315754078746889)Node24:0,PEREGRINEFALCON_HK_D0028_2004:0.005310132566451705)Node23:0)Node11:0.005324661683629481)Node1:0.01075596871747344,(SWINE_ANHUI_2004:0.01612941661930485,DUCK_VIETNAM_272_2005:0.005332918418299072)Node28:0.005343017103483476,((HUMAN_VIETNAM_CL105_2005:0.01615610619083483,DUCK_VIETNAM_376_2005:0.01622595235090508)Node32:0.005403715853761346,(VI'..b'326_2005:0.05043188754547531)Node3:0.005435427857192481,CK_HK_WF157_2003:0)Node2:0.005367407450692448,((TREESPARROW_HENAN_4_2004:0.005403333058615208,(((CHICKEN_HONGKONG_915_97:0,(GOOSE_HONGKONG_W355_97:0,DUCK_HONGKONG_Y283_97:0)Node18:0.005303893633648769)Node16:0,HONGKONG_97_98:0)Node15:0,HONGKONG_538_97:0)Node14:0.0160450339521957)Node12:0.01065583332597595,((DUCK_GUANGZHOU_20_2005:0.02156617952480121,GOOSE_SHANTOU_2216_2005:0.005315754078746889)Node24:0,PEREGRINEFALCON_HK_D0028_2004:0.005310132566451705)Node23:0)Node11:0.005324661683629481)Node1:0.01075596871747344,(SWINE_ANHUI_2004:0.01612941661930485,DUCK_VIETNAM_272_2005:0.005332918418299072)Node28:0.005343017103483476,((HUMAN_VIETNAM_CL105_2005:0.01615610619083483,DUCK_VIETNAM_376_2005:0.01622595235090508)Node32:0.005403715853761346,(VIETNAM_3062_2004:0,(MALLARD_VIETNAM_16_2003:0.01066500623689393,CHICKEN_THAILAND_KANCHANABURI_CK_160_2005:0)Node37:0.005316372444723995)Node35:0)Node31:0)"\n-    },\n-   "2":{\n-     "newickString":"(((((((TREESPARROW_HENAN_1_2004:0.001954008187578416,TREESPARROW_HENAN_3_2004:0.003896316686452931)Node6:0.001938506890703636,TREESPARROW_HENAN_4_2004:0)Node5:0.005890390504717593,((CHICKEN_HEBEI_326_2005:0.0118071209154035,(SWINE_ANHUI_2004:0,TREESPARROW_HENAN_2_2004:0.001939618450410721)Node13:0.01383276182388132)Node11:0.009974147725569997,((CHICKEN_HONGKONG_915_97:0,GOOSE_HONGKONG_W355_97:0)Node17:0.001941847153562197,((DUCK_HONGKONG_Y283_97:0,HONGKONG_538_97:0)Node21:0,HONGKONG_97_98:0.005840445012127324)Node20:0.001942715072940801)Node16:0.01987446960589405)Node10:0.005913850464875198)Node4:0.005886750382502456,(((DUCK_GUANGZHOU_20_2005:0.003893604042162395,GOOSE_SHANTOU_2216_2005:0.003905435534305824)Node27:0.001941929040845811,PEREGRINEFALCON_HK_D0028_2004:0.003899037680959845)Node26:0.001951590047476092,CK_HK_WF157_2003:0.005864286435734769)Node25:0.001925781157666145)Node3:0.009803928769016717,DUCK_VIETNAM_272_2005:0.01580156952635168)Node2:0,CHICKEN_THAILAND_KANCHANABURI_CK_160_2005:0.003890598353658968)Node1:0,((HUMAN_VIETNAM_CL105_2005:0.003905371265795895,DUCK_VIETNAM_376_2005:0.0137665342243657)Node35:0.001932867193851192,MALLARD_VIETNAM_16_2003:0.003887459835724038)Node34:0,VIETNAM_3062_2004:0.00194102772512571)"\n-    },\n-   "3":{\n-     "newickString":"(TREESPARROW_HENAN_1_2004:0,((TREESPARROW_HENAN_3_2004:0.005713571370335148,TREESPARROW_HENAN_4_2004:0)Node3:0.00189607429178233,(CHICKEN_HEBEI_326_2005:0.005712316903324198,(SWINE_ANHUI_2004:0.005701998731046575,(((CHICKEN_HONGKONG_915_97:0,HONGKONG_97_98:0)Node12:0.001895621658683188,(GOOSE_HONGKONG_W355_97:0,HONGKONG_538_97:0)Node15:0)Node11:0.01534803500036186,((DUCK_HONGKONG_Y283_97:0.0001941011375241317,(((HUMAN_VIETNAM_CL105_2005:0.001928158212296621,DUCK_VIETNAM_376_2005:0.003824815472812549)Node23:0.001900556734363341,(VIETNAM_3062_2004:0.001945118440015516,MALLARD_VIETNAM_16_2003:0)Node26:0)Node22:0,(CHICKEN_THAILAND_KANCHANABURI_CK_160_2005:0.00380825293990655,DUCK_VIETNAM_272_2005:0.005749159769097136)Node29:0.001884324589239026)Node21:0.006136201619836659)Node19:0.006956231732180768,(((DUCK_GUANGZHOU_20_2005:0.003782134609270667,GOOSE_SHANTOU_2216_2005:0.005745814969410696)Node34:0.005728903718080784,PEREGRINEFALCON_HK_D0028_2004:0)Node33:0.001883474579815808,CK_HK_WF157_2003:0.005783922057126199)Node32:0.003686703559509533)Node18:0.002359028413665987)Node10:0.003781653021946557)Node8:0.003836413220384994)Node6:0.001894177236914572)Node2:0,TREESPARROW_HENAN_2_2004:0.00767492225272676)"\n-    }\n+   "newickString":"((BRO1_1:0,BRO2:0.007408926175080497)Node1:0,(BRO3_1:0.002455524404855022,BRO4:0.002457415365886074)Node4:0,(((BRO5:0.002464479112430666,(BRO6:0.004940850926103086,BRO7_1:0)Node11:0)Node9:0.002460542727794925,(BRO9:0.004931390913247289,((BRO10:0.002933649980980302,BRO11_1:0.01453409151831228)Node17:0.002436921831192641,BRO12:0.006987081661112447)Node16:0.005417113782285573)Node14:0)Node8:0,BRO8:0.002462026756884335)Node7:0)"\n   }\n }\n\\ No newline at end of file\n'
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/gard-out1.nex
--- a/test-data/gard-out1.nex Tue Apr 20 19:54:45 2021 +0000
+++ b/test-data/gard-out1.nex Tue Apr 27 18:06:02 2021 +0000
b
@@ -0,0 +1,37 @@
+#NEXUS
+
+BEGIN TAXA;
+ DIMENSIONS NTAX = 12;
+ TAXLABELS
+ 'BRO1_1' 'BRO2' 'BRO3_1' 'BRO4' 'BRO5' 'BRO6' 'BRO7_1' 'BRO9' 'BRO10' 'BRO11_1' 'BRO12' 'BRO8' ;
+END;
+
+BEGIN CHARACTERS;
+ DIMENSIONS NCHAR = 428;
+ FORMAT
+ DATATYPE = DNA
+ GAP=-
+ MISSING=?
+ ;
+
+MATRIX
+ 'BRO1_1'  ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACATCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO2'    ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACATCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAACTGCAACCCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTGCATAAAGCCATTTACCGTACATAGCACATTATA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO3_1'  ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACATCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAGCTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO4'    ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACATCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCATACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO5'    ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACAGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACACCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO6'    ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAATACATAAAAACCC-AATTC--ACACCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO7_1'  ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACACCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO9'    ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AACCC--ACATCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAATTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO10'   ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACATCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAATCAACCTTCAA-CTATCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGCACATAAAGCCATTTATCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO11_1' ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACTACCTGTAGTACATAAAAACCC-AACCC--ACATCAAAA----CCCTGCCCC-CATGCTTACAAGCAAGTACAGCAATCAACCTTCAA-CTGTCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGCACATAAAGTCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO12'   ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACAGTACCATAAA-TACTTGACTACCTGTAGTACATAAAAACCC-AATCC--ACATCAAAA----CCCCCTCCC-CATGCTTACAAGCAAGTACAGCAATCAACCTTCAA-CTATCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGGATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC
+ 'BRO8'    ATTCTAATTTAAACTATTCT-CTGTTCTTTCATGGGGAAGCAGATTTGGGTACCACCCAAGTATTGACTCACCCATCAACAACCGCTATGTATTTCGTACATTAC-TGCCAGCCACCATGAATATTGTACGGTACCATAAA-TACTTGACCACCTGTAGTACATAAAAACCC-AATCC--ACATCAAAA----CCCCCCCCC-CATGCTTACAAGCAAGTACAGCAACCAACCCTCAA-CTATCACACATCAACTGCAACTCCAAAGCCAC-CCCTCACCCAC-TAGAATACCAACAAACCTACCCACCC-TTAACAGTACATAGTACATAAAGCCATTTACCGTACATAGCACATTACA-GTCAAATCCCTTCTCGTCCCC-ATGG-ATGACCCCCC-TCAGAT-AGGGGTCCCTTGACCACCATCC;
+END;
+
+BEGIN ASSUMPTIONS;
+ CHARSET span_1 = 1-428;
+END;
+BEGIN TREES;
+ TREE tree_1 = (((BRO1_1:1e-10,BRO2:0.007446048662499141)Node2:1.699082282683926e-05,(((BRO5:0.002291051777868712,(BRO6:0.004965711012353363,BRO7_1:1e-10)Node9:0.0001955595843124865)Node7:0.002191267352110732,(BRO9:0.004469023573866088,((BRO10:0.002870586428367251,BRO11_1:0.01462820627459318)Node15:0.001454707803761366,BRO12:0.006027959855484483)Node14:0.005661145911839864)Node12:0.0003039715029625697)Node6:0.0002803682303139691,BRO8:0.002464847898104024)Node5:1.386828770154629e-05)Node1:3.347353150116037e-06,BRO3_1:0.002467834453766899,BRO4:0.002462875913900023);
+END;
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/relax-in2.nex
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/relax-in2.nex Tue Apr 27 18:06:02 2021 +0000
[
b"@@ -0,0 +1,86 @@\n+#NEXUS\n+\n+[\n+Generated by HYPHY 0.9920070619beta(MP) for MacOS(Carbon) on Thu Jun 21 16:59:04 2007\n+\n+]\n+\n+BEGIN TAXA;\n+\tDIMENSIONS NTAX = 10;\n+\tTAXLABELS\n+\t\t'Human' 'Chimp' 'Baboon' 'RhMonkey' 'Cow' 'Pig' 'Horse' 'Cat' 'Mouse' 'Rat' ;\n+END;\n+\n+BEGIN CHARACTERS;\n+\tDIMENSIONS NCHAR = 561;\n+\tFORMAT\n+\t\tDATATYPE = DNA\n+\n+\t\tGAP=-\n+\t\tMISSING=?\n+\t;\n+\n+MATRIX\n+\t'Human'     AAAGAGATTACGAATGCCTTGGAAACCTGGGGTGCCTTGGGTCAGGACATCAACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGACGATATAAAATGGGAAAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAAGAGACTTTCAAGGAAAAAGATACATATAAGCTATTTAAAAATGGAACTCTGAAAATTAAGCAT---CTGAAGACCGATGATCAGGATATCTACAAGGTATCAATATATGATACAAAAGGAAAAAATGTGTTGGAAAAAATATTTGATTTGAAGATTCAAGAGAGGGTCTCAAAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAGGTAATGAATGGAACTGACCCCGAATTAAACCTGTATCAAGATGGGAAACATCTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCAAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGTGTCGAGCCTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'Chimp'     GAAGAGATTACGAATGCCTTGGAAACCTGGGGTGCCTTGGGTCAGGACATCAACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGACGATATAAAATGGGAAAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAAGAGACTTTCAAGGAAAAAGATACATATAAGCTATTTAAAAATGGAACTCTGAAAATTAAGCAT---CTGAAGACCGATGATCAGGATATCTACAAGGTATCAATATATGATACAAAAGGAAAAAATGTGTTGGAAAAAATATTTGATTTGAAGATTCAAGAGAGGGTCTCAAAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAGGTAATGAATGGAACTGACCCCGAATTAAACCTGTATCAAGATGGGAAACATCTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCAAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGTGTCGAGCCTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'Baboon'    AAAGAGATTAGGAATGCTTTGGAAACCTGGGGAGCGCTGGGTCAGGACATCGACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGATGATATAAAATGGGAGAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAGGAGACTTTCGAGGAAAAAGATGCATATAAGCTATTTAAAAACGGAACTCTGAAAATTAAGCAT---CTGAAGATCCATGATCAGGATAGCTACAAGGTATCAATATACGATACAAAAGGAAAAAATGTGTTGGAAAAAACATTTGATTTGAAGATTCAAGAGAGGGTCTCAGAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAAGTAATGAATGGAACTGACCCCGAATTAAACCTGTATCAAGATGGGAAACATCTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCGAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGGATGGAGACTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'RhMonkey'  AAAGAGATTAGGAATGCTTTGGAAACCTGGGGAGCGCTGGGTCAGGACATCGACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGATGATATAAGATGGGAAAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAGGAGACTTTCGAGGAAAAAGATGCATATAAGCTATTTAAAAACGGAACTCTGAAAAYTAAGCAT---CTGAAGATCCATGATCAGGATAGCTACAAGGTATCAATATACGATACAAAAGGAAAAAATGTGTTGGAAAAAACATTTGATTTGAAGATTCAAGAGAGGGTCTCAGAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAAGTAATGAATGGAACTGRCCCCGAATTAAACCTGTATCAAGATGGGAAACATGTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCGAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGGATGGAGACTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'Cow'       ------------GAAAGCATTGTCGTCTGGGGTGCCCTGGATCATGACCTCAACCTGGACATTCCTGGTTTTCCAAGAAGTGATATAGTGGCAGATATAAAATGGAACAGA------AACAAAAACAAGATTGCACGAATAAAGAAAGATATGCCACTTCACAATGAAATGGACAAATATGATATGTTTACAAATGGAACTCTGAAAATTAAAACT---CTGATGAGAAACGATAGTGGTCTCTATGAGGTAGAGGTTTATGATTCAAATGGAGTAAACCTACTGAGCAAAAAATTTGATTTGAAGATTCAAGAGATGCTCTCAGGACCTGAAATTAACTGGATCTGTACCAACAGAACTGTGAGCTGCAAGGTAGAAAATGGAAGTAATCCTAAATTACAACTGTTTTTAAATACGACCCGTGTCAAACAAGATCATGGGAAGCTCATCACCTACACGTGGAACACCAGATGGAATAAAACATTCAAGTGCGTGGCGAGTAACCATGTCGATAGCAAAGTCAGCATAGAGATCGCCGTGTGTCCAGATGAAGGTCTGGAT\n+\t'Pig'       ---------------ACTGAGGTTGTCTGGGGCATCGTGGATCAAGACATCAACCTGGACATTCCTGAACTTTCAAAACATGATAACGTAGATCATATACGATGGCAGAAG------AATGAAAACAAGATCGCAGAATTTAAAAAAAACAAAGAAACTCACCCTGTGAAAGACACATACATGATGTTACCAAATGGAACTCTGAGAATTAAAGAT---CTGAAGAGAGATGATGAGGGTATCTACAAGGTAACTGTCTATGCTACGGATGGAAAACACATGCTGGAGAGAAAATTTGATTTGCCGATTCTAGATGGGGTCTCAAAACCTGTAATCTCCTGGAGCTGTGCCGACAAAACGGTGACCTGTGAGGTAGCAGAAGGAAGTGACCCTAAGTTAAAACTGTATGTAAATAAGTCCACTGCCAGAGAAGGTCGTCAGAAGGTCATCCTGTGGAAGTGGAACACCAAATGGAGCACATTATTCAAGTGTGTGGCCAGTAACAACGCCAGTGAGCAAATCAGCATGGTGACCATCAGTTGTACGGGGCAAGGTCTGGAT\n+\t'Horse'     ------------AAGAATATCACCATCTTGGGTGCCCTGGAACGTGATATCAACCTGGACATTCCTGCTTTTCAAATGAGTGAGCATGTAGAAGATATACAATGGAGCAAA------GGAAAAACCAAGATTGCAAAATTCAAAAATGGCAGTATGA"..b'GCAAATGATGATATCGTCTGGGGTACCCTGGGTCAGGACATCAACCTGGACATTCCTGATTCTCAA---GGGATTAATATAGATGATATACACTGGGAAAAA------GGCAAGAAGAAGGTGGCGAGGTTCCAAATTAGCAACAAGCCTAAGAATCCAGATGAAAAATATAATGTGTCAATGAATGGAACTCTGAAAATTAAACAT---CTGATGCTAGAAGACTGCGATACCTACAAGGTTGTTATATACGATAAGGATGGAAAGAATGTGTTGGATAAAACATTTCAGCTGAAGATTCAAGAGAAGGTCTCAACGCCTAACATCGACTGGAATTGTATCAACAAAACCCTGGTCTGTAAGGTATCAAATGGAACAGACCCTGAATTAAAACTGTACGTAAATGGGACCAGTATCAAGCCCGTTTCTTCGAAGTTCAGCACATACAGGTTTATAAACAAGCAGAAGATATTAGTCAACTGCACGGCAGAAAACAAAGTCAGCAAGGAAAGCGACGTGAAGATGATCACTTGTTCAGAGAAGGGTCTGGAC\n+\t\'Mouse\'     ---------AGAGACAATGAGACCATCTGGGGTGTCTTGGGTCATGGCATCACCCTGAACATCCCCAACTTTCAAATGACTGATGATATTGATGAGGTGCGATGGGTAAGG------AGGGGCACCCTGGTCGCAGAGTTTAAAAGGAAGAAGCCACCTTTTTTGATATCAGAAACGTATGAGGTCTTAGCAAACGGATCCCTGAAGATAAAGAAGCCGATGATGAGAAACGACAGTGGCACCTATAATGTAATGGTGTATGGCACAAATGGGATGACTAGGCTGGAGAAGGACCTGGACGTGAGGATTCTGGAGAGGGTCTCAAAGCCCATGATCCACTGGGAATGCCCCAACACAACCCTGACCTGTGCGGTCTTGCAAGGGACAGATTTTGAACTGAAGCTGTATCAAGGGGAAACACTACTCAATAGTCTCCCCCAGAAGAACATGAGTTACCAGTGG---ACCAACCTGAACGCACCATTCAAGTGTGAGGCGATAAACCCGGTCAGCAAGGAGTCTAAGATGGAAGTTGTTAACTGTCCAGAGAAAGGTCTGTCC\n+\t\'Rat\'       ---------AGAGACAGTGGGACCGTCTGGGGTGCCCTGGGTCATGGCATCAACCTGAACATCCCTAACTTTCAAATGACTGATGATATTGATGAGGTGCGATGGGAGAGG------GGGAGCACCCTGGTTGCCGAGTTTAAAAGGAAGATGAAGCCTTTTTTGAAATCGGGAGCATTTGAGATCTTAGCAAATGGAGACTTGAAGATAAAGAAT---CTGACAAGAGATGACAGTGGCACCTATAATGTAACGGTATACAGCACAAATGGGACACGTATCCTGGACAAGGCACTGGACTTGAGGATTCTAGAGATGGTCTCAAAGCCGATGATCTACTGGGAGTGCAGCAACGCAACCCTGACCTGTGAGGTCTTGGAAGGAACAGATGTTGAACTAAAGCTGTACCAAGGAAAGGAGCATCTCAGGAGCCTCCGTCAGAAGACCATGAGTTACCAGTGG---ACCAACCTGAGAGCACCGTTTAAGTGCAAGGCGGTAAACAGGGTCAGCCAGGAGTCTGAGATGGAAGTTGTCAACTGTCCAGAGAAAGGTCTGCCC;\n+END;\n+\n+BEGIN TREES;\n+\tTREE tree = ((((Pig:0.147969,Cow:0.213430):0.085099,Horse:0.165787,Cat:0.264806):0.058611,((RhMonkey{PR}:0.002015,Baboon{PR}:0.003108){PR}:0.022733,(Human{PR}:0.004349,Chimp{PR}:0.000799){PR}:0.011873){PR}:0.101856){PR}:0.340802,Rat:0.050958,Mouse:0.097950);\n+END;\n+\n+BEGIN HYPHY;\n+\n+DATA_PANEL_SOURCE_PATH="";\n+DataSet CD2 = ReadDataFile (USE_NEXUS_FILE_DATA);\n+DataSetFilter CD2_part = CreateFilter (CD2,3,"0-560","(null)","TAA,TAG,TGA");\n+Tree CD2_tree=((((Pig,Cow)Node3,Horse,Cat)Node2,((RhMonkey,Baboon)Node9,(Human,Chimp)Node12)Node8)Node1,Rat,Mouse);\n+OpenDataPanel(CD2,"(null)","0,0,0,0,2,0,10","MG94xREV_3x4,6,0,14443523,CD2_tree",CD2_LF);\n+\n+global CD2_part_Shared_AC=0.5710714109234631;\n+global CD2_part_Shared_AT=0.2756961373814887;\n+global CD2_part_Shared_CG=0.4456308760366715;\n+global CD2_part_Shared_CT=0.9454291760502562;\n+global CD2_part_Shared_GT=0.3224080259059006;\n+CD2_tree.Node12.synRate=0.08168540352543764;\n+CD2_tree.Node8.nonSynRate=0.2938890030949185;\n+CD2_tree.Human.synRate=1e-30;\n+CD2_tree.Baboon.nonSynRate=0;\n+CD2_tree.Baboon.synRate=0.01716650886630832;\n+CD2_tree.Node12.nonSynRate=0.02962786133931161;\n+CD2_tree.Chimp.synRate=1e-30;\n+CD2_tree.Chimp.nonSynRate=0.005834136861061837;\n+CD2_tree.Node9.nonSynRate=0.0433794394076959;\n+CD2_tree.Human.nonSynRate=0;\n+CD2_tree.Node9.synRate=0.1139808847086709;\n+CD2_tree.Node8.synRate=0.1586154947684889;\n+CD2_tree.Rat.synRate=0.1533251749656354;\n+CD2_tree.Mouse.nonSynRate=0.2332402459120736;\n+CD2_tree.Mouse.synRate=0.4604267094040578;\n+CD2_tree.Rat.nonSynRate=0.1585006932449144;\n+CD2_tree.Node1.nonSynRate=0.5097883291072822;\n+CD2_tree.Node1.synRate=1.227004762734591;\n+CD2_tree.RhMonkey.synRate=0;\n+CD2_tree.Node3.synRate=0.1877287558837553;\n+CD2_tree.Pig.nonSynRate=0.4826366809736332;\n+CD2_tree.Horse.nonSynRate=0.518627673056701;\n+CD2_tree.Cow.nonSynRate=0.681831503936737;\n+CD2_tree.Cat.synRate=0.4672127128189641;\n+CD2_tree.Node3.nonSynRate=0.2648957141533021;\n+CD2_tree.Horse.synRate=0.430860769501357;\n+CD2_tree.Pig.synRate=0.3721230153643161;\n+CD2_tree.Node2.synRate=0.2137858662796737;\n+CD2_tree.RhMonkey.nonSynRate=0.01161863878954486;\n+CD2_tree.Cat.nonSynRate=0.7203773170684773;\n+CD2_tree.Node2.nonSynRate=0.1379570235930709;\n+CD2_tree.Cow.synRate=0.3646221318370894;\n+\n+END;\n'
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/relax-out1.json
--- a/test-data/relax-out1.json Tue Apr 20 19:54:45 2021 +0000
+++ b/test-data/relax-out1.json Tue Apr 27 18:06:02 2021 +0000
[
b'@@ -1,486 +1,174 @@\n {\n- "analysis":{\n-   "info":"RELAX (a random effects test of selection relaxation) uses a random effects branch-site model framework to test whether a set of \'Test\' branches evolves under relaxed selection relative to a set of \'Reference\' branches (R), as measured by the relaxation parameter (K).\\n                                                Version 2.1 adds a check for stability in K estimates to try to mitigate convergence problems",\n-   "version":"2.1",\n-   "citation":"RELAX: Detecting Relaxed Selection in a Phylogenetic Framework (2015). Mol Biol Evol 32 (3): 820-832",\n-   "authors":"Sergei L Kosakovsky Pond, Ben Murrell, Steven Weaver and Temple iGEM / UCSD viral evolution group",\n-   "contact":"spond@temple.edu",\n-   "requirements":"in-frame codon alignment and a phylogenetic tree, with at least two groups of branches defined using the {} notation (one group can be defined as all unlabeled branches)"\n-  },\n- "input":{\n-   "file name":"/home/galaxystaging/galaxy/database/jobs_directory/000/79/working/relax_input.fa",\n-   "number of sequences":10,\n-   "number of sites":187,\n-   "partition count":1,\n-   "trees":{\n-     "0":"((((Pig:0.147969,Cow:0.21343)Node3:0.085099,Horse:0.165787,Cat:0.264806)Node2:0.058611,((RhMonkey:0.002015,Baboon:0.003108)Node9:0.022733,(Human:0.004349,Chimp:0.000799)Node12:0.011873)Node8:0.101856)Node1:0.340802,Rat:0.050958,Mouse:0.09795)"\n-    }\n-  },\n- "fits":{\n-   "Nucleotide GTR":{\n-     "Log Likelihood":-3531.964032193222,\n-     "estimated parameters":24,\n-     "AIC-c":7112.142925621896,\n-     "Equilibrium frequencies":[\n-      [0.3563279857397504],\n-      [0.1837789661319073],\n-      [0.2402852049910873],\n-      [0.2196078431372549]\n-      ],\n-     "Rate Distributions":{\n-       "Substitution rate from nucleotide A to nucleotide C":0.5465757107930838,\n-       "Substitution rate from nucleotide A to nucleotide G":1,\n-       "Substitution rate from nucleotide A to nucleotide T":0.264131123736109,\n-       "Substitution rate from nucleotide C to nucleotide G":0.4861812992450388,\n-       "Substitution rate from nucleotide C to nucleotide T":1.017887589835735,\n-       "Substitution rate from nucleotide G to nucleotide T":0.3021773696987065\n-      },\n-     "display order":0\n-    },\n-   "MG94xREV with separate rates for branch sets":{\n-     "Log Likelihood":-3460.931903299487,\n-     "estimated parameters":32,\n-     "AIC-c":6987.013507197776,\n-     "Equilibrium frequencies":[\n-      [0.04785168648208005],\n-      [0.03911346079312019],\n-      [0.04842054060317301],\n-      [0.04069665666515385],\n-      [0.01817010744778396],\n-      [0.01485205303123071],\n-      [0.01838611113045968],\n-      [0.01545321970309948],\n-      [0.0190138898893362],\n-      [0.01554175184587827],\n-      [0.01923992433904378],\n-      [0.01617083546230154],\n-      [0.02796925133965702],\n-      [0.02286176927319543],\n-      [0.02830174586718859],\n-      [0.02378714529482007],\n-      [0.01687870337957799],\n-      [0.01379647304851128],\n-      [0.01707935503225849],\n-      [0.01435491300079667],\n-      [0.006409133648843443],\n-      [0.005238757839512795],\n-      [0.006485324528555602],\n-      [0.005450807083377185],\n-      [0.006706760641638982],\n-      [0.005482034985409045],\n-      [0.006786489669196017],\n-      [0.005703931360295163],\n-      [0.009865581170012171],\n-      [0.008064021368173151],\n-      [0.009982861811889655],\n-      [0.008390428826966174],\n-      [0.03287895020631181],\n-      [0.02687490503171972],\n-      [0.033269810543681],\n-      [0.0279627207822247],\n-      [0.01248470225271473],\n-      [0.01020486315060615],\n-      [0.01263311863778268],\n-      [0.01061792547971484],\n-      [0.01306446616325414],\n-      [0.01067875601941112],\n-      [0.01321977470017749],\n-      [0.01111100011404158],\n-      [0.01921770557551369],\n-      [0.01570834862514353],\n-      [0.01944616295744266],\n-      [0.01634417558076203],\n-      [0.01975257768318274],\n-      [0.02055210293142003],\n-      [0.00917603432671'..b'    "omega":2.19971015101232,\n+           "proportion":0.4620912608347691\n+          }\n+        },\n+       "Test":{\n+         "0":{\n+           "omega":1.054840660215593e-44,\n+           "proportion":0.4504132080156107\n+          },\n+         "1":{\n+           "omega":1.054840660215593e-44,\n+           "proportion":0.08749553114962015\n+          },\n+         "2":{\n+           "omega":6.350765030542244,\n+           "proportion":0.4620912608347691\n+          }\n+        }\n+      },\n+     "display order":2,\n+     "estimated parameters":36\n+    },\n+   "RELAX null":{\n+     "AIC-c":6946.046672507471,\n+     "Log Likelihood":-3437.336313352972,\n+     "Rate Distributions":{\n+       "Reference":{\n+         "0":{\n+           "omega":0.1150006822953211,\n+           "proportion":0.542405144031307\n+          },\n+         "1":{\n+           "omega":0.2401508033814633,\n+           "proportion":0.08223353569636478\n+          },\n+         "2":{\n+           "omega":3.69263720947284,\n+           "proportion":0.3753613202723282\n+          }\n+        },\n+       "Test":{\n+         "0":{\n+           "omega":0.1150006822953211,\n+           "proportion":0.542405144031307\n+          },\n+         "1":{\n+           "omega":0.2401508033814633,\n+           "proportion":0.08223353569636478\n+          },\n+         "2":{\n+           "omega":3.69263720947284,\n+           "proportion":0.3753613202723282\n+          }\n+        }\n+      },\n+     "display order":3,\n+     "estimated parameters":35\n+    },\n+   "RELAX partitioned descriptive":{\n+     "AIC-c":6937.04707353565,\n+     "Log Likelihood":-3427.626871923648,\n+     "Rate Distributions":{\n+       "Reference":{\n+         "0":{\n+           "omega":0.1943306885168281,\n+           "proportion":0.3359458135547653\n+          },\n+         "1":{\n+           "omega":0.2607558478587815,\n+           "proportion":0.4332550437039296\n+          },\n+         "2":{\n+           "omega":4.134271224300374,\n+           "proportion":0.2307991427413052\n+          }\n+        },\n+       "Test":{\n+         "0":{\n+           "omega":0,\n+           "proportion":0.4753488415962343\n+          },\n+         "1":{\n+           "omega":0,\n+           "proportion":0.03858318603767072\n+          },\n+         "2":{\n+           "omega":6.017221959238516,\n+           "proportion":0.486067972366095\n+          }\n+        }\n+      },\n+     "display order":5,\n+     "estimated parameters":40\n+    }\n+  },\n+ "input":{\n+   "file name":"/tmp/tmpbauyxvdv/job_working_directory/000/3/working/./input.fasta",\n+   "number of sequences":10,\n+   "number of sites":187,\n+   "partition count":1,\n+   "trees":{\n+     "0":"((((Pig:0.147969,Cow:0.21343)Node3:0.08509899999999999,Horse:0.165787,Cat:0.264806)Node2:0.058611,((RhMonkey:0.002015,Baboon:0.003108)Node9:0.022733,(Human:0.004349,Chimp:0.000799)Node12:0.011873)Node8:0.101856)Node1:0.340802,Rat:0.050958,Mouse:0.09795)"\n+    }\n+  },\n+ "test results":{\n+   "LRT":16.20300116404087,\n+   "p-value":5.690389451429301e-05,\n+   "relaxation or intensification parameter":2.344938789085876\n+  },\n+ "tested":{\n+   "0":{\n+     "Baboon":"Reference",\n+     "Cat":"Test",\n+     "Chimp":"Reference",\n+     "Cow":"Test",\n+     "Horse":"Test",\n+     "Human":"Reference",\n+     "Mouse":"Reference",\n+     "Node1":"Reference",\n+     "Node12":"Reference",\n+     "Node2":"Reference",\n+     "Node3":"Reference",\n+     "Node8":"Reference",\n+     "Node9":"Reference",\n+     "Pig":"Reference",\n+     "Rat":"Reference",\n+     "RhMonkey":"Reference"\n+    }\n+  },\n+ "timers":{\n+   "General descriptive model fitting":{\n+     "order":2,\n+     "timer":132\n+    },\n+   "Overall":{\n+     "order":0,\n+     "timer":324\n+    },\n+   "Preliminary model fitting":{\n+     "order":1,\n+     "timer":8\n+    },\n+   "RELAX alternative model fitting":{\n+     "order":3,\n+     "timer":65\n+    },\n+   "RELAX null model fitting":{\n+     "order":4,\n+     "timer":66\n+    },\n+   "RELAX partitioned descriptive":{\n+     "order":5,\n+     "timer":53\n+    }\n   }\n }\n\\ No newline at end of file\n'
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/relax-out2-1.json
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/relax-out2-1.json Tue Apr 27 18:06:02 2021 +0000
[
b'@@ -0,0 +1,525 @@\n+{\n+ "Synonymous site-posteriors":  [\n+[0.2717348381268465, 0.5635158143073038, 0.5580958985234344, 0.4419267392327456, 0.1775513314393873, 0.01303903284026798, 0.00125071980641974, 0.2376006215377312, 0.8837687913991302, 0.697488732253998, 0.76440666757549, 0.1859131352278527, 0.4676839084660329, 0.5320161033357819, 0.4513687349306337, 0.8514990846254493, 0.9101555767748879, 0.9112112654445478, 0.9624478616388369, 0.9472735837633918, 0.8715894173292342, 0.9248557783401865, 0.3852870058193366, 0.7504386709677574, 0.7275694213063354, 0.1979029930878776, 0.2585255313006612, 0.426812371619399, 0.03564414482172501, 0.3292941826787388, 0.1427864928730872, 0.8324822227281018, 0.9311167984186458, 0.1209539167468769, 0.8726642109588801, 1.798577785061191e-05, 0.584102863529593, 0.5747348035616634, 0.5814406651776627, 0.02251840852563282, 0.1607002230623373, 0.499402751628698, 0.9327364905381349, 0.1164578922551683, 0.7414294602970922, 0.04004321943657974, 0.1252919499200681, 0.7016909689123084, 0.3099147090319344, 0.1594271385761372, 0.001721088532959928, 0.0002106542138058178, 0.7670245595319366, 0.01765606757145387, 0.08528617939935416, 0.01185048796937841, 0.008218049898468071, 0.294285815926654, 0.4088819101657737, 0.1819814701794264, 0.06741006997792491, 0.1648869485561294, 0.6641637464726119, 0.1696421815033135, 0.5337345685031449, 0.9785329157950278, 0.6089117121305408, 0.9214198863478258, 0.7040511824729966, 0.8748810706893103, 0.7603207301908244, 0.1587648197710778, 0.28009015626025, 0.9473769917649164, 0.03667497185975328, 0.6646432746680839, 0.005707942217110617, 0.4172795698582943, 0.002426641321239203, 0.9036479890682162, 0.8165937919459629, 0.5380686752231248, 0.8348947378138256, 0.8825623566861087, 0.003194607014896987, 0.1034323185116836, 0.1251183560640448, 0.6731739154042389, 0.02049494631105097, 0.7272701824195644, 0.9376578101739967, 0.03161216349993828, 0.08060038603464466, 0.05713766811791781, 0.7236799166647727, 0.03049686563145122, 0.4831167820690108, 0.3259152163804719, 0.9355990871535625, 0.2598255955335131, 0.649043721535527, 0.1395744847880144, 0.8655465087299661, 0.5535165112383126, 0.8759601400318501, 0.1331847979385494, 0.9385179832315519, 0.9590294197009215, 0.1177725177949566, 0.5310119897564728, 0.04436932195281112, 0.877273324169318, 0.07309972300058731, 0.8726642109588801, 0.0823189228969842, 0.674523501474146, 0.4093792912564727, 0.8837869026711356, 0.1208432028324516, 0.7875716367935024, 0.9342090047974666, 0.7281976470761374, 0.6335560301454017, 0.2642816632435872, 0.7601746562441594, 0.3971532429331721, 0.07967254820654134, 0.9205611063156326, 0.4126590346311199, 0.4718488374006657, 0.2324587146851173, 0.2530879180545083, 0.6735650752750825, 0.6778765556227805, 0.9393601651000962, 0.1757530496502532, 0.3931715696771638, 0.7913247469383754, 0.05849180455884406, 0.007027787474985664, 0.002980098295172601, 0.188566078048928, 0.02837364230845984, 0.008078062694748547, 0.4539148488391334, 0.2698860362225194, 0.03514662617223248, 0.8484422182215356, 0.5225500090930265, 0.8755046006217432, 0.01843684138990366, 0.6656855504934923, 0.7055653078864003, 0.236066638653311, 0.1567944657161159, 0.8781825877891227, 0.4362059051173987, 0.163327088806649, 0.004489982631214889, 0.03278376720971195, 0.03145887036879619, 0.6158371404365816, 0.8645373560796752, 0.6351698115277034, 0.04852978149651612, 0.2563374460202912, 0.07868312568368593, 0.9507071563814276, 0.007940480658871889, 0.9455733391892027, 0.4761025393051341, 0.06329409664274364, 0.486108446887078, 0.6554607928409875, 0.3294403859657761, 0.02718898835050054, 0.2632599441236705, 0.09462189267290307, 0.2747185708643275, 0.01971954733787938, 0.9591589529807978, 0.445066915720239, 0.7380186230528726, 0.3686112558749921, 0.9765948224615983, 0.9830372572797934, 0.3428744949444181],\n+  [0.6626894690463668, 0.4168101145890417, 0.4227716136375318, 0.5460469376257381, 0.7920802349484392, 0.9148946082433389, 0.6025519507878634, 0.757695330723024'..b'+           "proportion":0\n+          },\n+         "2":{\n+           "omega":1.00002621592498,\n+           "proportion":0.6185938027364208\n+          }\n+        }\n+      },\n+     "display order":2,\n+     "estimated parameters":41\n+    },\n+   "RELAX null":{\n+     "AIC-c":6910.777215627525,\n+     "Log Likelihood":-3414.491942969586,\n+     "Rate Distributions":{\n+       "Reference":{\n+         "0":{\n+           "omega":0,\n+           "proportion":0.04600287725570656\n+          },\n+         "1":{\n+           "omega":0,\n+           "proportion":0.3432466763494891\n+          },\n+         "2":{\n+           "omega":1.896477292815783,\n+           "proportion":0.6107504463948045\n+          }\n+        },\n+       "Synonymous site-to-site rates":{\n+         "0":{\n+           "proportion":0.2641095265325168,\n+           "rate":0.2744837794048847\n+          },\n+         "1":{\n+           "proportion":0.6330324513232739,\n+           "rate":0.9729462345446451\n+          },\n+         "2":{\n+           "proportion":0.1028580221442092,\n+           "rate":3.029415427453047\n+          }\n+        },\n+       "Test":{\n+         "0":{\n+           "omega":0,\n+           "proportion":0.04600287725570656\n+          },\n+         "1":{\n+           "omega":0,\n+           "proportion":0.3432466763494891\n+          },\n+         "2":{\n+           "omega":1.896477292815783,\n+           "proportion":0.6107504463948045\n+          }\n+        }\n+      },\n+     "display order":3,\n+     "estimated parameters":40\n+    },\n+   "RELAX partitioned descriptive":{\n+     "AIC-c":6909.42794543134,\n+     "Log Likelihood":-3408.579104294618,\n+     "Rate Distributions":{\n+       "Reference":{\n+         "0":{\n+           "omega":0.05035684771569529,\n+           "proportion":0.3805224005150026\n+          },\n+         "1":{\n+           "omega":0.07013486723985318,\n+           "proportion":0.003462124163764896\n+          },\n+         "2":{\n+           "omega":2.344944921144001,\n+           "proportion":0.6160154753212325\n+          }\n+        },\n+       "Test":{\n+         "0":{\n+           "omega":0,\n+           "proportion":0.03160127773820932\n+          },\n+         "1":{\n+           "omega":0,\n+           "proportion":0.4079248249859382\n+          },\n+         "2":{\n+           "omega":1.11269394570475,\n+           "proportion":0.5604738972758525\n+          }\n+        }\n+      },\n+     "display order":5,\n+     "estimated parameters":45\n+    }\n+  },\n+ "input":{\n+   "file name":"/tmp/tmpbauyxvdv/job_working_directory/000/5/working/./input.nex",\n+   "number of sequences":10,\n+   "number of sites":187,\n+   "partition count":1,\n+   "trees":{\n+     "0":"((((PIG:0.147969,COW:0.21343)Node3:0.08509899999999999,HORSE:0.165787,CAT:0.264806)Node2:0.058611,((RHMONKEY:0.002015,BABOON:0.003108)Node9:0.022733,(HUMAN:0.004349,CHIMP:0.000799)Node12:0.011873)Node8:0.101856)Node1:0.340802,RAT:0.050958,MOUSE:0.09795)"\n+    }\n+  },\n+ "test results":{\n+   "LRT":11.47546300090198,\n+   "p-value":0.0007052107640865835,\n+   "relaxation or intensification parameter":2.986519568748539e-05\n+  },\n+ "tested":{\n+   "0":{\n+     "BABOON":"Test",\n+     "CAT":"Reference",\n+     "CHIMP":"Test",\n+     "COW":"Reference",\n+     "HORSE":"Reference",\n+     "HUMAN":"Test",\n+     "MOUSE":"Reference",\n+     "Node1":"Test",\n+     "Node12":"Test",\n+     "Node2":"Reference",\n+     "Node3":"Reference",\n+     "Node8":"Test",\n+     "Node9":"Test",\n+     "PIG":"Reference",\n+     "RAT":"Reference",\n+     "RHMONKEY":"Test"\n+    }\n+  },\n+ "timers":{\n+   "General descriptive model fitting":{\n+     "order":2,\n+     "timer":127\n+    },\n+   "Overall":{\n+     "order":0,\n+     "timer":645\n+    },\n+   "Preliminary model fitting":{\n+     "order":1,\n+     "timer":8\n+    },\n+   "RELAX alternative model fitting":{\n+     "order":3,\n+     "timer":178\n+    },\n+   "RELAX null model fitting":{\n+     "order":4,\n+     "timer":115\n+    },\n+   "RELAX partitioned descriptive":{\n+     "order":5,\n+     "timer":216\n+    }\n+  }\n+}\n\\ No newline at end of file\n'
b
diff -r b6509de48a62 -r 70fab57f3c05 test-data/relax-out2-2.json
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/relax-out2-2.json Tue Apr 27 18:06:02 2021 +0000
[
b"@@ -0,0 +1,3427 @@\n+#NEXUS\n+\n+BEGIN TAXA;\n+\tDIMENSIONS NTAX = 10;\n+\tTAXLABELS\n+\t\t'HUMAN' 'CHIMP' 'BABOON' 'RHMONKEY' 'COW' 'PIG' 'HORSE' 'CAT' 'MOUSE' 'RAT' ;\n+END;\n+\n+BEGIN CHARACTERS;\n+\tDIMENSIONS NCHAR = 561;\n+\tFORMAT\n+\t\tDATATYPE = DNA\n+\t\tGAP=-\n+\t\tMISSING=?\n+\t;\n+\n+MATRIX\n+\t'HUMAN'    AAAGAGATTACGAATGCCTTGGAAACCTGGGGTGCCTTGGGTCAGGACATCAACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGACGATATAAAATGGGAAAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAAGAGACTTTCAAGGAAAAAGATACATATAAGCTATTTAAAAATGGAACTCTGAAAATTAAGCAT---CTGAAGACCGATGATCAGGATATCTACAAGGTATCAATATATGATACAAAAGGAAAAAATGTGTTGGAAAAAATATTTGATTTGAAGATTCAAGAGAGGGTCTCAAAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAGGTAATGAATGGAACTGACCCCGAATTAAACCTGTATCAAGATGGGAAACATCTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCAAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGTGTCGAGCCTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'CHIMP'    GAAGAGATTACGAATGCCTTGGAAACCTGGGGTGCCTTGGGTCAGGACATCAACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGACGATATAAAATGGGAAAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAAGAGACTTTCAAGGAAAAAGATACATATAAGCTATTTAAAAATGGAACTCTGAAAATTAAGCAT---CTGAAGACCGATGATCAGGATATCTACAAGGTATCAATATATGATACAAAAGGAAAAAATGTGTTGGAAAAAATATTTGATTTGAAGATTCAAGAGAGGGTCTCAAAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAGGTAATGAATGGAACTGACCCCGAATTAAACCTGTATCAAGATGGGAAACATCTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCAAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGTGTCGAGCCTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'BABOON'   AAAGAGATTAGGAATGCTTTGGAAACCTGGGGAGCGCTGGGTCAGGACATCGACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGATGATATAAAATGGGAGAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAGGAGACTTTCGAGGAAAAAGATGCATATAAGCTATTTAAAAACGGAACTCTGAAAATTAAGCAT---CTGAAGATCCATGATCAGGATAGCTACAAGGTATCAATATACGATACAAAAGGAAAAAATGTGTTGGAAAAAACATTTGATTTGAAGATTCAAGAGAGGGTCTCAGAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAAGTAATGAATGGAACTGACCCCGAATTAAACCTGTATCAAGATGGGAAACATCTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCGAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGGATGGAGACTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'RHMONKEY' AAAGAGATTAGGAATGCTTTGGAAACCTGGGGAGCGCTGGGTCAGGACATCGACTTGGACATTCCTAGTTTTCAAATGAGTGATGATATTGATGATATAAGATGGGAAAAAACTTCAGACAAGAAAAAGATTGCACAATTCAGAAAAGAGAAGGAGACTTTCGAGGAAAAAGATGCATATAAGCTATTTAAAAACGGAACTCTGAAAAYTAAGCAT---CTGAAGATCCATGATCAGGATAGCTACAAGGTATCAATATACGATACAAAAGGAAAAAATGTGTTGGAAAAAACATTTGATTTGAAGATTCAAGAGAGGGTCTCAGAACCAAAGATCTCCTGGACTTGTATCAACACAACCCTGACCTGTGAAGTAATGAATGGAACTGRCCCCGAATTAAACCTGTATCAAGATGGGAAACATGTAAAA---CTTTCTCAGAGGGTCATCACACACAAGTGGACCACCAGCCTGAGTGCGAAATTCAAGTGCACAGCAGGGAACAAAGTCAGCAAGGAATCCAGGATGGAGACTGTCAGCTGTCCAGAGAAAGGTCTGGAC\n+\t'COW'      ------------GAAAGCATTGTCGTCTGGGGTGCCCTGGATCATGACCTCAACCTGGACATTCCTGGTTTTCCAAGAAGTGATATAGTGGCAGATATAAAATGGAACAGA------AACAAAAACAAGATTGCACGAATAAAGAAAGATATGCCACTTCACAATGAAATGGACAAATATGATATGTTTACAAATGGAACTCTGAAAATTAAAACT---CTGATGAGAAACGATAGTGGTCTCTATGAGGTAGAGGTTTATGATTCAAATGGAGTAAACCTACTGAGCAAAAAATTTGATTTGAAGATTCAAGAGATGCTCTCAGGACCTGAAATTAACTGGATCTGTACCAACAGAACTGTGAGCTGCAAGGTAGAAAATGGAAGTAATCCTAAATTACAACTGTTTTTAAATACGACCCGTGTCAAACAAGATCATGGGAAGCTCATCACCTACACGTGGAACACCAGATGGAATAAAACATTCAAGTGCGTGGCGAGTAACCATGTCGATAGCAAAGTCAGCATAGAGATCGCCGTGTGTCCAGATGAAGGTCTGGAT\n+\t'PIG'      ---------------ACTGAGGTTGTCTGGGGCATCGTGGATCAAGACATCAACCTGGACATTCCTGAACTTTCAAAACATGATAACGTAGATCATATACGATGGCAGAAG------AATGAAAACAAGATCGCAGAATTTAAAAAAAACAAAGAAACTCACCCTGTGAAAGACACATACATGATGTTACCAAATGGAACTCTGAGAATTAAAGAT---CTGAAGAGAGATGATGAGGGTATCTACAAGGTAACTGTCTATGCTACGGATGGAAAACACATGCTGGAGAGAAAATTTGATTTGCCGATTCTAGATGGGGTCTCAAAACCTGTAATCTCCTGGAGCTGTGCCGACAAAACGGTGACCTGTGAGGTAGCAGAAGGAAGTGACCCTAAGTTAAAACTGTATGTAAATAAGTCCACTGCCAGAGAAGGTCGTCAGAAGGTCATCCTGTGGAAGTGGAACACCAAATGGAGCACATTATTCAAGTGTGTGGCCAGTAACAACGCCAGTGAGCAAATCAGCATGGTGACCATCAGTTGTACGGGGCAAGGTCTGGAT\n+\t'HORSE'    ------------AAGAATATCACCATCTTGGGTGCCCTGGAACGTGATATCAACCTGGACATTCCTGCTTTTCAAATGAGTGAGCATGTAGAAGATATACAATGGAGCAAA------GGAAAAACCAAGATTGCAAAATTCAAAAATGGCAGTATGACTTTCCAGAAAGATAAAACATACGAGGTATTAAAAAATGGAACTCTGAAAATTAAACAT---CTGGAGAGAATTCATGAAGGTACCTACAAGGTAGACGCATAT"..b'023*t*relax._shared_srv.rv_gdd*relax.test.theta_CT*relax.test.omega3;\n+relax.test_Q_component_3[60][56]:=0.1659124550130927*t*relax._shared_srv.rv_gdd*relax.test.theta_GT*relax.test.omega3;\n+relax.test_Q_component_3[60][57]:=0.2714859010999351*t*relax._shared_srv.rv_gdd*relax.test.theta_AT*relax.test.omega3;\n+relax.test_Q_component_3[60][58]:=0.2217183243979596*t*relax._shared_srv.rv_gdd*relax.test.theta_CT;\n+relax.test_Q_component_3[60][59]:=0.2755785876068215*t*relax._shared_srv.rv_gdd*relax.test.theta_GT*relax.test.omega3;\n+\n+\n+relax.test_pi={\n+{0.04876086397632926} \n+{0.03982224127745371} \n+{0.0494959405650269} \n+{0.04152830682371345} \n+{0.01810345074594235} \n+{0.01478480742075003} \n+{0.01837636270304883} \n+{0.01541821854326533} \n+{0.01892910005547028} \n+{0.01545910240515689} \n+{0.01921445878706794} \n+{0.01612140169177368} \n+{0.02829747763884322} \n+{0.02311011106415996} \n+{0.02872406592369408} \n+{0.02410019507229243} \n+{0.01663805943190806} \n+{0.01358804506433061} \n+{0.01688888041768248} \n+{0.01417018446135155} \n+{0.006177213955434187} \n+{0.005044834822351987} \n+{0.006270336287397383} \n+{0.00526096577332031} \n+{0.006458939937330408} \n+{0.005274916061901906} \n+{0.006556309326397028} \n+{0.005500904159606446} \n+{0.009655594186286731} \n+{0.007885574003572761} \n+{0.009801153568494354} \n+{0.008223407979974132} \n+{0.03327785448210915} \n+{0.02717750758119989} \n+{0.03377952261833742} \n+{0.02834184710178207} \n+{0.01235507229404208} \n+{0.01009019590244688} \n+{0.01254132667213156} \n+{0.01052248035033209} \n+{0.01291855364640413} \n+{0.01055038238273567} \n+{0.01311330258173711} \n+{0.01100238215235291} \n+{0.01931219560697892} \n+{0.0157719705998612} \n+{0.01960332955558842} \n+{0.01644767379420349} \n+{0.01965132601881859} \n+{0.02049322866376102} \n+{0.008933621043463411} \n+{0.007295949736387736} \n+{0.009068296583349953} \n+{0.007608523013862737} \n+{0.007628698224327205} \n+{0.009481877006094248} \n+{0.00795552712159237} \n+{0.01396413011303752} \n+{0.01140428846505027} \n+{0.01417464125436193} \n+{0.01189287130232129} \n+};\n+Model relax.test=(\'Exp(relax.test_Q_component_1)*1*relax.test.bsrel_mixture_aux_0+Exp(relax.test_Q_component_2)*1*(1-relax.test.bsrel_mixture_aux_0)*relax.test.bsrel_mixture_aux_1+Exp(relax.test_Q_component_3)*1*(1-relax.test.bsrel_mixture_aux_0)*(1-relax.test.bsrel_mixture_aux_1)\',relax.test_pi,EXPLICIT_FORM_MATRIX_EXPONENTIAL);\n+TRY_NUMERIC_SEQUENCE_MATCH=0;\n+ACCEPT_ROOTED_TREES=0;\n+Tree O_SJmzEX.tree_id_0=((((PIG{relax.reference},COW{relax.reference})Node3{relax.reference},HORSE{relax.reference},CAT{relax.reference})Node2{relax.reference},((RHMONKEY{relax.test},BABOON{relax.test})Node9{relax.test},(HUMAN{relax.test},CHIMP{relax.test})Node12{relax.test})Node8{relax.test})Node1{relax.test},RAT{relax.reference},MOUSE{relax.reference});\n+\n+O_SJmzEX.tree_id_0.PIG.t=0.4536172365685577;\n+O_SJmzEX.tree_id_0.COW.t=0.6103212983838633;\n+O_SJmzEX.tree_id_0.Node3.t=0.2484758643179366;\n+O_SJmzEX.tree_id_0.HORSE.t=0.4753747079725157;\n+O_SJmzEX.tree_id_0.CAT.t=0.6775352409130856;\n+O_SJmzEX.tree_id_0.Node2.t=0.1183276378467358;\n+O_SJmzEX.tree_id_0.RHMONKEY.t=0.01420592464486799;\n+O_SJmzEX.tree_id_0.BABOON.t=0.00485857082125907;\n+O_SJmzEX.tree_id_0.Node9.t=0.1015257435189256;\n+O_SJmzEX.tree_id_0.HUMAN.t=0;\n+O_SJmzEX.tree_id_0.CHIMP.t=0.006303088015069631;\n+O_SJmzEX.tree_id_0.Node12.t=0.05697754732758539;\n+O_SJmzEX.tree_id_0.Node8.t=0.4948729770972781;\n+O_SJmzEX.tree_id_0.Node1.t=1.445991686991002;\n+O_SJmzEX.tree_id_0.RAT.t=0.1315597576696033;\n+O_SJmzEX.tree_id_0.MOUSE.t=0.2736332745970835;SetParameter (DEFER_CONSTRAINT_APPLICATION, 1, 0);\n+SetParameter (DEFER_CONSTRAINT_APPLICATION, 0, 0);\n+\n+DataSet relax.codon_data = ReadDataFile(USE_NEXUS_FILE_DATA);\n+DataSetFilter relax.filter.default = CreateFilter(relax.codon_data,3,"0-560","5,4,6,7,3,2,0,1,9,8","TAA,TAG,TGA");\n+ASSUME_REVERSIBLE_MODELS=1;\n+USE_LAST_RESULTS=0;\n+LikelihoodFunction O_SJmzEX.likelihoodFunction = (relax.filter.default,O_SJmzEX.tree_id_0);\n+\n+END;\n\\ No newline at end of file\n'