#!/bin/bash -euo pipefail
set -euxo pipefail

seqtool saturation -region 'chr20' \
    "rna_s1221_S42.md.bam" > rna_s1221_S42_full_bam_saturation_values.txt