samtools view -h sfm-2-best-alignment.bam  | pmdtools --threshold {threshold} --writesamfield --header | '
        f'samtools view -Sb - > "{output_path}"'