| 1 | item:FastTree2 (nucleotide) |
|---|
| 2 | itemmeta:F |
|---|
| 3 | seqtype:N |
|---|
| 4 | itemmethod:(GOTO_LOCAL_DIR; \ |
|---|
| 5 | arb_sed 's/#/>/' <TmpInputFile | \ |
|---|
| 6 | arb_sed 's/%/>/' | \ |
|---|
| 7 | arb_sed 's/"/>/'| \ |
|---|
| 8 | arb_sed 's/@/>/' | \ |
|---|
| 9 | arb_sed '/^[^>]/ y/./-/' > fasttree_in; \ |
|---|
| 10 | RUN_IN_WINDOW(( ($ARBHOME/bin/FastTree $FAST $MODEL -nt fasttree_in > fasttree_out || STOP_ON_FAILURE(FastTree)) && \ |
|---|
| 11 | $DISPLAY_FUNC));)& |
|---|
| 12 | |
|---|
| 13 | |
|---|
| 14 | arg:DISPLAY_FUNC |
|---|
| 15 | argtype:choice_menu |
|---|
| 16 | arglabel:What to do with the tree? |
|---|
| 17 | argchoice:Load into ARB ('tree_fasttree_//'):$ARBHOME/bin/arb_read_tree tree_fasttree_$$ fasttree_out "PRG=FastTree$FAST FILTER=$FILTER MODEL=!$MODEL" |
|---|
| 18 | argchoice:Treetool:treetool fasttree_out |
|---|
| 19 | |
|---|
| 20 | |
|---|
| 21 | arg:FAST |
|---|
| 22 | argtype:choice_menu |
|---|
| 23 | arglabel:speed up the neighbor joining phase\n(recommended for >50,000 sequences) |
|---|
| 24 | argchoice:No: |
|---|
| 25 | argchoice:Yes:-fastest |
|---|
| 26 | |
|---|
| 27 | arg:MODEL |
|---|
| 28 | argtype:choice_menu |
|---|
| 29 | arglabel:ML Model |
|---|
| 30 | argchoice:Jukes-Cantor: |
|---|
| 31 | argchoice:GTR:-gtr |
|---|
| 32 | |
|---|
| 33 | |
|---|
| 34 | in:TmpInputFile |
|---|
| 35 | informat:flat |
|---|
| 36 | insave: |
|---|
| 37 | |
|---|
| 38 | //---------------------- |
|---|
| 39 | |
|---|
| 40 | #define PREPARE_MRBAYES_INFILE(INFILE) \ |
|---|
| 41 | ( \ |
|---|
| 42 | arb_sed -i '/LABELPOS/d' INFILE `#MrBayes does not understand the LABELPOS command, therefore it is removed.` && \ |
|---|
| 43 | arb_sed -i '/EQUATE/d' INFILE `#MrBayes does not support EQUATE macros, remove them.` && \ |
|---|
| 44 | arb_sed -i '/OPTIONS/,/;/d' INFILE `#MrBayes does not support the OPTIONS block, remove it. (delete everything between OPTIONS and ';')` \ |
|---|
| 45 | ) |
|---|
| 46 | |
|---|
| 47 | #define NEXUS2NEWICK(FROM,APPENDTO) ($ARBHOME/bin/arb_nexus2newick.awk FROM >> APPENDTO || STOP_ON_FAILURE(arb_nexus2newick)) |
|---|
| 48 | |
|---|
| 49 | item:MrBayes |
|---|
| 50 | itemmeta:B |
|---|
| 51 | seqtype:* |
|---|
| 52 | itemmethod:(GOTO_LOCAL_DIR; \ |
|---|
| 53 | GEN_NEXUS(mrbayes_in) && \ |
|---|
| 54 | PREPARE_MRBAYES_INFILE(mrbayes_in) && \ |
|---|
| 55 | ( \ |
|---|
| 56 | `# append the mrbayes command block to the nexus file`; \ |
|---|
| 57 | COMMAND="begin mrbayes; \ |
|---|
| 58 | set autoclose=yes nowarn=yes; \ |
|---|
| 59 | lset nst=$NST rates=$RATES ngammacat=$NGAMMACAT; \ |
|---|
| 60 | mcmc ngen=$NGEN nchains=$NCHAINS temp=$TEMP samplefreq=$SAMPLEFREQ burninfrac=$BURNINFRAC; \ |
|---|
| 61 | sump; \ |
|---|
| 62 | sumt Outputname=mrbayes_out Conformat=simple; \ |
|---|
| 63 | end;" \ |
|---|
| 64 | && echo $COMMAND >> mrbayes_in; \ |
|---|
| 65 | ) && \ |
|---|
| 66 | RUN_IN_WINDOW(( ($ARBHOME/bin/mb mrbayes_in || STOP_ON_FAILURE(MrBayes)) && \ |
|---|
| 67 | NEXUS2NEWICK(mrbayes_out.con.tre,mr_bayes_out.newick) && \ |
|---|
| 68 | $DISPLAY_FUNC));)& |
|---|
| 69 | |
|---|
| 70 | itemhelp:MrBayes.help |
|---|
| 71 | |
|---|
| 72 | arg:DISPLAY_FUNC |
|---|
| 73 | argtype:choice_menu |
|---|
| 74 | arglabel:What to do with the tree? |
|---|
| 75 | argchoice:Load into ARB ('tree_mrBayes_//'):$ARBHOME/bin/arb_read_tree tree_mrbayes_$$ mr_bayes_out.newick "PRG=MrBayes FILTER=$FILTER RATE=$RATES" |
|---|
| 76 | argchoice:Treetool:treetool mr_bayes_out.newick |
|---|
| 77 | |
|---|
| 78 | arg:NST |
|---|
| 79 | argtype:choice_menu |
|---|
| 80 | arglabel:Number of substitution types |
|---|
| 81 | argchoice:6:6 |
|---|
| 82 | argchoice:mixed:mixed |
|---|
| 83 | argchoice:1:1 |
|---|
| 84 | argchoice:2:2 |
|---|
| 85 | |
|---|
| 86 | arg:RATES |
|---|
| 87 | argtype:choice_menu |
|---|
| 88 | arglabel:Model for among-site rate variation |
|---|
| 89 | argchoice:A proportion of the sites are invariable:propinv |
|---|
| 90 | argchoice:No rate variation across sites:equal |
|---|
| 91 | argchoice:Gamma-distributed rates across sites:gamma |
|---|
| 92 | argchoice:Autocorrelated rates across sites:adgamma |
|---|
| 93 | argchoice:A proportion of the sites are invariable while the rate for the remaining sites are drawn from a gamma distribution:invgamma |
|---|
| 94 | |
|---|
| 95 | |
|---|
| 96 | arg:NGAMMACAT |
|---|
| 97 | argtype:text |
|---|
| 98 | arglabel:Number of rate categories for the gamma distribution\n(Range [0 .. n]) |
|---|
| 99 | argtext:4 |
|---|
| 100 | |
|---|
| 101 | arg:NGEN |
|---|
| 102 | argtype:text |
|---|
| 103 | arglabel:Number of cycles for the MCMC algorithm |
|---|
| 104 | argtext:100000 |
|---|
| 105 | |
|---|
| 106 | arg:NCHAINS |
|---|
| 107 | argtype:text |
|---|
| 108 | arglabel:Number of chains |
|---|
| 109 | argtext:4 |
|---|
| 110 | |
|---|
| 111 | arg:TEMP |
|---|
| 112 | argtype:text |
|---|
| 113 | arglabel:Temperature parameter for heating the chains |
|---|
| 114 | argtext:0.5 |
|---|
| 115 | |
|---|
| 116 | arg:SAMPLEFREQ |
|---|
| 117 | argtype:text |
|---|
| 118 | arglabel:Markov chain sample frequency |
|---|
| 119 | argtext:500 |
|---|
| 120 | |
|---|
| 121 | arg:BURNINFRAC |
|---|
| 122 | argtype:text |
|---|
| 123 | arglabel:Fraction of samples that will be discarded\nwhen convergence diagnostics are calculated [0..1] |
|---|
| 124 | argtext:0.25 |
|---|
| 125 | |
|---|
| 126 | |
|---|
| 127 | in:TmpInputFile |
|---|
| 128 | informat:genbank |
|---|
| 129 | insave: |
|---|
| 130 | |
|---|
| 131 | //------------------------------------------------------------------------------------------ |
|---|
| 132 | |
|---|
| 133 | |
|---|
| 134 | item:MrBayes (custom) |
|---|
| 135 | itemmeta:c |
|---|
| 136 | seqtype:* |
|---|
| 137 | itemmethod:(GOTO_LOCAL_DIR; \ |
|---|
| 138 | GEN_NEXUS(mrbayes_in) && \ |
|---|
| 139 | PREPARE_MRBAYES_INFILE(mrbayes_in) && \ |
|---|
| 140 | ( \ |
|---|
| 141 | `# append the mrbayes command block to the nexus file`; \ |
|---|
| 142 | cat $COMMAND_FILE >> mrbayes_in; \ |
|---|
| 143 | ) && \ |
|---|
| 144 | RUN_IN_WINDOW(( ($ARBHOME/bin/mb mrbayes_in || STOP_ON_FAILURE(MrBayes)) && \ |
|---|
| 145 | NEXUS2NEWICK(mrbayes_out.con.tre,mr_bayes_out.newick) && \ |
|---|
| 146 | $DISPLAY_FUNC)))& |
|---|
| 147 | |
|---|
| 148 | itemhelp:MrBayesCustom.help |
|---|
| 149 | |
|---|
| 150 | arg:DISPLAY_FUNC |
|---|
| 151 | argtype:choice_menu |
|---|
| 152 | arglabel:What to do with the tree? |
|---|
| 153 | argchoice:Load into ARB ('tree_mrBayes_//'):$ARBHOME/bin/arb_read_tree tree_mrbayes_$$ mr_bayes_out.newick "PRG=MrBayes(custom) FILTER=$FILTER RATE=$RATES" |
|---|
| 154 | argchoice:Treetool:treetool mr_bayes_out.newick |
|---|
| 155 | |
|---|
| 156 | arg:COMMAND_FILE |
|---|
| 157 | argtype:filename |
|---|
| 158 | arglabel:MrBayes command file |
|---|
| 159 | |
|---|
| 160 | in:TmpInputFile |
|---|
| 161 | informat:genbank |
|---|
| 162 | insave: |
|---|
| 163 | |
|---|
| 164 | //------------------------------------------------------------------------------------------ |
|---|
| 165 | |
|---|
| 166 | item:TREEPUZZLE ... |
|---|
| 167 | itemmeta: T |
|---|
| 168 | seqtype:* |
|---|
| 169 | itemmethod:(GOTO_LOCAL_DIR; GEN_PHYLIP(infile) && \ |
|---|
| 170 | RUN_IN_WINDOW(((puzzle || STOP_ON_FAILURE(puzzle)) && $DISPLAY_FUNC));\ |
|---|
| 171 | RM_LOCAL_FILES(TmpInputFile outdist) )& |
|---|
| 172 | itemhelp:treepuzzle.help |
|---|
| 173 | |
|---|
| 174 | arg:DISPLAY_FUNC |
|---|
| 175 | argtype:choice_menu |
|---|
| 176 | arglabel:What to do with the tree? |
|---|
| 177 | argchoice:ARB ('tree_puz_//'):arb_read_tree tree_puz_$$ outtree "PRG=PUZZLE FILTER=$FILTER PKG=strimmer/v. haeseler" |
|---|
| 178 | argchoice:Treetool:treetool outtree |
|---|
| 179 | |
|---|
| 180 | in:TmpInputFile |
|---|
| 181 | informat:genbank |
|---|
| 182 | insave: |
|---|
| 183 | |
|---|