********************************************************************************
INDELible V1.03 by Will Fletcher : Simulation began at: Fri Jun 12 01:59:44 2020
********************************************************************************

  * Block 1 was completed in 0.000462 seconds.

********************************************************************************


  * Block 2 was completed in 0.00052 seconds.

********************************************************************************


  * Block 3 was completed in 0.000428 seconds.

********************************************************************************


  * Block 4 was completed in 0.000395 seconds.

********************************************************************************


  * Block 5 was completed in 0.000303 seconds.

********************************************************************************

  * Simulation completed. Whole batch took: 0.003156 seconds.

INDELible V1.03 Simulations completed at: Fri Jun 12 01:59:44 2020

********************************************************************************


 Original Control File 

-----------------------


/////////////////////////////////////////////////////////////////////////////////////
//                                                                                 //
//  INDELible V1.03 control file - rates.txt                                       //
//                                                                                 //
//      Adding proportion of invariable sites or gamma rate heterogeneity.         //
//                                                                                 //
/////////////////////////////////////////////////////////////////////////////////////

/*
  Nucleotide models are used as an example but the principle is exactly the 
  same for amino-acid models. You just add the [rates] command to a model:

    [rates]  pinv  alpha  ngamcat

  This command has no effect in CODON simulations.

  * pinv is the proportion of invariable sites (0<=pinv<=1).
  * alpha is the shape parameter for the gamma distribution (alpha>0).
  * If alpha=0 then there will be no gamma rate variation.
  * ngamcat is the number of categories to use in the discrete gamma approximation.
  * If ngamcat=0 then continuous gamma distribution will be used for rate variation.
  * If the command [rates] is not specified all rates will be constant.

*/


[TYPE] NUCLEOTIDE 1	//  nucleotide simulation using algorithm from method 1

[SETTINGS]
  [randomseed] 1568746 // seed for random generator so outputs can be checked

[MODEL] mymodel1  [submodel]  JC         //   JC69 with pinv=0.2 and no gamma
                  [rates]     0.2 0 0       
 									
[MODEL] mymodel2  [submodel]  JC         //   JC69 with no proportion invariable,
                  [rates]     0 2 0      //   continuous gamma with alpha=2     

[MODEL] mymodel3  [submodel]  JC         //   JC69 with no proportion invariable,
                  [rates]     0 1 10     //   10 category discrete gamma with alpha=1
                                                                    
[MODEL] mymodel4  [submodel]  JC         //   JC69 with pinv=0.2 
                  [rates]     0.2 0.5 0  //   continuous gamma with alpha=0.5     

[MODEL] mymodel5  [submodel]  JC         //   JC69 with pinv=0.1
                  [rates]     0.1 1.5 5  //   5 category discrete gamma with alpha=1.5


[TREE] T1  (A:0.1,B:0.1);             


[PARTITIONS] Pname1  [T1 mymodel1 20]   // tree T1, model mymodel1, root length 20
[PARTITIONS] Pname2  [T1 mymodel2 20]   // tree T1, model mymodel2, root length 20
[PARTITIONS] Pname3  [T1 mymodel3 20]   // tree T1, model mymodel3, root length 20
[PARTITIONS] Pname4  [T1 mymodel4 20]   // tree T1, model mymodel4, root length 20
[PARTITIONS] Pname5  [T1 mymodel5 20]   // tree T1, model mymodel5, root length 20



 
[EVOLVE]     Pname1  5  outputname1     // 5 replicates generated from partition Pname1
             Pname2  5  outputname2     // 5 replicates generated from partition Pname2
             Pname3  5  outputname3     // 5 replicates generated from partition Pname3
             Pname4  5  outputname4     // 5 replicates generated from partition Pname4
             Pname5  5  outputname5     // 5 replicates generated from partition Pname5
