docs: def load_sequences_to_strains(self, joblib=False, cores=1, force_rerun=False):
          """Wrapper function for _load_sequences_to_strain"""
          log.info('Loading sequences to strain GEM-PROs...')
          if joblib:
              result = DictList(Parallel(n_jobs=cores)(delayed(self._load_sequences_to_strain)(s, force_rerun) for s in self.strain_ids))
          else:
              result = []
              for s in tqdm(self.strain_ids):
                  result.append(self._load_sequences_to_strain(s, force_rerun))
  
          for strain_id, gp_seqs_path in result:
              self.strain_infodict[strain_id]['gp_seqs_path'] = gp_seqs_path