We consider the dynamics in infinite population evolution models with a general symmetric fitness landscape. We find shock waves, i.e., discontinuous transitions in the mean fitness, in evolution dynamics even with smooth fitness landscapes, which means that the search for the optimal evolution trajectory is more complicated. These shock waves appear in the case of positive epistasis and can be used to represent punctuated equilibria in biological evolution during long geological time scales. We find exact analytical solutions for discontinuous dynamics at the large-genome-length limit and derive optimal mutation rates for a fixed fitness landscape to send the population from the initial configuration to some final configuration in the fastest way.