code: def _add_segmetrics_to_output(out, data):
      
      out_file = "%s-segmetrics.txt" % os.path.splitext(out["cns"])[0]
      if not utils.file_exists(out_file):
          with file_transaction(data, out_file) as tx_out_file:
              cmd = [os.path.join(os.path.dirname(sys.executable), "cnvkit.py"), "segmetrics",
                     "--median", "--iqr", "--ci", "--pi",
                     "-s", out["cns"], "-o", tx_out_file, out["cnr"]]
              # Use less fine grained bootstrapping intervals for whole genome runs
              if dd.get_coverage_interval(data) == "genome":
                  cmd += ["--alpha", "0.1", "--bootstrap", "50"]
              else:
                  cmd += ["--alpha", "0.01", "--bootstrap", "500"]
              do.run(cmd, "CNVkit segmetrics")
      out["segmetrics"] = out_file
      return out
docs: Add metrics for measuring reliability of CNV estimates.