code: def make_bcbiornaseq_object(data):
      
      if "bcbiornaseq" not in dd.get_tools_on(data):
          return data
      upload_dir = tz.get_in(("upload", "dir"), data)
      report_dir = os.path.join(upload_dir, "bcbioRNASeq")
      safe_makedir(report_dir)
      organism = dd.get_bcbiornaseq(data).get("organism", None)
      groups = dd.get_bcbiornaseq(data).get("interesting_groups", None)
      loadstring = create_load_string(upload_dir, groups, organism)
      r_file = os.path.join(report_dir, "load_bcbioRNAseq.R")
      with file_transaction(r_file) as tmp_file:
          memoize_write_file(loadstring, tmp_file)
      rcmd = Rscript_cmd()
      with chdir(report_dir):
          do.run([rcmd, "--no-environ", r_file], "Loading bcbioRNASeq object.")
      make_quality_report(data)
      return data
docs: load the initial bcb.rda object using bcbioRNASeq