wrote e-step function for theta, mstep function for sequencing error rate