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 $MLMODEL -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=!$MLMODEL" |
---|
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:MLMODEL |
---|
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:Mixed invariable/gamma: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 | |
---|