An asymptotic method for stability analysis of quasilinear functional differential equations, with small perturbations dependent on phase coordinates and an ergodic Markov process, is presented. The proposed method is based on an averaging procedure with respect to: 1) time along critical solutions of the linear equation; and 2) the invariant measure of the Markov process. For asymptotic analysis of the initial random equation with delay, it is proved that one can approximate its solutions (which are stochastic processes) by corresponding solutions of a specially constructed averaged, deterministic ordinary differential equation. Moreover, it is proved that exponential stability of the resulting deterministic equation is sufficient for exponential p-stability of the initial random system for all positive numbers p, and for sufficiently small perturbation terms.