let somatic_variant_caller t bam_pair =
    Somatic_variant_caller (t, bam_pair)