docs: def _load_sequences_to_reference_gene(self, g_id, force_rerun=False):
          """Load orthologous strain sequences to reference Protein object, save as new pickle"""
          protein_seqs_pickle_path = op.join(self.sequences_by_gene_dir, '{}_protein_withseqs.pckl'.format(g_id))
  
          if ssbio.utils.force_rerun(flag=force_rerun, outfile=protein_seqs_pickle_path):
              protein_pickle_path = self.gene_protein_pickles[g_id]
              protein_pickle = ssbio.io.load_pickle(protein_pickle_path)
  
              for strain, info in self.strain_infodict.items():
                  strain_sequences = SeqIO.index(info['genome_path'], 'fasta')
                  strain_gene_functional = info['functional_genes'][g_id]
                  if strain_gene_functional:
                      # Pull the gene ID of the strain from the orthology matrix
                      strain_gene_key = self.df_orthology_matrix.at[g_id, strain]
                      new_id = '{}_{}'.format(g_id, strain)
                      if protein_pickle.sequences.has_id(new_id):
                          continue
                      protein_pickle.load_manual_sequence(seq=strain_sequences[strain_gene_key],
                                                          ident=new_id,
                                                          set_as_representative=False)
              protein_pickle.save_pickle(outfile=protein_seqs_pickle_path)
  
          return g_id, protein_seqs_pickle_path