docs: def store_disorder(self, sc=None, force_rerun=False):
          """Wrapper for _store_disorder"""
          log.info('Loading sequences to reference GEM-PRO...')
          from random import shuffle
          g_ids = [g.id for g in self.reference_gempro.functional_genes]
          shuffle(g_ids)
  
          def _store_disorder_sc(g_id, outdir=self.sequences_by_gene_dir,
                                 g_to_pickle=self.gene_protein_pickles, force_rerun=force_rerun):
              """Load orthologous strain sequences to reference Protein object, save as new pickle"""
              import ssbio.utils
              import ssbio.io
              import os.path as op
  
              protein_seqs_pickle_path = op.join(outdir, '{}_protein_withseqs_dis.pckl'.format(g_id))
  
              if ssbio.utils.force_rerun(flag=force_rerun, outfile=protein_seqs_pickle_path):
                  protein_pickle_path = g_to_pickle[g_id]
                  protein_pickle = ssbio.io.load_pickle(protein_pickle_path)
                  protein_pickle.get_all_disorder_predictions(representative_only=False)
                  protein_pickle.save_pickle(outfile=protein_seqs_pickle_path)
  
              return g_id, protein_seqs_pickle_path
  
          if sc:
              genes_rdd = sc.parallelize(g_ids)
              result = genes_rdd.map(_store_disorder_sc).collect()
          else:
              result = []
              for g in tqdm(g_ids):
                  result.append(self._load_sequences_to_reference_gene(g, force_rerun))
  
          log.info('Storing paths to new Protein objects in self.gene_protein_pickles...')
          for g_id, protein_pickle in result:
              self.gene_protein_pickles[g_id] = protein_pickle