We demonstrate that the nonlinear optimization of a finite-amplitude disturbance over a freely evolving and possibly even turbulent flow, can successfully identify subcritical dynamo branches as well as the structure and amplitude of their critical perturbations. As this approach does not require prior knowledge of the magnetic field amplification mechanisms, it opens a new avenue for systematically probing subcritical dynamo flows.