code: def read_hlas(fasta_fai):
      
      out = []
      with open(fasta_fai) as in_handle:
          for line in in_handle:
              if line.startswith("HLA"):
                  out.append(line.split()[0])
      return out
docs: Get HLA alleles from the hg38 fasta fai file.