We study the inspiral, merger and ringdown of unequal mass black hole binaries by analyzing a catalogue of numerical simulations for seven different values of the mass ratio (from q = M2/M1 = 1 to q = 4). We compare numerical and Post-Newtonian results by projecting the waveforms onto spin-weighted spherical harmonics, characterized by angular indices (l , m). We find that the PostNewtonian equations predict remarkably well the relation between the wave amplitude and the orbital frequency for each (l , m), and that the convergence of the Post-Newtonian series to the numerical results is non-monotonic. To leading order the total energy emitted in the merger phase scales like η 2 and the spin of the final black hole scales like η, where η = q/(1 + q) 2 is the symmetric mass ratio. We study the multipolar distribution of the radiation, finding that odd-l multipoles are suppressed in the equal mass limit. Higher multipoles carry a larger fraction of the total energy as q increases. We introduce and compare three different definitions for the ringdown starting time. Applying linear estimation methods (the so-called Prony methods) to the ringdown phase, we find resolution-dependent time variations in the fitted parameters of the final black hole. By crosscorrelating information from different multipoles we show that ringdown fits can be used to obtain precise estimates of the mass and spin of the final black hole, which are in remarkable agreement with energy and angular momentum balance calculations.