  We study a 3-3-1 model based on non-Abelian discrete symmetry group $T_{13}$
which accommodates lepton mixing with non-zero $\theta_{13}$ and CP violation
phase. The neutrinos get small masses and mixing with CP violation phase from
$SU(3)_L$ antisextets which are all in triplets under $T_{13}$. If both
breakings $T_{13}\rightarrow Z_3$ and $Z_{3}\rightarrow \{\mathrm{Identity}\}$
are taken palce in neutrino sector and $T_{13}$ is broken into $Z_3$ in lepton
sector, the realistic neutrino mixing form is obtained as a natural consequence
of $P_l$ and $T_{13}$ symmetries. The model predicts the lepton mixing with
non-zero $\theta_{13}$, and also gives a remarkable prediction of Dirac CP
violation $\delta_{CP}=292.5^\circ$ in the normal spectrum, and
$\delta_{CP}=303.161^\circ$ in the inverted spectrum which is still missing in
the neutrino mixing matrix. There exist some regions of model parameters that
can fit the experimental data in 2014 on neutrino masses and mixing without
perturbation.
