nicely outputs results
calculate peak binding in each pass... this doesn't make any sense until EM has converged