INCLUDE dna.reform

# Generate a Root consisting of only one segment:
# The segment receives a random uniform sequence of length 100.
# The default Model (unless "overwritten" later) is JC.
ROOT r DNA
{
  SEGMENT TheOnlySegment
    { SEQUENCE  1000 * uniform   MODEL JC }
}


# Have a "human" and a "chimp" sequence evolve from the root
NODE human { PARENT r  TIME 25 }
NODE chimp { PARENT r  TIME 35 }

# Output the leaf sequences
OUTPUT FASTA { human chimp }