Technical/Simulation/Codon/MCcodonNSbranchsites.dat

   1 0          * 0,1:seqs or patters in paml format (mc.paml); 2:paup format (mc.nex)\r
   2 113147       * random number seed (odd number)\r
   3 8  10000  5  * <# seqs>  <# codons>  <# replicates>\r
   4 \r
   5 -1           * <tree length, use -1 if tree has absolute branch lengths>\r
   6 \r
   7 (((a:0.1, b:0.2):0.12, (c:0.3, d:0.4):0.34):0.1234, ((e:0.5, f:0.6):0.56, (g:0.7, h:0.8):0.78):0.5678);\r
   8 \r
   9 4                 * number of site classes, & frequencies for site classes\r
  10   0.2  0.2  0.3  0.3  * must sum to 1.\r
  11 \r
  12 (((a #0.1, b #0.1) #0.1, (c #0.1, d #0.1) #0.1) #0.1, ((e #0.1, f #0.1) #0.1, (g #0.1, h #0.1) #0.1) #0.1);\r
  13 (((a #1.0, b #1.0) #1.0, (c #1.0, d #1.0) #1.0) #1.0, ((e #1.0, f #1.0) #1.0, (g #1.0, h #1.0) #1.0) #1.0);\r
  14 (((a #0.1, b #0.1) #3.0, (c #0.1, d #0.1) #0.1) #0.1, ((e #0.1, f #0.1) #3.0, (g #0.1, h #0.1) #0.1) #0.1);\r
  15 (((a #1.0, b #1.0) #3.0, (c #1.0, d #1.0) #1.0) #1.0, ((e #1.0, f #1.0) #3.0, (g #1.0, h #1.0) #1.0) #1.0);\r
  16 \r
  17 2.5    * kappa\r
  18 \r
  19 0.01639344 0.01639344 0.01639344 0.01639344 \r
  20 0.01639344 0.01639344 0.01639344 0.01639344 \r
  21 0.01639344 0.01639344 0          0\r
  22 0.01639344 0.01639344 0          0.01639344 \r
  23 \r
  24 0.01639344 0.01639344 0.01639344 0.01639344 \r
  25 0.01639344 0.01639344 0.01639344 0.01639344 \r
  26 0.01639344 0.01639344 0.01639344 0.01639344 \r
  27 0.01639344 0.01639344 0.01639344 0.01639344 \r
  28 \r
  29 0.01639344 0.01639344 0.01639344 0.01639344 \r
  30 0.01639344 0.01639344 0.01639344 0.01639344 \r
  31 0.01639344 0.01639344 0.01639344 0.01639344 \r
  32 0.01639344 0.01639344 0.01639344 0.01639344 \r
  33 \r
  34 0.01639344 0.01639344 0.01639344 0.01639344 \r
  35 0.01639344 0.01639344 0.01639344 0.01639344 \r
  36 0.01639344 0.01639344 0.01639344 0.01639344 \r
  37 0.01639344 0.01639344 0.01639344 0.01639344 \r
  38 \r
  39 0    * genetic code (0:universal; 1:mammalian mt; 2-10:see below)\r
  40 \r
  41 \r
  42 //end of file\r
  43 \r
  44 ============================================================================\r
  45 MCcodonNSbranches.dat\r
  46 Notes for using evolver to simulate codon sequences under the branch model, \r
  47 with different omega's for different branches. \r
  48 \r
  49 Change values of parameters, but do not delete them.  You can add\r
  50 empty lines, but do not break one line into several lines or merge several \r
  51 lines into one line.\r
  52 \r
  53 Use : to indicate the branch length in the tree topology.  Use # to\r
  54 indicate the omega ratio for the branch.\r
  55 \r
  56 64 codon freqs are in fixed order TTT, TTC, TTA, TTG, TCT, TCC, ..., GGG.\r
  57 Stop codons have frequency 0.\r
  58 \r
  59 =================!! Check screen output carefully!! ========================\r
  60 \r
  61 nice trees:\r
  62 \r
  63 (((a:0.2, b:0.2):0.2, (c:0.2, d:0.2):0.2):0.2, ((e:0.2, f:0.2):0.2, (g:0.2, h:0.2):0.2):0.2);\r