@jbrowse/plugin-alignments 2.16.1 → 2.17.0

package/dist/PileupRenderer/renderModifications.js

@@ -2,47 +2,44 @@
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.renderModifications = renderModifications;
 const util_1 = require("@jbrowse/core/util");
-const MismatchParser_1 = require("../MismatchParser");
-const util_2 = require("../util");
-const util_3 = require("./util");
-const colord_1 = require("@jbrowse/core/util/colord");
+// locals
+const util_2 = require("./util");
+const util_3 = require("../shared/util");
+const getMaximumModificationAtEachPosition_1 = require("../shared/getMaximumModificationAtEachPosition");
 // render modifications stored in MM tag in BAM
-//
-// ML stores probabilities as array of numerics and MP is scaled phred scores
-// https://github.com/samtools/hts-specs/pull/418/files#diff-e765c6479316309f56b636f88189cdde8c40b854c7bdcce9ee7fe87a4e76febcR596
-//
-// if we have ML or Ml, it is an 8bit probability, divide by 255
-//
-// if we have MP or Mp it is phred scaled ASCII, which can go up to 90 but
-// has very high likelihood basecalls at that point, we really only care
-// about low qual calls <20 approx
-function renderModifications({ ctx, feat, region, bpPerPx, renderArgs, canvasWidth, }) {
+function renderModifications({ ctx, feat, region, bpPerPx, renderArgs, canvasWidth, cigarOps, }) {
+    var _a, _b, _c;
     const { feature, topPx, heightPx } = feat;
-    const { modificationTagMap = {} } = renderArgs;
+    const { colorBy, visibleModifications = {} } = renderArgs;
     const seq = feature.get('seq');
     if (!seq) {
         return;
     }
-    const mm = (0, util_2.getTagAlt)(feature, 'MM', 'Mm') || '';
-    const cigar = feature.get('CIGAR');
     const start = feature.get('start');
-    const strand = feature.get('strand');
-    const cigarOps = (0, MismatchParser_1.parseCigar)(cigar);
-    const probabilities = (0, MismatchParser_1.getModificationProbabilities)(feature);
-    const modifications = (0, MismatchParser_1.getModificationPositions)(mm, seq, strand);
-    // probIndex applies across multiple modifications e.g.
-    let probIndex = 0;
-    for (const { type, positions } of modifications) {
-        const col = modificationTagMap[type] || 'black';
-        const base = (0, colord_1.colord)(col);
-        for (const readPos of (0, MismatchParser_1.getNextRefPos)(cigarOps, positions)) {
-            const r = start + readPos;
-            const [leftPx, rightPx] = (0, util_1.bpSpanPx)(r, r + 1, region, bpPerPx);
-            const prob = (probabilities === null || probabilities === void 0 ? void 0 : probabilities[probIndex]) || 0;
-            const c = prob !== 1 ? base.alpha(prob + 0.1).toHslString() : col;
+    const isolatedModification = (_a = colorBy === null || colorBy === void 0 ? void 0 : colorBy.modifications) === null || _a === void 0 ? void 0 : _a.isolatedModification;
+    const twoColor = (_b = colorBy === null || colorBy === void 0 ? void 0 : colorBy.modifications) === null || _b === void 0 ? void 0 : _b.twoColor;
+    (_c = (0, getMaximumModificationAtEachPosition_1.getMaxProbModAtEachPosition)(feature, cigarOps)) === null || _c === void 0 ? void 0 : _c.forEach(({ allProbs, prob, type }, pos) => {
+        const r = start + pos;
+        const [leftPx, rightPx] = (0, util_1.bpSpanPx)(r, r + 1, region, bpPerPx);
+        const mod = visibleModifications[type];
+        if (!mod) {
+            console.warn(`${type} not known yet`);
+            return;
+        }
+        if (isolatedModification && mod.type !== isolatedModification) {
+            return;
+        }
+        const col = mod.color || 'black';
+        const s = 1 - (0, util_1.sum)(allProbs);
+        if (twoColor && s > (0, util_1.max)(allProbs)) {
+            const c = (0, util_3.alphaColor)('blue', s);
             const w = rightPx - leftPx + 0.5;
-            (0, util_3.fillRect)(ctx, leftPx, topPx, w, heightPx, canvasWidth, c);
-            probIndex++;
+            (0, util_2.fillRect)(ctx, leftPx, topPx, w, heightPx, canvasWidth, c);
         }
-    }
+        else {
+            const c = (0, util_3.alphaColor)(col, prob);
+            const w = rightPx - leftPx + 0.5;
+            (0, util_2.fillRect)(ctx, leftPx, topPx, w, heightPx, canvasWidth, c);
+        }
+    });
 }

package/dist/PileupRenderer/renderPerBaseLettering.d.ts

@@ -1,6 +1,6 @@
 import { Region } from '@jbrowse/core/util';
 import { LayoutFeature } from './util';
-export declare function renderPerBaseLettering({ ctx, feat, region, bpPerPx, colorForBase, contrastForBase, charWidth, charHeight, canvasWidth, }: {
+export declare function renderPerBaseLettering({ ctx, feat, region, bpPerPx, colorForBase, contrastForBase, charWidth, charHeight, canvasWidth, cigarOps, }: {
     ctx: CanvasRenderingContext2D;
     feat: LayoutFeature;
     region: Region;
@@ -10,4 +10,5 @@ export declare function renderPerBaseLettering({ ctx, feat, region, bpPerPx, col
     charWidth: number;
     charHeight: number;
     canvasWidth: number;
+    cigarOps: string[];
 }): void;

package/dist/PileupRenderer/renderPerBaseLettering.js

@@ -2,13 +2,11 @@
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.renderPerBaseLettering = renderPerBaseLettering;
 const util_1 = require("@jbrowse/core/util");
-const MismatchParser_1 = require("../MismatchParser");
 const util_2 = require("./util");
-function renderPerBaseLettering({ ctx, feat, region, bpPerPx, colorForBase, contrastForBase, charWidth, charHeight, canvasWidth, }) {
+function renderPerBaseLettering({ ctx, feat, region, bpPerPx, colorForBase, contrastForBase, charWidth, charHeight, canvasWidth, cigarOps, }) {
     const heightLim = charHeight - 2;
     const { feature, topPx, heightPx } = feat;
     const seq = feature.get('seq');
-    const cigarOps = (0, MismatchParser_1.parseCigar)(feature.get('CIGAR'));
     const w = 1 / bpPerPx;
     const start = feature.get('start');
     let soffset = 0;

package/dist/PileupRenderer/renderPerBaseQuality.d.ts

@@ -1,9 +1,10 @@
 import { Region } from '@jbrowse/core/util';
 import { LayoutFeature } from './util';
-export declare function renderPerBaseQuality({ ctx, feat, region, bpPerPx, canvasWidth, }: {
+export declare function renderPerBaseQuality({ ctx, feat, region, bpPerPx, canvasWidth, cigarOps, }: {
     ctx: CanvasRenderingContext2D;
     feat: LayoutFeature;
     region: Region;
     bpPerPx: number;
     canvasWidth: number;
+    cigarOps: string[];
 }): void;

package/dist/PileupRenderer/renderPerBaseQuality.js

@@ -2,13 +2,11 @@
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.renderPerBaseQuality = renderPerBaseQuality;
 const util_1 = require("@jbrowse/core/util");
-const MismatchParser_1 = require("../MismatchParser");
 const util_2 = require("./util");
-function renderPerBaseQuality({ ctx, feat, region, bpPerPx, canvasWidth, }) {
+function renderPerBaseQuality({ ctx, feat, region, bpPerPx, canvasWidth, cigarOps, }) {
     const { feature, topPx, heightPx } = feat;
     const qual = feature.get('qual') || '';
     const scores = qual.split(' ').map(val => +val);
-    const cigarOps = (0, MismatchParser_1.parseCigar)(feature.get('CIGAR'));
     const w = 1 / bpPerPx;
     const start = feature.get('start');
     let soffset = 0; // sequence offset

package/dist/PileupRenderer/renderSoftClipping.js

@@ -3,9 +3,8 @@ Object.defineProperty(exports, "__esModule", { value: true });
 exports.renderSoftClipping = renderSoftClipping;
 const configuration_1 = require("@jbrowse/core/configuration");
 const util_1 = require("@jbrowse/core/util");
-// locals
-const MismatchParser_1 = require("../MismatchParser");
 const util_2 = require("./util");
+const MismatchParser_1 = require("../MismatchParser");
 function renderSoftClipping({ ctx, feat, renderArgs, config, theme, colorForBase, canvasWidth, }) {
     const { feature, topPx, heightPx } = feat;
     const { regions, bpPerPx } = renderArgs;
@@ -19,12 +18,13 @@ function renderSoftClipping({ ctx, feat, renderArgs, config, theme, colorForBase
         return;
     }
     const heightLim = charHeight - 2;
-    const CIGAR = (0, MismatchParser_1.parseCigar)(feature.get('CIGAR'));
     let seqOffset = 0;
     let refOffset = 0;
-    for (let i = 0; i < CIGAR.length; i += 2) {
-        const op = CIGAR[i + 1];
-        const len = +CIGAR[i];
+    const CIGAR = feature.get('CIGAR');
+    const cigarOps = (0, MismatchParser_1.parseCigar)(CIGAR);
+    for (let i = 0; i < cigarOps.length; i += 2) {
+        const op = cigarOps[i + 1];
+        const len = +cigarOps[i];
         if (op === 'S') {
             for (let k = 0; k < len; k++) {
                 const base = seq[seqOffset + k];

package/dist/PileupRenderer/sortUtil.d.ts

@@ -1,8 +1,3 @@
 import { Feature } from '@jbrowse/core/util';
-interface SortObject {
-    pos: number;
-    type: string;
-    tag?: string;
-}
-export declare const sortFeature: (features: Map<string, Feature>, sortedBy: SortObject) => Map<string, Feature>;
-export {};
+import { SortedBy } from '../shared/types';
+export declare const sortFeature: (features: Map<string, Feature>, sortedBy: SortedBy) => Map<string, Feature>;

package/dist/PileupRenderer/sortUtil.js

@@ -43,29 +43,29 @@ const sortFeature = (features, sortedBy) => {
         // first sort all mismatches, then all reference bases at the end
         case 'Base pair': {
             const baseSortArray = [];
-            featuresInCenterLine.forEach(feature => {
+            for (const feature of featuresInCenterLine) {
                 const mismatches = feature.get('mismatches');
-                mismatches.forEach(mismatch => {
+                for (const m of mismatches) {
                     const start = feature.get('start');
-                    const offset = start + mismatch.start + 1;
-                    const consuming = mismatch.type === 'insertion' || mismatch.type === 'softclip';
-                    const len = consuming ? 0 : mismatch.length;
+                    const offset = start + m.start + 1;
+                    const consuming = m.type === 'insertion' || m.type === 'softclip';
+                    const len = consuming ? 0 : m.length;
                     if (pos >= offset && pos < offset + len) {
-                        baseSortArray.push([feature.id(), mismatch]);
+                        baseSortArray.push([feature.id(), m]);
                     }
-                });
-            });
+                }
+            }
             const baseMap = new Map(baseSortArray);
             featuresInCenterLine.sort((a, b) => {
                 const aMismatch = baseMap.get(a.id());
                 const bMismatch = baseMap.get(b.id());
                 const acode = bMismatch === null || bMismatch === void 0 ? void 0 : bMismatch.base.toUpperCase();
                 const bcode = aMismatch === null || aMismatch === void 0 ? void 0 : aMismatch.base.toUpperCase();
-                if (acode === bcode && acode === '*') {
-                    // @ts-expect-error
-                    return aMismatch.length - bMismatch.length;
-                }
-                return ((acode ? acode.charCodeAt(0) : 0) - (bcode ? bcode.charCodeAt(0) : 0));
+                return acode === bcode && acode === '*'
+                    ? // @ts-expect-error
+                        aMismatch.length - bMismatch.length
+                    : (acode ? acode.charCodeAt(0) : 0) -
+                        (bcode ? bcode.charCodeAt(0) : 0);
             });
             break;
         }

package/dist/SNPCoverageAdapter/SNPCoverageAdapter.js

@@ -9,8 +9,8 @@ const rxjs_1 = require("@jbrowse/core/util/rxjs");
 const operators_1 = require("rxjs/operators");
 const rxjs_2 = require("rxjs");
 // locals
-const generateCoverageBins_1 = __importDefault(require("./generateCoverageBins"));
 const util_1 = require("../util");
+const generateCoverageBins_1 = require("./generateCoverageBins");
 class SNPCoverageAdapter extends BaseAdapter_1.BaseFeatureDataAdapter {
     async configure() {
         var _a, _b;
@@ -38,14 +38,19 @@ class SNPCoverageAdapter extends BaseAdapter_1.BaseFeatureDataAdapter {
     getFeatures(region, opts = {}) {
         return (0, rxjs_1.ObservableCreate)(async (observer) => {
             const { subadapter } = await this.configure();
-            const feats = await (0, rxjs_2.firstValueFrom)(subadapter.getFeatures(region, opts).pipe((0, operators_1.toArray)()));
-            const { bins, skipmap } = await (0, generateCoverageBins_1.default)(feats, region, opts, arg => this.fetchSequence(arg));
+            const features = await (0, rxjs_2.firstValueFrom)(subadapter.getFeatures(region, opts).pipe((0, operators_1.toArray)()));
+            const { bins, skipmap } = await (0, generateCoverageBins_1.generateCoverageBins)({
+                features,
+                region,
+                opts,
+                fetchSequence: (region) => this.fetchSequence(region),
+            });
             bins.forEach((bin, index) => {
                 const start = region.start + index;
                 observer.next(new simpleFeature_1.default({
                     id: `${this.id}-${start}`,
                     data: {
-                        score: bin.total,
+                        score: bin.depth,
                         snpinfo: bin,
                         start,
                         end: start + 1,
@@ -63,7 +68,7 @@ class SNPCoverageAdapter extends BaseAdapter_1.BaseFeatureDataAdapter {
                         end: skip.end,
                         strand: skip.strand,
                         score: skip.score,
-                        xs: skip.xs,
+                        effectiveStrand: skip.effectiveStrand,
                     },
                 }));
             });

package/dist/SNPCoverageAdapter/generateCoverageBins.d.ts

@@ -1,13 +1,17 @@
 import { AugmentedRegion as Region } from '@jbrowse/core/util/types';
-import { Feature } from '@jbrowse/core/util/simpleFeature';
-import { Bin, SkipMap } from './util';
-export default function generateCoverageBins(features: Feature[], region: Region, opts: {
+import { Feature } from '@jbrowse/core/util';
+import { ColorBy, PreBaseCoverageBin, SkipMap } from '../shared/types';
+interface Opts {
     bpPerPx?: number;
-    colorBy?: {
-        type: string;
-        tag?: string;
-    };
-}, fetchSequence: (arg: Region) => Promise<string>): Promise<{
-    bins: Bin[];
+    colorBy?: ColorBy;
+}
+export declare function generateCoverageBins({ fetchSequence, features, region, opts, }: {
+    features: Feature[];
+    region: Region;
+    opts: Opts;
+    fetchSequence: (arg: Region) => Promise<string>;
+}): Promise<{
+    bins: PreBaseCoverageBin[];
     skipmap: SkipMap;
 }>;
+export {};