We present a new theoretical method to determine and visualize the average tunneling route of the electron transfer (ET) in protein media. In this, we properly took into account the fluctuation of the tunneling currents and the quantum-interference effect. The route was correlated with the electronic factor < T-DA(2)> in the case of ET by the elastic tunneling mechanism. We expanded < T-DA(2)> by the interatomic tunneling currents J(ab)'s. Incorporating the quantum-interference effect into the mean-square interatomic tunneling currents, denoted as <(J) over tilde (2)(ab)>, we could express < T-DA(2)> as a sum of h(2)<(J) over tilde (2)(ab)>. Drawing the distribution of <(J) over tilde (2)(ab)> on the protein structure, we obtain the <(J) over tilde (2)(ab)> map which visually represents which parts of bonds and spaces most significantly contribute to < T-DA(2)>. We applied this method to the ET from the bacteriopheophytin anion to the primary quinone in the bacterial photosynthetic reaction center of Rhodobacter sphaeroides. We obtained J(ab)'s by a combined method of molecular dynamics simulations and quantum chemical calculations. In calculating <(J) over tilde (2)(ab)>, we found that much destructive interference works among the interatomic tunneling currents even after taking the average. We drew the <(J) over tilde (2)(ab)> map by a pipe model where atoms a and b are connected by a pipe with width proportional to the magnitude Of <(J) over tilde (2)(ab)>. We found that two groups of <(J) over tilde (2)(ab)>'s which are mutually coupled with high correlation in each group, have broad pipes and form the average tunneling routes, called Trp route and Met route. Each of the two average tunneling routes is composed of a few major pathways in the Pathways model which are fused at considerable part to each other. We also analyzed the average tunneling route for the ET by the inelastic tunneling mechanism.