/////////////////////////////////////////////////////////////////////////////////////
//                                                                                 //
//  INDELible V1.03 control file - basicnucleotide.txt                             //
//                                                                                 //
//      A basic introduction to the structure of the INDELible control file.       //
//                                                                                 //
/////////////////////////////////////////////////////////////////////////////////////

// It is useful to know that anything on a line after two forward slashes is ignored.

/*
   Another useful thing to know is that anything after a forward slash and star
   is ignored until INDELible sees a star followed by a forward slash later on.
*/     

[TYPE] NUCLEOTIDE 1	//  EVERY control file must begin with a [TYPE] command.
			//  The number after "NUCLEOTIDE" can be 1 or 2 and chooses the 
			//  algorithm that INDELible uses (see manuscript). Both give 
			//  identical results but in some cases one is quicker.
			//  Other blocks and commands following this statement
			//  can come in any order you like.

[MODEL]    modelname        //  Evolutionary models are defined in [MODEL] blocks.
  [submodel]     JC         //  Here the substitution model is simply set as JC69.
  [indelmodel]   NB  0.4 1  //  Geometric indel length distribution (q=0.4, r=1)
  [insertrate]   0.08       //  insertion rate = 0.08 relative to substitution rate of 1
  [deleterate]   0.12       //  deletion rate = 0.12 relative to substitution rate of 1
  
[TREE] treename  (A:0.1,B:0.1);        //  User trees are defined here

[PARTITIONS] partitionname             //  [PARTITIONS] blocks say which models go with
  [treename modelname 1000]            //  which trees and define the length of the
                                       //  sequence generated at the root (1000 here).

[EVOLVE] partitionname 100 outputname  //  This will generate 100 replicate datasets 
                                       //  from the [PARTITIONS] block named above.

// The true alignment will be output in a file named outputname_TRUE.phy
// The unaligned sequences will be output in a file named outputname.fas
// To learn how to implement more complicated simulations (or different 
// models) please consult the manual or the other example control files.

---

Please click here to return to the tutorial menu page for more examples.