script uses pysam now