@sjcrh/proteinpaint-client 2.223.0 → 2.224.1

package/dist/chunk-O3XTTT4K.js

@@ -0,0 +1,217 @@
+import {
+  dofetch,
+  dofetch2
+} from "./chunk-KBTFT76L.js";
+import {
+  contigNameNoChr2
+} from "./chunk-WJLRU5OF.js";
+
+// tracks/hic/data/parseData.ts
+async function hicParseFile(hic, debugmode, errList = []) {
+  if (debugmode) window["hic"] = hic;
+  if (hic.tklst) {
+    const lst = [];
+    for (const t of hic.tklst) {
+      if (!t.type) {
+        errList.push("type missing from one of the tracks accompanying HiC");
+      } else {
+        t.iscustom = true;
+        lst.push(t);
+      }
+    }
+    if (lst.length) {
+      hic.tklst = lst;
+    } else {
+      delete hic.tklst;
+    }
+  }
+  if (hic.enzyme) {
+    if (hic.genome.hicenzymefragment) {
+      let frag = null;
+      for (const f of hic.genome.hicenzymefragment) {
+        if (f.enzyme == hic.enzyme) {
+          frag = f;
+          break;
+        }
+      }
+      if (frag) {
+        hic.enzymefile = frag.file;
+      } else {
+        errList.push("unknown enzyme: " + hic.enzyme);
+        delete hic.enzyme;
+      }
+    } else {
+      errList.push("no enzyme fragment information available for this genome");
+      delete hic.enzyme;
+    }
+  }
+  try {
+    if (hic.sv && hic.sv.file) {
+      const re = await dofetch(hic.hostURL + "/textfile", {
+        method: "POST",
+        body: JSON.stringify({ file: hic.sv.file, jwt: hic.jwt })
+      });
+      const data2 = re.json();
+      const [err2, header, items] = parseSV(data2.text);
+      if (err2) throw { message: "Error parsing SV: " + err2 };
+      hic.sv.header = header;
+      hic.sv.items = items;
+    }
+    const data = await dofetch2("hicstat?" + (hic.file ? "file=" + hic.file : "url=" + hic.url));
+    if (data.error) {
+      errList.push(data.error);
+      return;
+    }
+    const err = hicparsestat(hic, data.out);
+    if (err) throw { message: err };
+  } catch (err) {
+    errList.push(err.message || err);
+    if (err.stack) {
+      console.log(err.stack);
+    }
+  }
+  return hic;
+}
+function parseSV(txt) {
+  const lines = txt.trim().split(/\r?\n/);
+  const [err, header] = parseSVheader(lines[0]);
+  if (err) return ["header error: " + err];
+  const items = [];
+  for (let i = 1; i < lines.length; i++) {
+    const line = lines[i];
+    if (line[0] == "#") continue;
+    const [e, m] = parseSVline(line, header);
+    if (e) return ["line " + (i + 1) + " error: " + e];
+    items.push(m);
+  }
+  return [null, header, items];
+}
+function parseSVheader(line) {
+  const header = line.toLowerCase().split("	");
+  if (header.length <= 1) return "invalid file header for fusions";
+  const htry = (...lst) => {
+    for (const a of lst) {
+      const j = header.indexOf(a);
+      if (j != -1) return j;
+    }
+    return -1;
+  };
+  let i = htry("chr_a", "chr1", "chra");
+  if (i == -1) return "chr_A missing from header";
+  header[i] = "chr1";
+  i = htry("chr_b", "chr2", "chrb");
+  if (i == -1) return "chr_B missing from header";
+  header[i] = "chr2";
+  i = htry("pos_a", "position_a", "position1", "posa");
+  if (i == -1) return "pos_a missing from header";
+  header[i] = "position1";
+  i = htry("pos_b", "position_b", "position2", "posb");
+  if (i == -1) return "pos_b missing from header";
+  header[i] = "position2";
+  i = htry("strand_a", "orta", "orienta");
+  if (i == -1) return "strand_a missing from header";
+  header[i] = "strand1";
+  i = htry("strand_b", "ortb", "orientb");
+  if (i == -1) return "strand_b missing from header";
+  header[i] = "strand2";
+  i = htry("numreadsa");
+  if (i != -1) header[i] = "reads1";
+  i = htry("numreadsb");
+  if (i != -1) header[i] = "reads2";
+  return [null, header];
+}
+function parseSVline(line, header) {
+  const lst = line.split("	");
+  const m = {};
+  for (let j = 0; j < header.length; j++) {
+    m[header[j]] = lst[j];
+  }
+  if (!m.chr1) return ["missing chr1"];
+  if (m.chr1.toLowerCase().indexOf("chr") != 0) {
+    m.chr1 = "chr" + m.chr1;
+  }
+  if (!m.chr2) return ["missing chr2"];
+  if (m.chr2.toLowerCase().indexOf("chr") != 0) {
+    m.chr2 = "chr" + m.chr2;
+  }
+  if (!m.position1) return ["missing position1"];
+  let v = Number.parseInt(m.position1);
+  if (Number.isNaN(v) || v <= 0) return ["position1 invalid value"];
+  m.position1 = v;
+  if (!m.position2) return ["missing position2"];
+  v = Number.parseInt(m.position2);
+  if (Number.isNaN(v) || v <= 0) return ["position2 invalid value"];
+  m.position2 = v;
+  if (m.reads1) {
+    v = Number.parseInt(m.reads1);
+    if (Number.isNaN(v)) return ["reads1 invalid value"];
+    m.reads1 = v;
+  }
+  if (m.reads2) {
+    v = Number.parseInt(m.reads2);
+    if (Number.isNaN(v)) return ["reads2 invalid value"];
+    m.reads2 = v;
+  }
+  return [null, m];
+}
+function hicparsestat(hic, j) {
+  if (!j) return "cannot stat hic file";
+  hic.normalization = j.normalization;
+  hic.version = j.version;
+  if (!j.Chromosomes) return "Chromosomes not found in file stat";
+  if (!Array.isArray(j.chrorder)) return ".chrorder[] missing";
+  if (j.chrorder.length == 0) return ".chrorder[] empty array";
+  hic.chrorder = j.chrorder;
+  if (!j["Base pair-delimited resolutions"]) return "Base pair-delimited resolutions not found in file stat";
+  if (!Array.isArray(j["Base pair-delimited resolutions"])) return "Base pair-delimited resolutions should be array";
+  hic.bpresolution = j["Base pair-delimited resolutions"];
+  if (!j["Fragment-delimited resolutions"]) return "Fragment-delimited resolutions not found in file stat";
+  if (!Array.isArray(j["Fragment-delimited resolutions"])) return "Fragment-delimited resolutions is not array";
+  hic.fragresolution = j["Fragment-delimited resolutions"];
+  const chrlst = [];
+  for (const chr in j.Chromosomes) {
+    chrlst.push(chr);
+  }
+  const [nochrcount, haschrcount] = contigNameNoChr2(hic.genome, chrlst);
+  if (nochrcount + haschrcount == 0) return "chromosome names do not match with genome build";
+  if (nochrcount > 0) {
+    hic.nochr = true;
+    for (let i = 0; i < hic.chrorder.length; i++) hic.chrorder[i] = "chr" + hic.chrorder[i];
+  }
+  hic.chrlst = [];
+  for (const chr of hic.genome.majorchrorder) {
+    const c2 = hic.nochr ? chr.replace("chr", "") : chr;
+    if (chrlst.indexOf(c2) != -1) {
+      hic.chrlst.push(chr);
+    }
+  }
+}
+function hicparsefragdata(items) {
+  const id2coord = /* @__PURE__ */ new Map();
+  let min = null, max;
+  for (const i of items) {
+    if (!i.rest || !i.rest[0]) {
+      return ["items[].rest data problem"];
+    }
+    const id = Number.parseInt(i.rest[0]);
+    if (Number.isNaN(id)) {
+      return [i.start + "." + i.stop + " invalid fragment id: " + i.rest[0]];
+    }
+    id2coord.set(id, [i.start, i.stop]);
+    if (min == null) {
+      min = id;
+      max = id;
+    } else {
+      min = Math.min(min, id);
+      max = Math.max(max, id);
+    }
+  }
+  return [null, id2coord, min, max];
+}
+
+export {
+  hicParseFile,
+  hicparsestat,
+  hicparsefragdata
+};
+//# sourceMappingURL=chunk-O3XTTT4K.js.map

package/dist/chunk-PGMHFGW7.js

@@ -0,0 +1,102 @@
+import {
+  getPlotConfig
+} from "./chunk-QQH32U6O.js";
+import {
+  fillTermWrapper
+} from "./chunk-3YQYWX4M.js";
+import {
+  dictionaryNumericTypes,
+  numericTypes
+} from "./chunk-JHLF2ZC2.js";
+import {
+  copyMerge
+} from "./chunk-HBNB5TRH.js";
+
+// plots/matrix/hierCluster.config.js
+async function getPlotConfig2(opts = {}, app) {
+  opts.chartType = "hierCluster";
+  if (dictionaryNumericTypes.has(opts.dataType) || opts.dataType == "numericDictTerm") {
+    const grp = opts.termgroups?.find((g) => g.type == "hierCluster");
+    for (const tw of grp?.lst || []) tw.q = { ...tw.q, mode: "continuous" };
+  }
+  const config = await getPlotConfig(opts, app);
+  delete config.genes;
+  if (config.dataType == "numericDictTerm") {
+    const lst = config.termgroups?.find((g) => g.type == "hierCluster")?.lst;
+    config.dataType = lst?.[0]?.term?.type || "float";
+  }
+  config.settings.hierCluster = {
+    /* type of data used for clustering
+    exciting todo: (to introduce new dt values)
+    - gene dependency
+    - numeric dic term
+    - non-gene genomic stuff that resolves into numeric quantities (cpg meth)
+    - metabolite
+    */
+    dataType: config.dataType,
+    // adjust the default group name based on automatically detected term types
+    // Done in matrix.cells.js: setHierClusterCellProps
+    // termGroupName: 'Expression',
+    clusterSamples: true,
+    clusterRows: true,
+    clusterMethod: "average",
+    // complete
+    distanceMethod: "euclidean",
+    zScoreCap: 5,
+    zScoreTransformation: true,
+    xDendrogramHeight: 100,
+    yDendrogramHeight: 200,
+    colorScale: "blueWhiteRed"
+  };
+  const overrides = app.vocabApi.termdbConfig.hierCluster || {};
+  const numericDictTermClusterOverrides = dictionaryNumericTypes.has(config.dataType) && app.vocabApi.termdbConfig.numericDictTermCluster ? app.vocabApi.termdbConfig.numericDictTermCluster : {};
+  copyMerge(
+    config.settings.hierCluster,
+    overrides.settings,
+    opts.settings?.hierCluster || {},
+    numericDictTermClusterOverrides.settings
+  );
+  {
+    const c = config.settings.hierCluster.colorScale;
+    if (!c) throw "colorScale missing";
+  }
+  config.settings.matrix.collabelpos = "top";
+  const termGroupName = config.settings.hierCluster.termGroupName;
+  const hcTermGroup = config.termgroups.find((g) => g.type == "hierCluster" || g.name == termGroupName) || {
+    name: termGroupName
+  };
+  hcTermGroup.type = "hierCluster";
+  if (!hcTermGroup.lst?.length) {
+    if (!Array.isArray(opts.terms)) throw "opts.terms[] not array (may show geneset edit ui)";
+    const promises = [];
+    for (const i of opts.terms) {
+      const tw = i.term ? i : { term: i };
+      if (!tw.term.type) {
+        if (config.dataType && numericTypes.has(config.dataType)) {
+          tw.term.type = config.dataType;
+        } else {
+          throw `term type missing and cannot be assigned by dataType '${config.dataType}'`;
+        }
+      } else if (!numericTypes.has(tw.term.type)) {
+        throw "term type is not numeric";
+      } else if (config.dataType && !canTermBeInHierGrp(config.dataType, tw.term.type)) {
+        throw `cannot have term type ${tw.term.type} in ${config.dataType} term group`;
+      }
+      if (dictionaryNumericTypes.has(tw.term.type)) tw.q = { ...tw.q, mode: "continuous" };
+      promises.push(fillTermWrapper(tw, app.vocabApi));
+    }
+    hcTermGroup.lst = await Promise.all(promises);
+    if (config.termgroups.indexOf(hcTermGroup) == -1) config.termgroups.unshift(hcTermGroup);
+  }
+  config.settings.matrix.maxSample = 1e5;
+  return config;
+}
+function canTermBeInHierGrp(grpType, twType) {
+  if (dictionaryNumericTypes.has(grpType) && dictionaryNumericTypes.has(twType)) return true;
+  return twType == grpType;
+}
+
+export {
+  getPlotConfig2 as getPlotConfig
+};
+//# sourceMappingURL=chunk-PGMHFGW7.js.map

package/dist/chunk-PUIQ3UF7.js

@@ -0,0 +1,176 @@
+import {
+  appInit
+} from "./chunk-JUPLOOCO.js";
+import {
+  showErrorsWithCounter
+} from "./chunk-3YQYWX4M.js";
+import {
+  mclass
+} from "./chunk-WJLRU5OF.js";
+
+// plots/disco/launch.adhoc.ts
+async function launch(arg, genomeObj, holder) {
+  const [mlst, errors] = await getMlst(arg);
+  if (errors?.length) {
+    return showErrorsWithCounter(errors, holder);
+  }
+  const opts = {
+    holder,
+    vocabApi: {
+      // api is required by plot.app.js, so create a mock one for the adhoc data
+      vocab: { terms: [] },
+      main: () => {
+        return;
+      },
+      getTermdbConfig: () => {
+        return {};
+      }
+    },
+    state: {
+      args: {
+        data: mlst,
+        genome: genomeObj
+      },
+      plots: [
+        {
+          chartType: "Disco",
+          subfolder: "disco",
+          extension: "ts",
+          /** NOTE: Users should only override the settings in the default
+           * settings.Disco:{}, not the entire settings:{} object.*/
+          overrides: { Disco: arg?.settings || {} }
+        }
+      ]
+    }
+  };
+  const plotAppApi = await appInit(opts);
+  return plotAppApi;
+}
+async function getMlst(arg) {
+  if (Array.isArray(arg.mlst)) {
+    return [arg.mlst, null];
+  }
+  const mlst = [];
+  const errors = [];
+  if (arg.snvText) parseSnvText(arg.snvText, mlst, errors);
+  if (arg.svText) parseSvText(arg.svText, mlst, errors);
+  if (arg.cnvText) parseCnvText(arg.cnvText, mlst, errors);
+  return [mlst, errors];
+}
+function parseSnvText(text, mlst, errors) {
+  for (const line of text.trim().split("\n")) {
+    const l = line.trim().split("	");
+    if (![5, 7, 9].includes(l.length)) {
+      errors.push("snv input not equal to 5, 7, or 9 columns");
+      continue;
+    }
+    let m;
+    try {
+      m = {
+        dt: 1,
+        chr: l[0],
+        position: Number(l[1]),
+        gene: l[2],
+        mname: l[3],
+        class: validateMutation(l[4], errors)
+      };
+      const vafs = parseOptionalVafs(l, errors);
+      if (vafs.length) m.vafs = vafs;
+    } catch (e) {
+      errors.push(e);
+      continue;
+    }
+    mlst.push(m);
+  }
+}
+function parseOptionalVafs(line, errors) {
+  const vafs = [];
+  const addVaf = (id, totalIndex, altIndex) => {
+    if (line.length <= altIndex) return;
+    const totalCount = Number(line[totalIndex]);
+    const altCount = Number(line[altIndex]);
+    if (!Number.isInteger(totalCount) || !Number.isInteger(altCount) || totalCount <= 0 || altCount < 0 || altCount > totalCount) {
+      errors.push(`${id} total/alt counts must be integers with total > 0, alt >= 0, and alt cannot exceed total`);
+      return;
+    }
+    vafs.push({ id, totalCount, altCount });
+  };
+  addVaf("DNA", 5, 6);
+  addVaf("RNA", 7, 8);
+  return vafs;
+}
+function parseSvText(text, mlst, errors) {
+  for (const line of text.trim().split("\n")) {
+    const l = line.trim().split("	");
+    if (l.length < 4 || l.length > 6) {
+      errors.push("sv input not equal to 4 or 6 columns");
+      continue;
+    }
+    let m;
+    try {
+      const length = l.length;
+      if (length == 4) {
+        m = {
+          dt: 2,
+          chrA: l[0],
+          posA: Number(l[1]),
+          chrB: l[2],
+          posB: Number(l[3])
+        };
+      } else {
+        m = {
+          dt: 2,
+          chrA: l[0],
+          posA: Number(l[1]),
+          geneA: l[2],
+          chrB: l[3],
+          posB: Number(l[4]),
+          geneB: l[5]
+        };
+      }
+    } catch (e) {
+      errors.push(e);
+      continue;
+    }
+    mlst.push(m);
+  }
+}
+function parseCnvText(text, mlst, errors) {
+  for (const line of text.trim().split("\n")) {
+    const l = line.trim().split("	");
+    if (l.length != 4) {
+      errors.push("cnv input not equal to 4 columns");
+      continue;
+    }
+    let m;
+    try {
+      m = {
+        dt: 4,
+        chr: l[0],
+        start: Number(l[1]),
+        stop: Number(l[2]),
+        value: Number(l[3])
+      };
+    } catch (e) {
+      errors.push(e);
+      continue;
+    }
+    mlst.push(m);
+  }
+}
+function validateMutation(mutation, errors) {
+  const mut2check = mutation.toLowerCase();
+  const foundMutation = Object.values(mclass).find(
+    (m) => m.key.toLowerCase() === mut2check || m.label.toLowerCase() === mut2check
+  );
+  if (foundMutation) {
+    return foundMutation.key;
+  } else {
+    errors.push(`Invalid mutation class: ${mutation}`);
+  }
+}
+
+export {
+  launch
+};
+//# sourceMappingURL=chunk-PUIQ3UF7.js.map

package/dist/chunk-QBGMD5UV.js

@@ -0,0 +1,50 @@
+import {
+  runproteinpaint
+} from "./chunk-AFEPQY5B.js";
+
+// test/front.helpers.js
+var serverData = /* @__PURE__ */ Object.create(null);
+function getRunPp(appname = "", defaultArgs = {}, _host = "http://localhost:3000") {
+  const host = window.testHost || _host;
+  const arg = {
+    host: window.testHost || _host,
+    noheader: 1,
+    nobox: true,
+    debug: 1,
+    norecover: true
+  };
+  if (appname == "mass" || appname == "termdb") {
+    defaultArgs.debounceInterval = 0;
+  }
+  if (appname) arg[appname] = defaultArgs;
+  else copyMerge(arg, defaultArgs);
+  const argStr = JSON.stringify(arg);
+  return function runpp(overrides = {}) {
+    const argCopy = JSON.parse(argStr);
+    if (appname) copyMerge(argCopy[appname], overrides);
+    else copyMerge(argCopy, overrides);
+    if (appname && defaultArgs.fetchOpts) {
+      argCopy[appname].fetchOpts = defaultArgs.fetchOpts;
+    }
+    return runproteinpaint(Object.assign(argCopy, { serverData }));
+  };
+}
+function copyMerge(base, ...args) {
+  const target = typeof base == "string" ? JSON.parse(base) : base;
+  for (const arg of args) {
+    if (arg) {
+      const source = typeof base == "string" ? JSON.parse(JSON.stringify(arg)) : arg;
+      for (const key in source) {
+        if (key === "__proto__" || key === "constructor" || key === "prototype") continue;
+        if (!target[key] || Array.isArray(target[key]) || typeof target[key] !== "object") target[key] = source[key];
+        else copyMerge(target[key], source[key]);
+      }
+    }
+  }
+  return target;
+}
+
+export {
+  getRunPp
+};
+//# sourceMappingURL=chunk-QBGMD5UV.js.map