def index_bam(in_bam): if not os.path.isfile(in_bam + ".bai") or os.stat(in_bam).st_mtime >= os.stat(in_bam + ".bai").st_mtime: pysam.index(in_bam, catch_stdout=False)