Subcritical transition to turbulence requires finite-amplitude perturbations. Using a nonlinear optimisation technique in a periodic computational domain, we identify the perturbations of plane Couette flow transitioning with least initial kinetic energy for Re <= 3000. We suggest a new scaling law E-c = O(Re-2.7) for the energy threshold vs. the Reynolds number, in quantitative agreement with experimental estimates for pipe flow. The route to turbulence associated with such spatially localised perturbations is analysed in detail for Re = 1500. Several known mechanisms are found to occur one after the other: Orr mechanism, oblique wave interaction, lift-up, streak bending, streak breakdown, and spanwise spreading. The phenomenon of streak breakdown is analysed in terms of leading finite-time Lyapunov exponents of the associated edge trajectory.