High fidelity polymerases are efficient catalysts of phosphodiester bond formation during DNA replication or repair. We interpret molecular dynamics simulations of a polymerase bound to its substrate DNA and incoming nucleotide using a quasiharmonic model to study the effect of external forces applied to the bound DNA on the kinetics of phosphoryl transfer. The origin of the force dependence is shown to be an intriguing coupling between slow, delocalized polymerase-DNA modes and fast catalytic site motions. Using noncognate DNA substrates we show that the force dependence is context specific.