calculate JSD in each simulation; add noise to read number