code: def _report(data, reference):
      
      seqcluster = op.join(get_bcbio_bin(), "seqcluster")
      work_dir = dd.get_work_dir(data)
      out_dir = safe_makedir(os.path.join(work_dir, "seqcluster", "report"))
      out_file = op.join(out_dir, "seqcluster.db")
      json = op.join(work_dir, "seqcluster", "cluster", "seqcluster.json")
      cmd = ("{seqcluster} report -o {out_dir} -r {reference} -j {json}")
      if not file_exists(out_file):
          do.run(cmd.format(**locals()), "Run report on clusters")
      return out_file
docs: Run report of seqcluster to get browser options for results