sim.seqs function

Simulate DNA sequences according to DNA substitution models