Initially assumes all reads to contain only SNV's. Also added new argument to program max_variants which will ignore reads with more than that many variants observed.