seed =  -1

       seqfile = ../../t_lon_grp_127t_part1.phy
      Imapfile = ../../t_lon_grp_127t.imap
*constraintfile = ../../constraints.txt
       jobname = out

  speciesdelimitation = 0 * fixed species tree
* speciesdelimitation = 1 0 2    * species delimitation rjMCMC algorithm0 and finetune(e)
* speciesdelimitation = 1 1 2 1 * species delimitation rjMCMC algorithm1 finetune (a m)
*         speciestree = 1  0.1 0.3 0.123 0.1234  * speciestree pSNL pshrink plambda+ plambda-
	  speciestree = 0
*   speciesmodelprior = 1  * 0: uniform LH; 1:uniform rooted trees; 2: uniformSLH; 3: uniformSRooted

  species&tree = 11  B  A  V  Y  C  D  L  M  P  R  T
                     9  10 11 9  35 5  32 4  3  1  8
(R, ((A,V)q,((((((d[&phi=0.500000,tau-parent=no],T)c)b[&phi=0.500000,tau-parent=yes],(C,Y)r)s,e[&phi=0.500000])f,g[&phi=0.500000])h,((b[&phi=0.500000,tau-parent=no],(((M,(P,D)n)u,f[&phi=0.500000])e,((B,L)v)d[&phi=0.500000,tau-parent=yes])w)a,h[&phi=0.500000])g)x)y)z;
         phase =  0  0  0  0  0  0  0  0  0  0  0

       usedata = 1  * 0: no data (prior); 1:seq like
         nloci = 502  * number of loci in seqfile

     cleandata = 0    * remove sites with ambiguity data (1:yes, 0:no)?

    thetaprior = 3 0.02 e  # invgamma(a, b) for theta
      tauprior = 3 0.01  # invgamma(a, b) for root tau & Dirichlet(a) for other tau's
      phiprior = 1 1
*     heredity = 1 4 4
*    locusrate = 1 5

      finetune = 1

         print = 1 0 0 0   * MCMC samples, locusrate, heredityscalars, Genetrees
        burnin = 200000
      sampfreq = 8
       nsample = 800000
*       threads = 16 17 1
    checkpoint = 100000 100000