def index(file_name): pysam.index(str(file_name)) if not os.path.isfile(file_name + ".bai"): raise RuntimeError("samtools index failed, likely bam is not sorted")