docs: def mutate(self, chromosomes, p_mutate):
          """ 
          Call every chromosome's ``mutate`` method. 
          
          p_mutate:  probability of mutation in [0, 1]
          """
          assert 0 <= p_mutate <= 1
          
          for chromosome in chromosomes:
              chromosome.mutate(p_mutate)