docs: def simulate(self, nsites, transition_matrix, tree, ncat=1, alpha=1):
          """
          Return sequences simulated under the transition matrix's model 
          """
          sim = SequenceSimulator(transition_matrix, tree, ncat, alpha)
          return list(sim.simulate(nsites).items())