estimate_m_rd function

Estimate heterozygote miscall rate for different read depth categories (no nulls)