The Fermi bubbles are two large structures in the gamma-ray sky extending to 55degrees above and below the Galactic center. We analyze 50 months of Fermi Large Area Telescope data between 100 MeV and 500 GeV above 10degrees in Galactic latitude to derive the spectrum and morphology of the Fermi bubbles. We thoroughly explore the systematic uncertainties that arise when modeling the Galactic diffuse emission through two separate approaches. The gamma-ray spectrum is well described by either a log parabola or a power law with an exponential cutoff. We exclude a simple power law with more than 7sigma significance. The power law with an exponential cutoff has an index of 1.9 + or - 0.2 and a cutoff energy of 110 + or - 50 GeV. We find that the gamma-ray luminosity of the bubbles is (ProQuest: Formulae and/or non-USASCII text omitted)erg s super(-1). We confirm a significant enhancement of gamma-ray emission in the southeastern part of the bubbles, but we do not find significant evidence for a jet. No significant variation of the spectrum across the bubbles is detected. The width of the boundary of the bubbles is estimated to be (ProQuest: Formulae and/or non-USASCII text omitted) deg. Both inverse Compton (IC) models and hadronic models including IC emission from secondary leptons fit the gamma-ray data well. In the IC scenario, synchrotron emission from the same population of electrons can also explain the WMAP and Planck microwave haze with a magnetic field between 5 and 20 mu G.