Output 100 pairs per position

Changed subsample_reads method to output 100 pairs for each random position. This fixes a bias in insert size calculation. In addition, added a filter to make sure reads are a proper pair, not pcr dups, and not a qcfail.