The goal of this paper is to investigate the effect that a distribution of kinesin motor velocities could have on cytoskeletal element (CE) concentration waves in slow axonal transport. Previous models of slow axonal transport based on the stop-and-go hypothesis (P. Jung, A. Brown, Modeling the slowing of neurofilament transport along the mouse sciatic nerve, Physical Biology 6 (2009) 046002) assumed that in the anterograde running state all CEs move with one and the same velocity as they are propelled by kinesin motors. This paper extends the aforementioned theoretical approach by allowing for a distribution of kinesin motor velocities; the distribution is described by a probability density function (PDF). For a two kinetic state model (that accounts for the pausing and running populations of CEs) an analytical solution describing the propagation of the CE concentration wave is derived. Published experimental data are used to obtain an analytical expression for the PDF characterizing the kinesin velocity distribution; this analytical expression is then utilized as an input for computations. It is demonstrated that accounting for the kinesin velocity distribution increases the rate of spreading of the CE concentration waves, which is a significant improvement in the two kinetic state model.