Turbulent mixing layers are canonical flow in nature and engineering, and deserve comprehensive studies under various conditions using different methods. In this paper, turbulent mixing layers are investigated using large eddy simulation and dynamic mode decomposition. The accuracy of the computations is verified and validated. Standard dynamic mode decomposition is utilized to flow decomposition, reconstruction and prediction. It was found that the dominant-mode selection criterion based on mode amplitude is more suitable for turbulent mixing layer flow compared with the other three criteria based on singular value, modal energy and integral modal amplitude, respectively. For the mixing layer with random disturbance, the standard dynamic mode decomposition method could accurately reconstruct and predict the region before instability happens, but is not qualified in the regions after that, which implies that improved dynamic mode decomposition methods need to be utilized or developed for the future dynamic mode decomposition of turbulent mixing layers.