Quantum statistics of multimode parametric amplification