mass eigenvalues and $\gamma_{1,2,3}$ are the respective Majorana phases. Therefore the effective neutrino mixing matrix (with diagonal charged lepton sector), can be written as $$\label{unu}
 U_{\nu}=U_{\rm TBM}U_{13}U_{m}$$ where $U_m$ is the diagonal Majorana phase matrix given by $U_m={\rm diag} (e^{i\gamma_1},e^{i\gamma_2},e^{i\gamma_3})$. Comparing this with the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix $$U_{\text{PMNS}}=\left(\begin{array}{ccc}
c_{12}c_{13}& s_{12}c_{13}& s_{13}e^{-i\delta}\\
-s_{12}c_{23}-c_{12}s_{23}s_{13}e^{i\delta}& c_{12}c_{23}-s_{12}s_{23}s_{13}e^{i\delta} & s_{23}c_{13} \\
s_{12}s_{23}-c_{12}c_{23}s_{13}e^{i\delta} & -c_{12}s_{23}-s_{12}c_{23}s_{13}e^{i\delta}& c_{23}c_{13}
\end{array}\right)P,
\label{PMNS}$$ one can obtain the correlation between the model parameters and neutrino oscillation parameters. In the above PMNS mixing matrix, $P (={\rm diag} (1,e^{i\alpha_{21}/2},e^{i\alpha_{31}/2}))$ is the Majorana phase matrix where $\alpha_{21}=(\gamma_1-\gamma_2)$ and $\alpha_{31}=(\gamma_1-\gamma_3)$ are the physical Majorana phases after rotating away one common phase which is irrelevant. The rotation parameter $\theta$ in the rotation matrix can be obtained as $$\begin{aligned}
 \tan 2\theta=\pm \sqrt{3}\frac{x_1}{x_2}, \end{aligned}$$ with $\psi=0,\pi$ respectively yielding maximal value for the Dirac CP phase $\delta$. The construction of the present set-up is such that we only have considered the Planck scale suppressed terms only for the Dirac neutrino mass matrix. Once we allow correction terms for the charged lepton and/or right handed neutrino mass matrix, non-maximal values for the Dirac CP phase $\delta$ becomes allowed [@Karmakar:2014dva; @Karmakar:2015jza; @Karmakar:2016cvb]. Comparing the neutrino mixing matrix of our scenario given 