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)