Rate constants and the product branching ratio for hydrogen abstraction from CH3OH by O(P-3) were computed with multistructural variational transition-state theory including microcanonically optimized multidimensional tunneling. Benchmark calculations of the forward and reverse classical barrier heights and the reaction energetics have been carried out by using coupled cluster theory and multireference calculations to select the most reliable density functional method for direct dynamics computations of the rate constants. The dynamics calculations included the anharmonicity of the zero-point energies and partition functions, with specific-reaction-parameter scaling factors for reactants and transition states, and multistructural torsional anharmonicity was included for the torsion around the C-O bond in methanol and in the transition states. The resulting rate constants are presented over a wider range than they are available from experiment, but in the temperature range where experiments are available, they agree well with experimental values, which is encouraging for their reliability over the wider temperature range and for future computations of oxygen atom reaction rates. In contrast to a previous computational prediction, the branching ratio predicted by the present work shows that the formation of CH2OH + OH is the dominant channel over the whole range of temperature from 250 to 2000 K.