code: def get_hla(sample, cromwell_dir, hla_glob):
      
      hla_dir = glob.glob(os.path.join(cromwell_dir, hla_glob, "align", sample, "hla"))[0]
      fastq = os.path.join(hla_dir, "OptiType-HLA-A_B_C-input.fq")
      calls = os.path.join(hla_dir, "%s-optitype.csv" % sample)
      return fastq, calls
docs: Retrieve HLA calls and input fastqs for a sample.