Updated, now using pysam