Let $\mu$ denote the mean of the process, then the first order autocovariance is given by:
$$
E\left[(X_t - \mu)(X_{t-1} - \mu)\right] = E\left[\tilde{X}_t \tilde{X}_{t-1}\right] = \\
E\left[(\phi \tilde{X}_{t-1} + \epsilon_t) \tilde{X}_{t-1}\right] =
\phi \underbrace{E\left[\tilde{X}_{t-1}^2\right]}_{\sigma^2_X} +
\underbrace{E\left[\epsilon_t \tilde{X}_{t-1}\right]}_{0} = \phi\sigma^2_X \,.
\qquad (1)
$$
The second expectation is zero because $\tilde{X}_{t-1}$ dependes on the innovations $\epsilon_{t-1}$, $\epsilon_{t-2}$,... but not on $\epsilon_t$.
The first expectation is the variance of the AR(1) process, denoted $\sigma^2_X$. The expression for the variance of $X_t$ can be obtained as:
$$
\hbox{Var }(X_t) = \phi^2 \hbox{Var }(X_{t-1}) +
\underbrace{\hbox{Var }(\epsilon_t)}_{\sigma^2_\epsilon} \,.
$$
As the process is stationary $\hbox{Var }(X_t) = \hbox{Var }(X_{t-1})$. Thus, the above expression can be written as:
$$
\hbox{Var }(X_t) = \frac{\sigma^2_\epsilon}{1 - \phi^2} \,. \qquad (2)
$$
Substituting (2) in (1) gives the following expression for the first order autocovariance:
$$
\frac{\phi \sigma^2_\epsilon}{1 - \phi^2} \,.
$$