We present an analysis of parameter variance statistics for the TLS-Prony method applied to damped exponential signals. We derive the covariance matrix of the estimated parameters for this method. The parameters include the magnitudes and angles of the poles, and the magnitudes and angles of the amplitude coefficients. We verify the theoretical results using Monte-Carlo simulations studies. We also compare the variance results to the corresponding Cramer-Rao bounds for several cases.