A directional spherical wavelet analysis is performed to examine the Gaussianity of the WMAP 1-year data. Such an analysis is facilitated by the introduction of a fast directional continuous spherical wavelet transform. Significant deviations from Gaussianity are detected in the skewness and kurtosis of spherical elliptical Mexican hat and real Morlet wavelet coefficients for both the WMAP and Tegmark foreground removed maps. In particular, the skewness of spherical real Morlet wavelet coefficients on wavelet scale 550' (corresponding to an effective size on the sky of approximately 26 degrees) and at an azimuthal orientation of 72 degrees, deviates by 5.88 standard deviations for the WMAP map and by 6.72 standard deviations for the Tegmark map. These detections are made at the 99% significance level. Furthermore, the wavelet analysis allows one to localise on the sky the most pronounced deviations from Gaussianity. When these regions are removed from the analysis the data are consistent with Gaussianity, suggesting these localised regions are the source of the detected non-Gaussianity.