The problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem to find the initial condition function (analysis). The data contain errors (observation and background errors), hence there is an error in the analysis. For mildly nonlinear dynamics, the analysis error covariance can be approximated by the inverse Hessian of the cost functional in the auxiliary data assimilation problem, whereas for stronger nonlinearity - by the 'effective' inverse Hessian. However, it has been noticed that the analysis error covariance is not the posterior covariance from the Bayesian perspective. While these two are equivalent in the linear case, the difference may become significant in practical terms with the nonlinearity level rising. For the proper Bayesian posterior covariance a new approximation via the Hessian is derived and its 'effective' counterpart is introduced. An approach for computing the mentioned estimates in the matrix-free environment using Lanczos method with preconditioning is suggested. Numerical examples which validate the developed theory are presented for the model governed by Burgers equation with a nonlinear viscous term.