We present a coherent timing analysis of the 2003 outburst of the accreting millisecond pulsar (AMXP) XTEJ1807-294. We find a 95% confidence interval for the pulse frequency derivative of (+0.7, +4.7) x 10(-14) Hz s(-1) and (-0.6, +3.8) x 10(-14) Hz s(-1) for the fundamental and second harmonics, respectively. The sinusoidal fractional amplitudes of the pulsations are the highest observed among AMXPs and can reach values of up to 27% (2.5-30 keV). The pulse arrival time residuals of the fundamental frequency follow a linear anti-correlation with the fractional amplitudes that suggests hot spot motion both in longitude and latitude over the surface of the neutron star. An anti-correlation between residuals and X-ray flux suggests an influenc...