Take insertions/deletions into account and script to do systematic testing with the number of reads