code: def print_mhc_peptide(neoepitope_info, peptides, pepmap, outfile):
      
      allele, pept, pred, core = neoepitope_info
      peptide_names = [x for x, y in peptides.items() if pept in y]
      # For each peptide, append the ensembl gene
      for peptide_name in peptide_names:
          print(allele, pept, peptide_name, core, '0', pred, pepmap[peptide_name], sep='\t',
                file=outfile)
      return None
docs: To reduce code redundancy, this module will accept data from merge_mhc_peptide_calls for a given
      neoepitope and print it to outfile
      ARGUMENTS
      1. neoepitope_info: Tuple of (<allele>, <peptide_sequence>,
                                    <binding_prediction>)
      2. peptides: Dict of all IARS considered
             peptides
                |- 'neoepitope_1': <peptide_sequence>
                ..
                |- 'neoepitope_n': <peptide_sequence>
      3. pepmap: Info correlating neoepitope with the gene and transcript level
                 mutations.
             peptides
                |- 'neoepitope_1':
                |      'ensembl_gene\thugo_gene\tcomma_sep_transcript_mutations'
                ..
                +- 'neoepitope_n':
                       'ensembl_gene\thugo_gene\tcomma_sep_transcript_mutations'