jbrowse-plugin-graphgenomeviewer 4.0.19 → 4.0.21
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/README.md +2 -1
- package/dist/chunks/{GbzBaseSyntenyAdapter-KKSNCZXY.js → GbzBaseSyntenyAdapter-7ZKFAEJG.js} +2 -2
- package/dist/chunks/{GbzBaseSyntenyAdapter-KKSNCZXY.js.map → GbzBaseSyntenyAdapter-7ZKFAEJG.js.map} +2 -2
- package/dist/chunks/{GraphGenomeView-37F7MLFQ.js → GraphGenomeView-CLUMXH4K.js} +2 -2
- package/dist/chunks/{GraphTrackSettingsDialog-2QV6TPX7.js → GraphTrackSettingsDialog-QF5OMTCH.js} +2 -2
- package/dist/chunks/{LinearGraphDisplay-CJHA56RB.js → LinearGraphDisplay-4CMVCEM5.js} +2 -2
- package/dist/chunks/chunk-3FGJLRG2.js +2 -0
- package/dist/chunks/chunk-3FGJLRG2.js.map +7 -0
- package/dist/chunks/chunk-3SZJ4AFI.js +3 -0
- package/dist/chunks/chunk-3SZJ4AFI.js.map +7 -0
- package/dist/chunks/{chunk-5E7TVPCI.js → chunk-CLRXQ6BM.js} +2 -2
- package/dist/chunks/chunk-VIFXYTIS.js +2 -0
- package/dist/chunks/chunk-VIFXYTIS.js.map +7 -0
- package/dist/jbrowse-plugin-graphgenomeviewer.esm.js +2 -2
- package/dist/jbrowse-plugin-graphgenomeviewer.esm.js.map +3 -3
- package/package.json +2 -2
- package/src/GbzBaseSyntenyAdapter/GbzBaseSyntenyAdapter.ts +1 -1
- package/src/GetSubgraph.ts +3 -0
- package/src/GraphGenomeView/components/BubbleHalos.tsx +3 -4
- package/src/GraphGenomeView/components/GraphCanvas.tsx +51 -15
- package/src/GraphGenomeView/components/LabelLayer.tsx +5 -6
- package/src/GraphGenomeView/components/WalkSelect.tsx +26 -13
- package/src/GraphGenomeView/model.test.ts +48 -0
- package/src/GraphGenomeView/model.ts +75 -21
- package/src/GraphGenomeView/viewModel.ts +15 -0
- package/src/LinearGraphDisplay/LinearGraphDisplay.test.ts +23 -2
- package/src/LinearGraphDisplay/model.ts +7 -2
- package/src/version.ts +1 -1
- package/dist/chunks/chunk-BI3M6FZD.js +0 -2
- package/dist/chunks/chunk-BI3M6FZD.js.map +0 -7
- package/dist/chunks/chunk-IUIIYSGF.js +0 -2
- package/dist/chunks/chunk-IUIIYSGF.js.map +0 -7
- package/dist/chunks/chunk-MNOQXDMY.js +0 -3
- package/dist/chunks/chunk-MNOQXDMY.js.map +0 -7
- /package/dist/chunks/{GraphGenomeView-37F7MLFQ.js.map → GraphGenomeView-CLUMXH4K.js.map} +0 -0
- /package/dist/chunks/{GraphTrackSettingsDialog-2QV6TPX7.js.map → GraphTrackSettingsDialog-QF5OMTCH.js.map} +0 -0
- /package/dist/chunks/{LinearGraphDisplay-CJHA56RB.js.map → LinearGraphDisplay-4CMVCEM5.js.map} +0 -0
- /package/dist/chunks/{chunk-5E7TVPCI.js.map → chunk-CLRXQ6BM.js.map} +0 -0
|
@@ -1,2 +0,0 @@
|
|
|
1
|
-
import{d as ct}from"./chunk-UQEPPJEO.js";import{c as lt}from"./chunk-SCMIJAAX.js";import{i as Ne,n as ut,o as dt,p as yn,q as xn,t as ft,u as pt,y as ht,z as Mn}from"./chunk-IUIIYSGF.js";import{a as Ge,b as _,c as Qe,d as we,e as Se,f as it,h as Ze}from"./chunk-EFVAQ2KM.js";import{a as Ke,c as st,d as at}from"./chunk-KULNKXPT.js";import{a as Ue}from"./chunk-XEJ47DLD.js";import{a as Ve,b as ie}from"./chunk-FY3BTKY4.js";var bt=Ve((ol,mt)=>{mt.exports=JBrowseExports["@mui/material/SvgIcon"]});var gs=Ve((Vl,Ut)=>{Ut.exports=JBrowseExports.mobx});var pn=Ve((Jc,Xo)=>{Xo.exports=JBrowseExports["@jbrowse/core/util/colorBits"]});var Ji=Ve((Qu,hr)=>{hr.exports=JBrowseExports["@jbrowse/render-core/canvas2dUtils"]});var te=ie(bt(),1);function il(e){let n=t=>Math.round(t*100)/100;return e.map((t,r)=>`${r?"L":"M"}${n(t.x)},${n(t.y)}`).join("")}function gt(e,n,t,r,o){let s=Math.hypot(r-n,t-e);if(s===0)return[t,r];let a=(t-e)/s,l=(r-n)/s;return[t+o*a,r+o*l]}function Ln(e,n){let t=1-n,r=t*t*t,o=3*t*t*n,s=3*t*n*n,a=n*n*n;return{x:r*e.x0+o*e.cx0+s*e.cx1+a*e.x1,y:r*e.y0+o*e.cy0+s*e.cy1+a*e.y1}}function Nt(e){let n=e[Math.floor(e.length/2)];return n?Ln(n,.5):void 0}function Lt(e){let n=1/0,t=1/0,r=-1/0,o=-1/0;for(let s of e)n=Math.min(n,s.x0,s.cx0,s.cx1,s.x1),r=Math.max(r,s.x0,s.cx0,s.cx1,s.x1),t=Math.min(t,s.y0,s.cy0,s.cy1,s.y1),o=Math.max(o,s.y0,s.cy0,s.cy1,s.y1);return{minX:n,minY:t,maxX:r,maxY:o}}function Mr(e,n,t){let r=Ln(e,n),o=Ln(e,t),s=(t-n)/3,a=yt(e,n),l=yt(e,t);return{x0:r.x,y0:r.y,cx0:r.x+a.x*s,cy0:r.y+a.y*s,cx1:o.x-l.x*s,cy1:o.y-l.y*s,x1:o.x,y1:o.y}}function yt(e,n){let t=1-n,r=3*t*t,o=6*t*n,s=3*n*n;return{x:r*(e.cx0-e.x0)+o*(e.cx1-e.cx0)+s*(e.x1-e.cx1),y:r*(e.cy0-e.y0)+o*(e.cy1-e.cy0)+s*(e.y1-e.cy1)}}function Rt(e,n){let t=[];for(let r of e){let o=Math.hypot(r.cx0-r.x0,r.cy0-r.y0)+Math.hypot(r.cx1-r.cx0,r.cy1-r.cy0)+Math.hypot(r.x1-r.cx1,r.y1-r.cy1),s=Math.max(3,Math.round(o/n)|1);for(let a=0;a<s;a+=2)t.push([Mr(r,a/s,(a+1)/s)])}return t}function Et(e,n,t){let r=new Array(e.length);for(let o=0;o<e.length;o++){let s=e[o];r[o]={x0:s.x0+n,y0:s.y0+t,cx0:s.cx0+n,cy0:s.cy0+t,cx1:s.cx1+n,cy1:s.cy1+t,x1:s.x1+n,y1:s.y1+t}}return r}function wn(e,n){let t=e.length-1;return n==="end"?{at:e[t],inward:e[t-1]}:{at:e[0],inward:e[1]}}var wr=.5,Sr=80;function xt(e){let n=1/0,t=-1/0;for(let r of e)n=Math.min(n,r.x),t=Math.max(t,r.x);return{min:n,max:t}}function Mt(e,n){return Math.max(0,n.min-e,e-n.max)}function wt(e,n,t){let r=Mt(e[0].x,n),o=Mt(e[e.length-1].x,n);return r===o?t:r<o?"start":"end"}function Nr(e,n){return{from:wt(e,xt(n),"end"),to:wt(n,xt(e),"start")}}var Lr=.75,Rr=1.15,Er=.35,Ar=60;function Tr(e,n,t,r,o,s){let a=Math.hypot(t-e,r-n),l=a===0?1:(t-e)/a,i=a===0?0:(r-n)/a,u=0,f=1/0,c=-1/0;for(let m of o){let b=m.x-e,g=m.y-n,M=b*-i+g*l;Math.abs(M)>Math.abs(u)&&(u=M);let R=b*l+g*i;f=Math.min(f,R),c=Math.max(c,R)}let d=o.length?c-f:0,h=Math.max(Math.abs(u)*Rr/Lr,d*Er),p=s===void 0?h:Math.min(h,s);return u<0?-p:p}function Je({scaleX:e,scaleY:n}){return n/e}var Cr=3;function Sn(e,n,t){let r=n.x-e.x,o=(n.y-e.y)*t,s=Math.hypot(r,o);return s===0?void 0:{x:r/s,y:o/s}}function At(e,n,t,r,o=Cr){let s=Je(r),a=Sn(e[e.length-1],n[0],s)??Sn(e[0],e.at(-1),s)??Sn(n[0],n.at(-1),s)??{x:1,y:0},l=[];for(let i=0;i<t;i++){let u=(i-(t-1)/2)*o;l.push({x:-a.y*u,y:a.x*u/s})}return l}function Nn(e,n){return e.map(t=>({x:t.x,y:t.y*n}))}function Ir(e,n){for(let t of e)t.y0/=n,t.cy0/=n,t.cy1/=n,t.y1/=n;return e}var _r=12,vr=40;function St(e,n){let t=0;for(let s=1;s<e.length;s++)t+=Math.hypot(e[s].x-e[s-1].x,e[s].y-e[s-1].y);let r=t*n;return Math.min(vr,Math.max(_r,r/2))/Math.max(n,Number.EPSILON)}function Tt(e){return e.from!==e.to?!1:e.fromStrand&&e.toStrand&&e.fromStrand!==e.toStrand?e.fromStrand==="+"?"end":"start":!0}function Br(e,n,t,r){let{at:o,inward:s}=e,a=s?o.x-s.x:1,l=s?o.y-s.y:0,i=Math.hypot(a,l)||1,u=a/i,f=l/i,c=o.x+t,d=o.y+r,h=c+u*n/3,p=d+f*n/3,m=c+u*n,b=d+f*n,g=-f*n/2,M=u*n/2;return[{x0:c,y0:d,cx0:h,cy0:p,cx1:m+g,cy1:b+M,x1:m,y1:b},{x0:m,y0:b,cx0:m-g,cy0:b-M,cx1:h,cy1:p,x1:c,y1:d}]}function De(e,n,t,r,o,s,a=[],l){let i=s.scaleX,u=Je(s),f=l??(s.pixelRows?Ar:void 0);if(u!==1)return Ir(De(Nn(e,u),Nn(n,u),t,r,o*u,{scaleX:i,scaleY:i},Nn(a,u),f===void 0?void 0:f*u),u);if(typeof t=="string")return Br(wn(e,t),St(e,i),r,o);let c=t?{from:"end",to:"start"}:Nr(e,n),d=wn(e,c.from),h=wn(n,c.to),p=d.at,m=h.at,b=p.x+r,g=p.y+o,M=m.x+r,R=m.y+o;if(t){let N=1,S=0;if(d.inward){let U=p.x-d.inward.x,j=p.y-d.inward.y,V=Math.hypot(U,j);V>0&&(N=U/V,S=j/V)}let x=St(e,i),L=-S,A=N,C=(p.x+m.x)/2+r+L*x,P=(p.y+m.y)/2+o+A*x,D=b+N*x,$=g+S*x,W=M-N*x,H=R-S*x;return[{x0:b,y0:g,cx0:D,cy0:$,cx1:D+L*x,cy1:$+A*x,x1:C,y1:P},{x0:C,y0:P,cx0:W+L*x,cy0:H+A*x,cx1:W,cy1:H,x1:M,y1:R}]}else{let N=d.inward??{x:p.x-(m.x-p.x),y:p.y-(m.y-p.y)},S=h.inward??{x:m.x-(p.x-m.x),y:m.y-(p.y-m.y)},x=Math.hypot(M-b,R-g),L=(Z,J,re,_e)=>{let se=J.x-Z.x,ae=J.y-Z.y,ee=Math.hypot(se,ae),fe=ee===0?0:(re*se+_e*ae)/ee;return Math.max(0,Math.min(x*wr,Sr,fe*.5))},A=m.x-p.x,C=m.y-p.y,[P,D]=gt(N.x,N.y,p.x,p.y,L(N,p,A,C)),[$,W]=gt(S.x,S.y,m.x,m.y,L(S,m,-A,-C)),H=a.length?Tr(b,g,M,R,a,f):0,U=x===0?1:x,j=(M-b)/U,V=(R-g)/U,Xe=-V*H,B=j*H,q=j*Math.abs(H),K=V*Math.abs(H);return[{x0:b,y0:g,cx0:P+r+Xe-q,cy0:D+o+B-K,cx1:$+r+Xe+q,cy1:W+o+B+K,x1:M,y1:R}]}}var Or=1;function Rn(e,n){return n.flatMap(t=>e[t]??[])}function Ct(e,n,t){let r=e[n.from],o=e[n.to];return r?.length&&o?.length?De(r,o,n.from===n.to,0,0,t,Rn(e,n.bypassed)):void 0}function Pr(e){let n=new Map;for(let t of e.nodes)if(_(t)){let r=n.get(t.stable.refName);r?r.push(t):n.set(t.stable.refName,[t])}for(let t of n.values())t.sort((r,o)=>r.stable.start-o.stable.start);return n}function Gr(e,n,t){let r=[];if(e)for(let o=it(e,n);o<e.length;o++){let s=e[o];if(s.stable.start>=t)break;s.stable.start+s.length<=t&&r.push(s.id)}return r}function Dr(e,n,t,r){let{fromStrand:o,toStrand:s}=e;if(o===void 0||s===void 0)return!0;let a=o===(n.stable.strand??"+"),l=s===(t.stable.strand??"+");return r?a&&l:!a&&!l}function dl(e){let n=new Map(e.nodes.map(o=>[o.id,o])),t=Pr(e),r=[];for(let o=0;o<e.edges.length;o++){let s=e.edges[o],a=n.get(s.from),l=n.get(s.to);if(a!==void 0&&l!==void 0&&_(a)&&_(l)&&a.stable.refName===l.stable.refName){let i=a.stable.start<=l.stable.start,[u,f]=i?[a,l]:[l,a],c=u.stable.start+u.length,d=f.stable.start,h=d-c;if(h>=Or&&Dr(s,a,l,i)){let p=u.stable.refName;r.push({edgeIndex:o,from:s.from,to:s.to,refName:p,start:c,end:d,bp:h,bypassed:Gr(t.get(p),c,d)})}}}return r}function Le(e){return e.endsWith("-")?"-":"+"}function It(e){return e==="+"?"-":"+"}function $r(e){let n=(e.fromStrand??Le(e.from))===Le(e.from)?"end":"start",t=(e.toStrand??Le(e.to))===Le(e.to)?"start":"end";return{from:n,to:t}}function hl(e){let n=e.edges.map($r),t=new Map,r=(d,h,p)=>{let m=`${d}|${h}`,b=t.get(m);b?b.push(p):t.set(m,[p])};e.edges.forEach((d,h)=>{r(d.from,n[h].from,h),r(d.to,n[h].to,h)});let o=d=>{let h=t.get(`${d}|end`);if(h?.length!==1)return;let p=h[0],m=e.edges[p];return m.from===d&&m.to!==d&&n[p].from==="end"&&n[p].to==="start"&&t.get(`${m.to}|start`)?.length===1?p:void 0},s=new Set,a=new Map,l=new Map;for(let d of e.nodes){let h=t.get(`${d.id}|start`),p=h?.length===1?e.edges[h[0]].from:void 0;if(a.has(d.id)||p!==void 0&&o(p)===h[0])continue;let m=[d.id];for(;;){let b=o(m.at(-1)),g=b===void 0?void 0:e.edges[b].to;if(g===void 0||a.has(g))break;s.add(b),m.push(g)}for(let b of m)a.set(b,d.id);m.length>1&&l.set(d.id,m)}for(let d of e.nodes)a.has(d.id)||a.set(d.id,d.id);let i=new Map(e.nodes.map(d=>[d.id,d])),u=[];for(let d of e.nodes){if(a.get(d.id)!==d.id)continue;let h=l.get(d.id);if(!h){u.push(d);continue}let p=h.map(g=>i.get(g)),m=p.filter(_),b=m.length?m.reduce((g,M)=>g.stable.start<=M.stable.start?g:M):p[0];u.push({id:d.id,name:d.name,length:p.reduce((g,M)=>g+M.length,0),depth:d.depth,stable:b.stable})}let f=new Set,c=[];return e.edges.forEach((d,h)=>{let p=a.get(d.from),m=a.get(d.to),b=n[h],g=`${p}|${b.from}>${m}|${b.to}`;s.has(h)||f.has(g)||(f.add(g),c.push({from:p,to:m,...b.from==="start"?{fromStrand:It(Le(p))}:{},...b.to==="end"?{toStrand:It(Le(m))}:{}}))}),{graph:{...e,nodes:u,edges:c,paths:void 0,pathVisits:void 0},runs:l}}function ml(e,n,t){let r={...e};for(let[o,s]of n){let a=e[o];if(!a?.length)continue;delete r[o];let l=An(a),i=l.at(-1),u=s.map(d=>Math.max(t(d),0)),f=u.reduce((d,h)=>d+h,0)||s.length,c=0;s.forEach((d,h)=>{let p=f?u[h]/f*i:i/s.length,m=h===s.length-1?i:c+p;r[d]=vt(a,l,c,m),c=m})}return r}function bl(e,n,t){if(e.length===0)return[];let r=An(e),o=r.at(-1);return vt(e,r,n*o,t*o)}function _t(e){let n=An(e);return En(e,n,n.at(-1)/2)}function An(e){let n=[0];for(let t=1;t<e.length;t++)n.push(n[t-1]+Math.hypot(e[t].x-e[t-1].x,e[t].y-e[t-1].y));return n}function En(e,n,t){if(t<=0)return{...e[0]};let r=1;for(;r<n.length&&n[r]<t;)r++;if(r>=n.length)return{...e.at(-1)};let o=n[r]-n[r-1],s=o>0?(t-n[r-1])/o:0,a=e[r-1],l=e[r];return{x:a.x+(l.x-a.x)*s,y:a.y+(l.y-a.y)*s}}function vt(e,n,t,r){let o=[En(e,n,t)];for(let s=0;s<e.length;s++)n[s]>t&&n[s]<r&&o.push({...e[s]});return o.push(En(e,n,r)),o}function qr(e){return e.length*6.2+8}function Fr(e,n){return e.x0<n.x1&&e.x1>n.x0&&e.y0<n.y1&&e.y1>n.y0}function Bt(e,n){return e.x1<=0||e.x0>=n.width||e.y1<=0||e.y0>=n.height}function Wr(e,n){return e.x0>=0&&e.x1<=n.width&&e.y0>=0&&e.y1<=n.height}function Tn(e,n=[]){let t=[...n];return r=>Bt(r,e)||t.some(o=>Fr(o,r))?!1:(t.push(r),!0)}function en(e,n,t){let r=[];for(let s of e){let a=qr(s.text),l=s.y>n.height/2?-1:1;for(let i=0;i<=(s.stack??0);i++){let u=s.y+l*i*19,f={x0:s.x-a/2,x1:s.x+a/2,y0:u-11-4,y1:u+4};if(Bt(f,n))break;if(Wr(f,n)&&t(f)){r.push({...s,y:u,w:a});break}}}return r}var Hr=26,Xr=4e4,kr=8,Yr=8,zr=5.7,jr=6,Vr=13,$e=3,nn=2;function rn(e){return(e.length*zr+jr+$e)/2}var le=(Vr+$e)/2,Ur=rn("1 bp");function tn(e,n,t){return{x0:n-e,x1:n+e,y0:t-le,y1:t+le}}function Dt(e,n){return{x0:6,x1:6+e.length*6.2+8+$e,y0:n-8-$e/2,y1:n+8+$e/2}}function on(e,n,t){return t<n?e:Math.min(Math.max(e,n),t)}function sn(e){return e>=1e6?`${+(e/1e6).toFixed(1)} Mb`:e>=1e3?`${+(e/1e3).toFixed(1)} kb`:`${e} bp`}function Ot(e){return`${sn(e)} del`}function Kr(e,n,t,r){let o=e[0],s=e[e.length-1],a=(n.x-(o.x0+s.x1)/2)*t,l=(n.y-(o.y0+s.y1)/2)*r,i=Math.hypot(a,l);return i===0?{x:0,y:-1}:{x:a/i,y:l/i}}function Qr(e,n){return Math.abs(n.x)*(e.x1-e.x0)/2+Math.abs(n.y)*(e.y1-e.y0)/2}function Zr(e,n){let t=(e.x1-e.x0)/2,r=(e.y1-e.y0)/2;return Math.min(n.x===0?1/0:t/Math.abs(n.x),n.y===0?1/0:r/Math.abs(n.y))}function Jr(e,n,t){let r=0;for(let o=1;o<e.length;o++){let s=e[o-1],a=e[o];r+=Math.hypot((a.x-s.x)*n,(a.y-s.y)*t)}return r}function es(e,n){let t=1/0,r=1/0,o=-1/0,s=-1/0;for(let a of e)for(let l of n[a]??[])t=Math.min(t,l.x),o=Math.max(o,l.x),r=Math.min(r,l.y),s=Math.max(s,l.y);return t===1/0?void 0:{minX:t,minY:r,maxX:o,maxY:s}}var Pt=new WeakMap;function ns(e,n){let t=Pt.get(e);if(t?.nodeLengths===n)return t.order;let r=Object.entries(e).filter(([o,s])=>s.length>0&&n.has(o)).sort(([o],[s])=>n.get(s)-n.get(o));return Pt.set(e,{nodeLengths:n,order:r}),r}var Gt=new WeakMap;function ts(e,n,t,r){let{scaleX:o,scaleY:s}=t,a=Gt.get(e);if(a?.deletions===n&&a.scaleX===o&&a.scaleY===s&&a.version===r)return a.arcs;let l=[];for(let i of[...n].sort((u,f)=>f.bp-u.bp)){let u=Ct(e,i,t),f=u?Nt(u):void 0;if(u&&f){let{minX:c,minY:d,maxX:h,maxY:p}=Lt(u);l.push({deletion:i,curves:u,apex:f,extent:Math.max((h-c)*o,(p-d)*s)})}}return Gt.set(e,{deletions:n,scaleX:o,scaleY:s,version:r,arcs:l}),l}function $t({nodePositions:e,nodeLengths:n,deletions:t,alleleDeletions:r,axis:o,translateX:s,translateY:a,width:l,height:i,version:u=0}){let{scaleX:f,scaleY:c}=o,d=[],h=[];for(let{deletion:b,curves:g,apex:M,extent:R}of ts(e,t,o,u))if(R>=Hr){let N=Ot(b.bp),S=rn(N),x=tn(S,0,0),L=M.x*f+s,A=M.y*c+a,C=R>=S*2/nn,P=Kr(g,M,f,c),D=Qr(x,P)+Yr,$=C?L:on(L+P.x*D,S,l-S),W=C?A:on(A+P.y*D,le,i-le),H=Math.hypot($-L,W-A),U;if(H>0){let j={x:($-L)/H,y:(W-A)/H},V=Zr(x,j);U={arcX:L,arcY:A,labelX:$-j.x*V,labelY:W-j.y*V}}d.push({label:{key:`del:${b.edgeIndex}`,text:N,x:$,y:W,kind:"deletion",leader:U},box:tn(S,$,W)})}let p=new Set,m=[...r??[]].sort((b,g)=>g.bp-b.bp);for(let b of m){let g=es(b.nodeIds,e);if(!g)continue;let M=Math.max(g.minX*f+s,0),R=Math.min(g.maxX*f+s,l),N=Math.max(g.minY*c+a,0),S=Math.min(g.maxY*c+a,i);if(R<M||S<N)continue;let x=Math.max(R-M,S-N),L=Ot(b.bp),A=rn(L);if(x>=A*2/nn){for(let D of b.nodeIds)p.add(D);let C=on((M+R)/2,A,l-A),P=on((N+S)/2,le,i-le);d.push({label:{key:`alleledel:${b.nodeIds.join(",")}`,text:L,x:C,y:P,kind:"deletion"},box:tn(A,C,P)})}}for(let[b,g]of ns(e,n)){if(p.has(b))continue;let{x:M,y:R}=_t(g),N=R*c+a;if(N+le<=0||N-le>=i)continue;let S=Jr(g,f,c);if(S<Ur*2/nn)continue;let x=sn(n.get(b)),L=rn(x);if(S>=L*2/nn){let A=M*f+s;h.push({label:{key:`node:${b}`,text:x,x:A,y:N,kind:"node"},box:tn(L,A,N)})}}return{deletions:d,nodes:h}}function qt(e,n){return Math.max(kr,Math.floor(e*n/Xr))}function Cn(e,n,t=1/0){let r=[];for(let{label:o,box:s}of e){if(r.length>=t)break;n(s)&&r.push(o)}return r}var as=3.4,Ft=6,is=8,ls=18,cs=2,us=15;function Wt(e,n){return[...e].sort((t,r)=>n(r)-n(t))}function El(e){let{paneWidth:n,canvasHeight:t,translateX:r,translateY:o}=e,{scaleX:s,scaleY:a}=e.axisScale,l={width:n,height:t},i=x=>({x:x.x*s+r,y:x.y*a+o}),u=e.drawnRowLabels.map(({label:x,y:L})=>Dt(x,L*a+o)),f=e.referenceStripZonePx??0;f>0&&u.push({x0:0,x1:n,y0:0,y1:f}),e.legendSize.width>0&&u.push({x0:n-Ft-e.legendSize.width,x1:n,y0:0,y1:f+Ft+e.legendSize.height}),e.poppedFrom&&u.push({x0:0,x1:60+e.poppedFrom.label.length*6.2*1.2,y0:0,y1:44});let c=Tn(l,u),d=e.drawnRowLabels.length>0,h=d?Math.min(...e.drawnRowLabels.map(x=>x.y)):1/0,p=e.contigThickness*as,m=Wt(e.bubbleHalos,x=>x.members),b=en(m.map(x=>{let{x:L,y:A}=i({x:x.labelAt.x,y:Math.min(x.labelAt.y,h)}),C=A-p/2-6;return{item:x,x:L,y:d?Math.max(C,us):C,text:x.label}}),l,c),g=en(Wt(e.genePins,x=>x.gene.end-x.gene.start).map(x=>{let{x:L,y:A}=i(x.at);return{item:x,x:L,y:A+e.contigThickness+ls,text:x.covered<.98?`${x.gene.name} \u2026`:x.gene.name,stack:cs}}),l,c),M=e.nodePositions&&e.labelsNodeSizes?$t({nodePositions:e.nodePositions,nodeLengths:e.nodeLengths,deletions:e.showDeletionEdges?e.deletions:[],alleleDeletions:e.alleleDeletions,axis:e.axisScale,translateX:r,translateY:o,width:n,height:t,version:e.positionsVersion}):{deletions:[],nodes:[]},R=Cn(M.deletions,c),N=en(m.flatMap(x=>x.routes.map(L=>{let{x:A,y:C}=i(L.at);return{item:{halo:x,route:L},x:A,y:C+4,text:L.text,stack:is}})),l,c),S=Cn(M.nodes,c,qt(n,t));return{bubbles:b,genes:g,routes:N,sizes:[...R,...S]}}function Il(e){return Ge(e)?Ht(e):void 0}function Ht(e){return{sample:Ke(e.stable.refName),haplotype:st(e.stable.refName),refName:at(e.stable.refName),start:e.stable.start,end:e.stable.start+e.length}}function ds(e,n){let t=[...e].sort((s,a)=>s.start-a.start),r,o;for(let s of t)o=o&&s.start-o.end<=n?{start:o.start,end:Math.max(o.end,s.end)}:{start:s.start,end:s.end},(!r||o.end-o.start>r.end-r.start)&&(r=o);return r}function fs(e){return`${e.haplotype??e.sample}\0${e.refName}`}function _l(e,{maxGap:n=1/0}={}){let t=new Map;for(let o of e.nodes)if(Ge(o)){let s=Ht(o),a=fs(s),l=t.get(a);l?(l.spans.push({start:s.start,end:s.end}),l.rank=Math.min(l.rank,o.stable.rank)):t.set(a,{sample:s.sample,haplotype:s.haplotype,refName:s.refName,rank:o.stable.rank,spans:[{start:s.start,end:s.end}]})}let r=new Map;for(let o of t.values()){let s=ds(o.spans,n);if(s){let a={sample:o.sample,haplotype:o.haplotype,refName:o.refName,start:s.start,end:s.end,rank:o.rank,nodeCount:o.spans.length},l=o.haplotype??o.sample,i=r.get(l);(!i||a.nodeCount>i.nodeCount)&&r.set(l,a)}}return[...r.values()].sort((o,s)=>+(o.rank!==0)-+(s.rank!==0)||o.sample.localeCompare(s.sample)||(o.haplotype??"").localeCompare(s.haplotype??""))}function ps(e,n){return(n.haplotype===void 0?void 0:e(n.haplotype))??e(n.sample)}function vl(e,n){return e.flatMap(t=>{let r=ps(n,t);return r===void 0?[]:[{...t,sample:r}]})}function Bl(e){return`${e.refName}:${e.start+1}-${e.end}`}function In(e){return`${e.refName}:${(e.start+1).toLocaleString()}-${e.end.toLocaleString()}`}var Xt=ie(Ue(),1),kt=(0,te.createSvgIcon)((0,Xt.jsx)("path",{d:"M9.01 14H2v2h7.01v3L13 15l-3.99-4zm5.98-1v-3H22V8h-7.01V5L11 9z"}),"CompareArrows");var Yt=ie(Ue(),1),zt=(0,te.createSvgIcon)((0,Yt.jsx)("path",{d:"M17 5h-2V3h2zm-2 16h2v-2.59L19.59 21 21 19.59 18.41 17H21v-2h-6zm4-12h2V7h-2zm0 4h2v-2h-2zm-8 8h2v-2h-2zM7 5h2V3H7zM3 17h2v-2H3zm2 4v-2H3c0 1.1.9 2 2 2M19 3v2h2c0-1.1-.9-2-2-2m-8 2h2V3h-2zM3 9h2V7H3zm4 12h2v-2H7zm-4-8h2v-2H3zm0-8h2V3c-1.1 0-2 .9-2 2"}),"HighlightAlt");var jt=ie(Ue(),1),_n=(0,te.createSvgIcon)((0,jt.jsx)("path",{d:"M3 16h5v-2H3zm6.5 0h5v-2h-5zm6.5 0h5v-2h-5zM3 20h2v-2H3zm4 0h2v-2H7zm4 0h2v-2h-2zm4 0h2v-2h-2zm4 0h2v-2h-2zM3 12h8v-2H3zm10 0h8v-2h-8zM3 4v4h18V4z"}),"LineStyle");var hs="Linear genome view",ms="Linear synteny view",bs="No synteny track in this session aligns these assemblies";function Vt({label:e,icon:n,entries:t,entryLabel:r,onSelect:o,disabledHelpText:s,nameTheOnlyEntry:a=!0}){let[l]=t,i=s?[{label:e,icon:n,disabled:!0,disabledHelpText:s,onClick:()=>{}}]:[];return t.length===0?i:t.length===1&&l?[{label:a?`${e} \u2014 ${r(l)}`:e,icon:n,onClick:o(l)}]:[{label:e,icon:n,type:"subMenu",subMenu:t.map(u=>({label:r(u),onClick:o(u)}))}]}function Yl({contributors:e,syntenyTracks:n,onShowLinear:t,onShowSynteny:r}){let o=n[0]?.coverage??e.length;return[...Vt({label:hs,icon:_n,entries:e,entryLabel:s=>`${s.sample} ${In(s)}`,onSelect:s=>()=>{t({location:s,assembly:s.sample})}}),...e.length>=2?Vt({label:`${ms} (${o} assemblies)`,icon:kt,entries:n,entryLabel:s=>s.coverage<o?`${s.name} (${s.coverage} of ${o})`:s.name,nameTheOnlyEntry:!1,disabledHelpText:bs,onSelect:s=>()=>{r(s.trackId)}}):[]]}function zl({own:e,reference:n,highlight:t,onShowLinear:r,onHighlight:o}){let s=[e?{label:`Open in ${e.assembly} \u2014 ${In(e.location)}`,target:e}:void 0,n&&n.assembly!==e?.assembly?{label:`Open in ${n.assembly} \u2014 around this node`,target:n}:void 0].filter(a=>a!==void 0);return[...o&&t?[{label:`Highlight in ${t.assembly}`,icon:zt,onClick:()=>{o(t)}}]:[],...s.map(({label:a,target:l})=>({label:a,icon:_n,onClick:()=>{r(l)}}))]}var ys=10;function xs(e,n){let t=e.end-e.start,r=e.longestAlleleLength-e.shortestAlleleLength;return n.find(o=>r<=ys*(o.end-o.start)&&(t===0?o.start<=e.start&&e.start<=o.end:Math.min(e.end,o.end)-Math.max(e.start,o.start)>=t/2))}var Ms=40;function he(e){return e>=1e4?`${(e/1e3).toFixed(0)} kb`:e>=1e3?`${(e/1e3).toFixed(1)} kb`:`${e} bp`}function ws(e,n){let t=he(e),r=he(n),o=r.slice(r.indexOf(" "));return t.endsWith(o)?`${t.slice(0,-o.length)}\u2013${r}`:`${t}\u2013${r}`}function Kt(e){return e>=ct?"\u22652.1B routes":e>=1e5?`${(e/1e3).toFixed(0)}k routes`:`${e} routes`}function Kl(e,n=[]){let t=Ss(e,n);return e.partial?{...t,label:`${t.label}, partial`}:t}function Ss(e,n){let t=e.end-e.start,{shortestAlleleLength:r,longestAlleleLength:o,pathCount:s,segmentCount:a,inversion:l}=e,i=ws(r,o);if(a>=Ms)return{kind:"superbubble",label:`${i} superbubble, ${a} segments, ${Kt(s)}${l?", inv":""}`};if(l)return{kind:"inversion",label:`${he(o)} inv`};let u=o!==r?xs(e,n):void 0;if(u){let c=u.name?` (${u.name})`:"";return{kind:"repeat",label:`${i} repeat array${c}, ${Kt(s)}`}}let f=s>2?`, ${s} alleles`:"";return r===o&&r===t?t<=1?{kind:"snp",label:"SNP"}:{kind:"substitution",label:`${he(t)} sub`}:t===0||r===t?{kind:"insertion",label:`\u2264${he(o-t)} ins${f}`}:o===t?{kind:"deletion",label:`${he(t-r)} del${f}`}:{kind:"complex",label:`${he(t)} ref \u2192 ${i}, ${s} alleles`}}var Ql={snp:"#5b6b7a",substitution:"#5b6b7a",insertion:"#2f8fd6",deletion:"#c94040",inversion:"#e07b00",repeat:"#8e3fbf",superbubble:"#8e3fbf",complex:"#8e3fbf"},Zl={snp:"SNP",substitution:"substitution",insertion:"insertion",deletion:"deletion",inversion:"inversion",repeat:"repeat array",superbubble:"superbubble",complex:"site"};function Jl(e){return e.segments.split(",").filter(Boolean)}function vn(e){return e.id.endsWith("-")?"-":"+"}function Bn(e,n){return n===vn(e)?e.id:`-${e.id}`}function Ns(e,n){let t=e.map(()=>{}),r=0,o=0;for(;o<e.length&&e[o]!==void 0;o++)t[o]=r,r+=e[o];let s=n;for(let a=e.length-1;a>o&&e[a]!==void 0;a--)s-=e[a],t[a]=s;return t}function Ls(e,n){let t=e.path.map(()=>[]);if(e.cs===void 0)return t;let r=Ns(n,e.pathLength),o=(l,i=!1)=>{for(let u=0;u<r.length;u++){let f=r[u],c=n[u];if(f!==void 0&&c!==void 0&&f<=l&&(l<f+c||i&&l===f+c))return{s:u,local:l-f,room:f+c-l}}},s=(l,i)=>{let u=t[l],f=u.at(-1);i.type==="substitution"&&f?.type==="substitution"&&f.pos+f.seq.length===i.pos?f.seq+=i.seq:u.push(i)},a=e.pathStart;for(let l of lt(e.cs))if(l.op==="match")a+=l.length;else if(l.op==="substitution"){let i=o(a);i&&s(i.s,{type:"substitution",pos:i.local,seq:l.seq}),a+=1}else if(l.op==="insertion"){let i=o(a,a===e.pathEnd);i&&s(i.s,{type:"insertion",pos:i.local,seq:l.seq})}else{let i=l.length;for(;i>0;){let u=o(a),f=u?Math.min(i,u.room):i;u&&s(u.s,{type:"deletion",pos:u.local,length:f}),a+=f,i-=f}}return t}function Qt(e,n,t){let r=new Map(e.nodes.map(s=>[s.name,s])),o=[];for(let s of n){let a=s.path.map(f=>r.get(f.name)),l=a.flatMap((f,c)=>f?[c]:[]),i=l[0],u=l.at(-1);if(i!==void 0&&u!==void 0){let f=Ls(s,a.map(h=>h?.length)),c=l.map(h=>Bn(a[h],s.path[h].strand)),d=a[u].length;o.push({id:t+o.length,name:s.name,type:"read",sourceTrackID:1,sequence:c,sequenceNew:l.map((h,p)=>({nodeName:c[p],mismatches:f[h]})),firstNodeOffset:i===0?s.pathStart:0,finalNodeCoverLength:u===s.path.length-1?d-(s.pathLength-s.pathEnd):d,mapping_quality:s.mappingQuality===255?0:s.mappingQuality,is_secondary:s.secondary,is_reverse:s.strand==="-"})}}return o}var Rs=["#6baed6","#4292c6","#2171b5","#08519c"],Es=["#fb6a4a","#ef3b2c","#cb181d","#a50f15"];function Zt(e){let n=e.is_reverse?Es:Rs;return n[e.id%n.length]}var As=4;function qe(e){let n=new Map;function t(r,o){let s=n.get(r);s?s.push(o):n.set(r,[o])}for(let r of e.edges)t(r.from,r.to),t(r.to,r.from);return n}function Ts({nodeId:e,nodeById:n,neighbors:t}){let r=new Set([e]),o=[e],s=1/0,a=-1/0;for(let l=0;l<As&&o.length>0;l++){let i=[];for(let u of o){let f=t.get(u);if(f){for(let c of f)if(!r.has(c)){r.add(c);let d=n.get(c);d&&_(d)?(s=Math.min(s,d.stable.start+d.length),a=Math.max(a,d.stable.start)):i.push(c)}}}o=i}return s===1/0?void 0:{start:Math.min(s,a),end:Math.max(s,a)}}function Cs({nodeId:e,nodeById:n,neighbors:t}){let r=n.get(e);return r===void 0?void 0:_(r)?{start:r.stable.start,end:r.stable.start+r.length}:Ts({nodeId:e,nodeById:n,neighbors:t})}function Jt(e){let n=new Map(e.nodes.map(o=>[o.id,o])),t=qe(e),r=new Map;for(let o of e.nodes){let s=Cs({nodeId:o.id,nodeById:n,neighbors:t});s&&r.set(o.id,(s.start+s.end)/2)}return r}function an(e){let n=e.split("#"),t=n.length>=2?Number(n[1]):Number.NaN;return{sample:n[0]??e,haplotype:Number.isFinite(t)?t:void 0}}function Is(e,n,t,r){let o=[],s=[],a=[e];for(r.add(e);a.length>0;){let l=a.pop(),i=n.get(l);if(i){o.push(i);for(let u of t.get(l)??[]){let f=n.get(u);f&&(f.stable.rank===0?s.push(f):r.has(u)||(r.add(u),a.push(u)))}}}return{members:o,anchors:s}}function _s(e){if(new Set(e.map(l=>l.stable.refName)).size!==1||e.length<2)return;let t=[...e].sort((l,i)=>l.stable.start-i.stable.start),r=t[0],o=t.at(-1),s=r.stable.start+r.length,a=o.stable.start;return a>=s?{start:s,end:a,entry:r}:void 0}function vs(e,n,t){let r=new Map(e.map(a=>[a.id,a])),o=new Map,s=[];for(let a of t.get(n.id)??[])r.has(a)&&(o.set(a,0),s.push(a));if(s.length===0)for(let a of e)o.set(a.id,0),s.push(a.id);for(let a of s){let l=o.get(a)+r.get(a).length;for(let i of t.get(a)??[])if(r.has(i)){let u=o.get(i);(u===void 0||l<u)&&(o.set(i,l),s.push(i))}}return o}function Bs(e,n,t){let r=vs(e,n,t),o=[...e].sort((i,u)=>r.get(i.id)-r.get(u.id)),s=new Map(o.map((i,u)=>[i.id,u])),a=new Map(e.map(i=>[i.id,i])),l=new Map;for(let[i,u]of o.entries()){let f=0;for(let c of t.get(u.id)??[]){let d=a.get(c);d&&s.get(c)<i&&(f=Math.max(f,l.get(c)+d.length))}l.set(u.id,f)}return{offsets:l,ordered:o}}function no(e){let n=new Map(e.nodes.filter(Ge).map(s=>[s.id,s])),t=qe(e),r=new Set,o=[];for(let s of n.values())if(s.stable.rank!==0&&!r.has(s.id)){let{members:a,anchors:l}=Is(s.id,n,t,r),i=_s(l);if(i){let{offsets:u,ordered:f}=Bs(a,i.entry,t);o.push({start:i.start,refSpan:i.end-i.start,nodeIds:f.map(c=>c.id),nodeOffsets:f.map(c=>u.get(c.id)),pathLength:f.reduce((c,d)=>Math.max(c,u.get(d.id)+d.length),0)})}}return o.sort((s,a)=>s.start-a.start)}var Os=1,Ps=.015;function ln({graph:e,span:n,rowY:t,positions:r}){let o=n*Ps,s=no(e),a=new Map(e.nodes.map(f=>[f.id,f])),l=[];for(let f of s){let c=f.refSpan-f.pathLength;c>=Os&&l.push({nodeIds:f.nodeIds,bp:c});let d=Math.max(f.refSpan,o),h=f.pathLength>0?d/f.pathLength:0;f.nodeIds.forEach((p,m)=>{let b=a.get(p);if(b){let g=f.nodeOffsets[m],M=d/f.nodeIds.length,R=S=>f.start+S*h,N=t(b);r[p]=h?[{x:R(g),y:N},{x:R(g+b.length),y:N}]:[{x:f.start+m*M,y:N},{x:f.start+(m+1)*M,y:N}]}})}let i=qe(e),u=Object.keys(r);for(let f of u){let c=r[f].at(-1).x;for(let d of i.get(f)??[]){let h=a.get(d);if(h&&!(d in r)){let p=t(h);r[d]=[{x:c,y:p},{x:c+o,y:p}],u.push(d)}}}return l}function Gs(e){let n=new Set;for(let t of e.nodes)t.stable&&n.add(t.stable.rank);return new Map([...n].sort((t,r)=>t-r).map((t,r)=>[t,r]))}function to(e,n){let t=we(e);if(t.length===0)return;let r=Ze(t,n),o=Gs(e),s=Se(t),a=ln({graph:e,span:r,rowY:i=>(i.stable?o.get(i.stable.rank)??o.size:o.size)*20,positions:s}),l=[...o].map(([i,u])=>({label:i===0?"Reference (rank 0)":`Rank ${i}`,y:u*20}));return{nodePositions:s,rowLabels:l,referenceAxis:!0,pixelRows:!0,alleleDeletions:a}}var Ds=12,On=20;function oo(e){return 10+8*Math.log2(1+e.length)}function ce(e,n){return e<n?-1:e>n?1:0}function so(e,n){return ce(e.at,n.at)||ce(e.side,n.side)||ce(e.hops,n.hops)}var ro={at:1/0,side:0,hops:0};function $s(e){let n=new Map(e.nodes.map(r=>[r.id,r])),t=Ws(e);return e.nodes.map(r=>r.id).sort((r,o)=>so(t.get(r)??ro,t.get(o)??ro)||ce(n.get(r).name,n.get(o).name)||ce(r,o))}function qs(e,n){let t=new Set(e.nodes.map(f=>f.id)),r=new Map(n.map(f=>[f.id,f.stable.start])),o=new Map,s=new Map;n.forEach((f,c)=>{s.set(f.id,n[c-1]?.id),o.set(f.id,n[c+1]?.id)});let a=n[0]?.id,l=(f,c)=>{let d=o.get(f);o.set(f,c),s.set(c,f),o.set(c,d),d!==void 0&&s.set(d,c)},i=(f,c)=>{let d=s.get(f);d===void 0?(a=c,s.set(c,void 0),o.set(c,f),s.set(f,c)):l(d,c)};for(let f of e.paths??[]){let c=0,d=0,h;for(let g of f.nodeIds){let M=r.get(g);M!==void 0&&(h!==void 0&&M!==h&&(M>h?c++:d++),h=M)}let p=d>c?[...f.nodeIds].reverse():f.nodeIds,m,b=new Set;for(let g of p)if(o.has(g)){if(m===void 0)for(let M of b)i(g,M);m=g}else t.has(g)&&(m===void 0?b.add(g):(l(m,g),m=g))}let u=[];for(let f=a;f!==void 0;f=o.get(f))u.push(f);return u}function Fs(e,n){let t=[...n.entries()].sort(([a,l],[i,u])=>so(l,u)||ce(a,i)).map(([a])=>a),r=t.length;for(let a of t){let l=n.get(a);for(let i of e.get(a)??[])n.has(i)||(n.set(i,{at:l.at,side:1,hops:l.hops+1}),t.push(i))}let o=new Set(t.slice(r)),s=new Set;for(let a of o){if(s.has(a))continue;let{at:l}=n.get(a),i=[a];s.add(a);let u=!0;for(let f of i){let{hops:c}=n.get(f),d=!0,h=!1;for(let p of e.get(f)??[]){let m=n.get(p);m.at===l&&o.has(p)?(d&&=m.hops<=c,s.has(p)||(s.add(p),i.push(p))):m.at!==l&&(h||=m.at<l,u&&=m.at<l)}u&&=!d||h}if(u)for(let f of i){let c=n.get(f);n.set(f,{...c,side:-1,hops:-c.hops})}}}function Ws(e){let n=new Set(e.nodes.map(s=>s.id)),t=new Map;for(let s of e.edges)!n.has(s.from)||!n.has(s.to)||((t.get(s.from)??t.set(s.from,[]).get(s.from)).push(s.to),(t.get(s.to)??t.set(s.to,[]).get(s.to)).push(s.from));let r=e.nodes.filter(_).sort((s,a)=>ce(s.stable.start,a.stable.start)||ce(s.name,a.name)),o=new Map;return qs(e,r).forEach((s,a)=>{o.set(s,{at:a,side:0,hops:0})}),Fs(t,o),o}function Hs(e){let n=new Map(e.nodes.map(c=>[c.id,c])),t=$s(e),r=new Map(t.map((c,d)=>[c,d])),o=c=>_(n.get(c)),s=new Map(t.map(c=>[c,[]])),a=new Map(t.map(c=>[c,[]])),l=(c,d)=>{s.get(d).push(c),a.get(c).push(d)};for(let c of e.edges)c.from===c.to||!r.has(c.from)||!r.has(c.to)||(r.get(c.from)<r.get(c.to)?l(c.from,c.to):l(c.to,c.from));let i=t.filter(o);for(let c=1;c<i.length;c++)l(i[c-1],i[c]);let u=new Map;for(let c of t){let d=0;for(let h of s.get(c))d=Math.max(d,u.get(h)+1);u.set(c,d)}for(let c=t.length-1;c>=0;c--){let d=t[c],h=a.get(d);if(o(d)||h.length===0)continue;let p=1/0;for(let g of h)p=Math.min(p,u.get(g));let m=u.get(d),b=Math.max(m,p-1);u.set(d,Math.floor((m+b)/2))}let f=[];for(let c of t)(f[u.get(c)]??=[]).push(c);return{layers:f,layerOf:u,preds:s}}function ao(e){if(!e.nodes.some(_))return;let n=new Map(e.nodes.map(c=>[c.id,c])),t=c=>_(n.get(c)),{layers:r,layerOf:o,preds:s}=Hs(e),a=r.map(c=>{let d=0;for(let h of c)d=Math.max(d,oo(n.get(h)));return d}),l=[],i=0;for(let c=0;c<r.length;c++)l[c]=i,i+=a[c]+Ds;let u=new Map;for(let c of r){let d=[];for(let p of c){if(t(p)){u.set(p,0);continue}let m=s.get(p).filter(g=>u.has(g)),b=0;for(let g of m)b+=u.get(g)/m.length;d.push({id:p,ideal:b===0?On:b})}d.sort((p,m)=>p.ideal-m.ideal);let h=new Set([0]);for(let{id:p,ideal:m}of d){let b=Math.round(m/On);b===0&&(b=m>=0?1:-1);for(let g=1;h.has(b);g++)h.has(b+g)?h.has(b-g)||(b-=g):b+=g;h.add(b),u.set(p,b*On)}}let f={};for(let c of e.nodes){let d=o.get(c.id),h=oo(c),p=l[d]+a[d]/2,m=u.get(c.id);f[c.id]=[{x:p-h/2,y:m},{x:p+h/2,y:m}]}return{nodePositions:f,referenceAxis:!1,pixelRows:!0}}function Xs(e,n){let t=new Map;for(let o of e.nodes)if(Qe(o)){let{sample:s}=an(o.stable.refName);t.set(s,(t.get(s)??0)+o.length)}let r=new Map(n.map((o,s)=>[o,s]));return[...t.keys()].sort((o,s)=>(r.get(o)??1/0)-(r.get(s)??1/0)||t.get(s)-t.get(o)||o.localeCompare(s))}function io(e,n,t=[]){let r=we(e),o=Xs(e,t);if(r.length===0||o.length===0)return;let s=Ze(r,n),a=new Map(o.map((d,h)=>[d,h+1])),l=Se(r),i=o.length+1,u=ln({graph:e,span:s,rowY:d=>{let h=d.stable;return((h?a.get(an(h.refName).sample):void 0)??i)*20},positions:l}),c=[{label:an(r[0].stable.refName).sample,y:0},...o.map(d=>({label:d,y:a.get(d)*20}))];return{nodePositions:l,rowLabels:c,sampleRows:o,referenceAxis:!0,pixelRows:!0,alleleDeletions:u}}function Pn(){return{rectangles:[],curves:[],corners:[],verticalRectangles:[],featureRectangles:[]}}var He=!1,lo=["#1f77b4","#ff7f0e","#2ca02c","#d62728","#9467bd","#8c564b","#e377c2","#7f7f7f"];function Ro(e){return{mergeNodesFlag:e.mergeNodes??!0,showReads:e.showReads??!0,coarsenedReadView:e.coarsenedReadView??!1,ignoreStrand:e.ignoreStrand??!1,nodeWidthOption:e.nodeWidthOption??"normal",charWidth:e.charWidth??8.401,trackWidth:e.trackWidth??15,mappingQualityCutoff:e.mappingQualityCutoff??0,focusReadNames:e.focusReadNames??null,showExonsFlag:e.showExons??!1,trackColor:e.trackColor??(n=>lo[n.id%lo.length]),trackAlpha:e.trackAlpha??(()=>1)}}var v=Ro({}),Fn=null,y=[],w=[],E=[],I=new Map,Ae=[],T=[],me=[],Hn=[],Xn=[],Q=-1,G=Pn(),Wn={minX:0,maxX:0,minY:0,maxY:0},Fe,cn=new Map;function Eo(e,n,t=[],r={}){if(v=Ro(r),Fn=r.bed??null,G=Pn(),me=[],Hn=[],Xn=[],Ae=[],Wn={minX:0,maxX:0,minY:0,maxY:0},Fe=void 0,cn=new Map,e.length===0||n.length===0)return;let o=Gn(e);y=[],o.forEach((i,u)=>{y[u+1]=i}),w=Gn(n);let s=t.length>0;E=Gn(Ia(t));for(let i=w.length-1;i>=0;i-=1){let u=w[i];if(u.type===void 0&&(u.type="haplotype"),u.hidden===!0){w.splice(i,1);continue}u.indexOfFirstBase!==void 0&&(Fe=u.name)}if(w.length===0)return;Ys(0),I=fo(),Re(w),v.showReads&&E.length>0&&Re(E),Lo(),v.mergeNodesFlag&&(po(),$n(),v.showReads&&E.length>0&&Dn(),Ca(),I=fo(),Lo(),Re(w),v.showReads&&E.length>0&&Re(E)),po(),fa(),He&&console.log(`${y.length} nodes.`),$n(),Q=mo(),pa(),$n(),Q=mo();let a=v.showReads&&v.coarsenedReadView&&!s;if(a){let i=Fe===void 0?-1:w.findIndex(c=>c.name===Fe),u=i===-1?0:i,f=w.filter((c,d)=>d!==u);if(f.length>0){let c=Mo(f,"haplotype");c.length>0&&(w=[w[u]],E=c)}}Na(),ga(),v.showExonsFlag&&Fn!==null&&Sa();let l=v.showReads&&v.coarsenedReadView&&E.length>0&&!a;return v.showReads&&E.length>0?(zs(),Dn(),Re(E),l&&(E=Mo(E,"read"),Dn(),Re(E)),Vs(),w=w.concat(E)):y.forEach(i=>{i&&(i.incomingReads=[],i.outgoingReads=[],i.internalReads=[])}),ma(),Ea(),He&&(console.log("Tracks:"),console.log(w),console.log("Nodes:"),console.log(y),console.log("Lane assignment:"),console.log(me)),aa(),{nodes:y,tracks:w,reads:E,nodeMap:I,shapes:G,bounds:Wn,maxOrder:Q,trackForRuler:Fe,coarsenedEdgeMeta:cn}}function Ee(e,n="plain"){return v.trackColor(e,n)}function ks(e){return v.trackAlpha(e)}function Gn(e){return structuredClone(e)}function k(e){let n=String(e);return n.length>=1&&n.startsWith("-")}function O(e){return k(e)?String(e).substring(1):e}function ue(e){return k(e)?e:`-${e}`}function Ao(e){return k(e)?O(e):ue(e)}function un(e){return(e||1/e)>=0}function Ys(e){let n=new Set,t=w[e].sequence;for(let r=0;r<t.length;r+=1){let o=t[r];if(k(o)){let s=O(o),a=t.indexOf(s);(a===-1||a>r)&&n.add(s)}}for(let r of w){t=r.sequence;for(let o=0;o<t.length;o+=1){let s=t[o];k(s)?n.has(O(s))&&(t[o]=O(s)):n.has(s)&&(t[o]=ue(s))}}y.forEach(r=>{r&&n.has(r.name)&&(r.seq=r.seq.split("").reverse().join(""))})}function zs(){Ae=[];for(let n=0;n<=Q;n+=1)Ae[n]=[];let e=new Map;y.forEach(n=>{n.y!==void 0&&co(e,n.order,n.y+n.contentHeight)}),w.forEach(n=>{n.type==="haplotype"&&n.path.forEach(t=>{co(e,t.order,(t.y??0)+n.width)})}),y.forEach((n,t)=>{n.order>=0&&n.y===void 0&&(n.y=(e.get(n.order)??0)+25,n.contentHeight=0,Ae[n.order].push(t))})}function co(e,n,t){let r=e.get(n);r===void 0?e.set(n,t):e.set(n,Math.max(r,t))}var dn=7;function js(){y.forEach(e=>{e&&(e.incomingReads=[],e.outgoingReads=[],e.internalReads=[])}),E.forEach((e,n)=>{if((e.width===void 0||e.width===0)&&(e.width=dn),e.path.length===1){let t=e.path[0].node;t!==null&&y[t].internalReads.push(n)}else e.path.forEach((t,r)=>{if(r===0){let o=e.path[0].node;o!==null&&y[o].outgoingReads.push([n,r])}else{let o=t.node;o!==null&&y[o].incomingReads.push([n,r])}})})}function Vs(){ta(),js();let e=y.slice();e.sort(fn);let n=new Map;for(let a=0;a<E.length;a++){let l=E[a].sourceTrackID,i=n.get(l);i===void 0?n.set(l,[a]):i.push(a)}let t=Array.from(n.keys()).sort((a,l)=>a-l),r=t.length>1?dn:0;e.forEach(a=>{if(a)for(let l of t)Us(n.get(l),a,r)});let o=na(),s=[];E.forEach((a,l)=>{let i=a.path.length,u=new Array(i),f;for(let h=0;h<i;h++){h===0?u[h]=void 0:u[h]=f??a.path[0];let p=a.path[h];p.node!==null&&p.node&&(f=p)}let c=new Array(i),d;for(let h=i-1;h>=0;h--){h===i-1?c[h]=void 0:c[h]=d??a.path[i-1];let p=a.path[h];p.node!==null&&p.node&&(d=p)}a.path.forEach((h,p)=>{if(h.y===void 0){let m=u[p],b=c[p],g=!!(b&&m&&b.isForward!==void 0&&m.isForward!==void 0&&(m.order>b.order||b.order===m.order&&b.isForward===m.isForward));h.betweenCycleReverseTraversal=g,s.push({readIndex:l,pathIndex:p,previousY:m?.y,previousNode:m?.node})}})}),s.sort(Qs),s.forEach(a=>{let l=E[a.readIndex],i=l.path[a.pathIndex];i.y=o[i.order],o[i.order]+=l.width}),He&&(console.log("Reads:"),console.log(E))}function Us(e,n,t){let r=new Set(e),o=n.incomingReads.filter(([p])=>r.has(p)),s=n.outgoingReads.filter(([p])=>r.has(p)),a=n.internalReads.filter(p=>r.has(p));o.length===0&&s.length===0&&a.length===0&&(t=0);let l=n.y+n.contentHeight+t*2,i=o.map(p=>{let[m,b]=p,g=E[m].path,M=b+1,R;for(let N=1;N<=b;N++){let S=g[b-N]?.y;if(S!==void 0){M=N,R=S;break}}return{entry:p,decisionStep:M,y:R}});i.sort(Ks),o=i.map(p=>p.entry);let u=l,f=new Map;o.forEach(p=>{let m=E[p[0]];m.path[p[1]].y=u,uo(f,m,p[1],u,n),u+=m.width});let c=u;s.sort(Zs);let d=new Map;u=l,s.forEach(p=>{let m=E[p[0]],b=m.firstNodeOffset??0;m.path[p[1]].y=u,d.set(u,b);let g=f.get(u);g===void 0||g+1<=b?(u+=m.width,c=Math.max(c,u)):(f.set(u,0),o.forEach(M=>{let R=E[M[0]],N=R.path[M[1]];N.y!==void 0&&N.y>=u&&(N.y+=m.width,uo(f,R,M[1],N.y,n))}),u+=m.width,c+=m.width)}),a.sort(ea),a.forEach(p=>{let m=E[p],b=m.firstNodeOffset??0,g=m.finalNodeCoverLength??0;for(u=l;b<(f.get(u)??-1/0)+2||g>(d.get(u)??1/0)-3;)u+=m.width||dn;m.path[0].y=u,f.set(u,g),c=Math.max(c,u)});let h=c-n.y-n.contentHeight;n.contentHeight+=h,wa(n,h)}function Ks(e,n){return e.decisionStep!==n.decisionStep?e.decisionStep-n.decisionStep:e.y===void 0||n.y===void 0?(e.y===void 0?0:1)-(n.y===void 0?0:1):e.y-n.y}function uo(e,n,t,r,o){t===n.path.length-1?e.set(r,n.finalNodeCoverLength??0):e.set(r,o.sequenceLength)}function Qs(e,n){let t=E[e.readIndex],r=E[n.readIndex],o=t.path[e.pathIndex],s=r.path[n.pathIndex];if(o.order!==s.order)return o.order-s.order;if(t.sourceTrackID!==r.sourceTrackID)return t.sourceTrackID-r.sourceTrackID;if(e.previousNode&&n.previousNode){let i=y[e.previousNode],u=y[n.previousNode];if(i&&u&&i.order!==u.order)return i.order-u.order}let a=e.previousY??0,l=n.previousY??0;return o.betweenCycleReverseTraversal&&s.betweenCycleReverseTraversal?l-a:a-l}function Zs([e,n],[t,r]){let o=E[e],s=E[t],a,l,i=o.path.length-1,u=s.path.length-1;for(;!a&&i>=0;)a=o.path[i]?.y,i-=1;for(;!l&&u>=0;)l=s.path[u]?.y,u-=1;if(a&&l)return a-l;let f=o.path[n]?.node,c=s.path[r]?.node,d=f!=null?y[f]:null,h=c!=null?y[c]:null;for(;d!=null&&h!=null&&d===h;){if(n<o.path.length-1){for(n+=1;o.path[n]?.node===null;)n+=1;let m=o.path[n]?.node;d=m!=null?y[m]:null}else d=null;if(r<s.path.length-1){for(r+=1;s.path[r]?.node===null;)r+=1;let m=s.path[r]?.node;h=m!=null?y[m]:null}else h=null}if(d!=null)return h!=null?fn(d,h):1;if(h!=null)return-1;let p=(o.firstNodeOffset??0)-(s.firstNodeOffset??0);return p!==0?p:(o.finalNodeCoverLength??0)-(s.finalNodeCoverLength??0)}function Js(e,n){if(e.path===void 0||n.path===void 0)return(e.path===void 0?0:1)-(n.path===void 0?0:1);let t=e.path,r=n.path,o=0,s=0;for(;o<t.length&&s<r.length;){let a=t[o],l=r[s],i=a.order,u=l.order;if(i===u&&a.node&&l.node)return(l.y??0)-(a.y??0);i<u?o+=1:(i>u||(o+=1),s+=1)}return 0}function ea(e,n){let t=E[e],r=E[n],o=t.firstNodeOffset??0,s=r.firstNodeOffset??0;if(o<s)return-1;if(o>s)return 1;let a=t.finalNodeCoverLength??0,l=r.finalNodeCoverLength??0;return a<l?-1:a>l?1:0}function na(){let e=[];for(let n=0;n<=Q;n+=1)e.push(0);return y.forEach(n=>{n.y!==void 0&&(e[n.order]=Math.max(e[n.order],n.y+n.contentHeight+20))}),w.forEach(n=>{n.path.forEach(t=>{e[t.order]=Math.max(e[t.order],(t.y??0)+n.width)})}),e}function ta(){E.forEach(e=>{let n=Math.abs(e.indexSequence[0]),t=un(e.indexSequence[0]),r=y[n],o,s;e.path=[],e.path.push({order:r.order,isForward:t,node:n});for(let a=1;a<e.sequence.length;a+=1)if(o=r,s=t,n=Math.abs(e.indexSequence[a]),t=un(e.indexSequence[a]),r=y[n],r.order>o.order){s||e.path.push({order:o.order,isForward:!0,node:null});for(let l=o.order+1;l<r.order;l+=1)e.path.push({order:l,isForward:!0,node:null});t?e.path.push({order:r.order,isForward:!0,node:n}):(e.path.push({order:r.order,isForward:!0,node:null}),e.path.push({order:r.order,isForward:!1,node:n}))}else if(r.order<o.order){s&&e.path.push({order:o.order,isForward:!1,node:null});for(let l=o.order-1;l>r.order;l-=1)e.path.push({order:l,isForward:!1,node:null});t?(e.path.push({order:r.order,isForward:!1,node:null}),e.path.push({order:r.order,isForward:!0,node:n})):e.path.push({order:r.order,isForward:!1,node:n})}else t!==s?e.path.push({order:r.order,isForward:t,node:n}):(e.path.push({order:r.order,isForward:!t,node:null}),e.path.push({order:r.order,isForward:t,node:n}))})}function Dn(){E.forEach(e=>{let n=0;for(;n<e.sequence.length&&e.sequence[n].startsWith("-");)n+=1;if(n===e.sequence.length){e.is_reverse=!0,e.sequence=e.sequence.reverse();for(let a=0;a<e.sequence.length;a+=1)e.sequence[a]=O(e.sequence[a]);let t=e.sequenceNew??[];t.reverse(),e.sequenceNew=t;for(let a of t)a.nodeName=O(a.nodeName),oa(a.mismatches,oe(a.nodeName).sequenceLength);let r=e.firstNodeOffset??0,o=oe(e.sequence[0]).sequenceLength;e.firstNodeOffset=o-(e.finalNodeCoverLength??0);let s=oe(e.sequence[e.sequence.length-1]).sequenceLength;e.finalNodeCoverLength=s-r}})}function oa(e,n){e.forEach(t=>{t.type==="insertion"?t.pos=n-t.pos:t.type==="deletion"?t.pos=n-t.pos-(t.length??0):t.type==="substitution"&&(t.pos=n-t.pos-(t.seq?.length??0)),t.seq!==void 0&&(t.seq=To(t.seq).split("").reverse().join(""))})}var ra={A:"T",T:"A",C:"G",G:"C"};function To(e){let n="";for(let t=e.length-1;t>=0;t-=1)n+=ra[e.charAt(t)]??"N";return n}function Re(e){e.forEach(n=>{n.indexSequence=[],n.sequence.forEach(t=>{let r=t,o=I.get(O(r)),s=y[o]?.switched??!1;s&&(r=Ao(r));let a=I.get(O(r));if(a===0)throw new Error("Node "+O(r)+" has prohibited index 0");k(r)!==s?n.indexSequence.push(-a):n.indexSequence.push(a)})})}var sa=20;function aa(){let e={minX:99,maxX:-99,minY:99,maxY:-99};y.forEach(n=>{n&&(n.x!==void 0&&(e.minX=Math.min(e.minX,n.x),e.maxX=Math.max(e.maxX,n.x+sa+n.pixelWidth)),n.y!==void 0&&(e.minY=Math.min(e.minY,n.y-10),e.maxY=Math.max(e.maxY,n.y+n.contentHeight+10)))}),w.forEach(n=>{n.path.forEach(t=>{let r=t.y??0;e.maxY=Math.max(e.maxY,r+n.width),e.minY=Math.min(e.minY,r)})}),Wn=e}function oe(e){let n=I.get(O(e)),t=n===void 0?void 0:y[n];if(t===void 0)throw new Error(`Unknown node ${e}`);return t}function fo(){return I=new Map,y.forEach((e,n)=>{e&&I.set(e.name,n)}),I}function po(){let e=y.map(()=>new Set),n=y.map(()=>new Set),t=r=>{for(let o=0;o<r.indexSequence.length-1;o+=1){let s=Math.abs(r.indexSequence[o]),a=Math.abs(r.indexSequence[o+1]);e[s].add(a),n[a].add(s)}};w.forEach(t),v.showReads&&E.length>0&&E.forEach(t),y.forEach((r,o)=>{r&&(r.successors=Array.from(e[o]),r.predecessors=Array.from(n[o]))})}function ho(e){let n=0,t=0,r=0;e.forEach(o=>{let s=Math.abs(o);if(o<0){T[s]===void 0&&(T[s]=t);let a=T[s];a<r&&(r=a),n=a,t=a-1}else{T[s]===void 0&&(T[s]=n);let a=T[s];n=a+1,t=a}}),r<0&&kn(-r)}function ia(e){let n=0,t=0;for(;n<e.length&&T[Math.abs(e[n])]===void 0;)n+=1;if(n>=e.length)return null;let r=e[n],o,s;r>=0?(o=T[r]-1,s=-1):(o=T[-r]+1,s=1);for(let a=n-1;a>=0;a-=1){let l=Math.abs(e[a]);T[l]===void 0&&(T[l]=o,t=Math.min(t,o),o+=s)}return t<0&&kn(-t),n}var la=-1;function ca(e){let n=[];for(let t of e)n.push(t??la);return n}function $n(){let e,n,t,r,o=0,s=v.showReads&&E.length>0?w.concat(E):w,a={stamp:new Int32Array(y.length),generation:0};T=new Array(y.length).fill(void 0),y.forEach(i=>{i&&(i.order=void 0)}),ho(w[0].indexSequence);for(let i=1;i<s.length;i+=1){if(He&&console.log(`generating order for track ${i+1}`),t=ia(s[i].indexSequence),t===null){s[i].type==="haplotype"?ho(s[i].indexSequence):(s.splice(i,1),E.splice(i-w.length,1),i-=1);continue}for(e=da(s[i].indexSequence);t<e.length;){for(r=t,t+=1;t<e.length&&T[e[t]]===void 0;)t+=1;if(t<e.length){n=T[e[r]]+1;for(let u=r+1;u<t;u+=1)T[e[u]]=n,n+=1;if(T[e[t]]>T[e[r]])T[e[t]]<n&&bo(e[t],e[t-1],n);else if(s[i].indexSequence[t]>=0&&!ua(e[t],e[r],a))bo(e[t],e[t-1],n);else if(s[i].indexSequence[r]<0||y[e[r+1]].degree<2&&T[e[t]]<T[e[r]]){n=T[e[r]]-1;for(let u=r+1;u<t;u+=1)T[e[u]]=n,o=Math.min(o,n),n-=1}}else if(s[i].indexSequence[r]>=0){n=T[e[r]]+1;for(let u=r+1;u<e.length;u+=1){let f=e[u];T[f]===void 0&&(T[f]=n,n+=1)}}else{n=T[e[r]]-1;for(let u=r+1;u<e.length;u+=1){let f=e[u];T[f]===void 0&&(T[f]=n,o=Math.min(o,n),n-=1)}}}}o<0&&kn(-o);let l=ca(T);T=l,l.forEach((i,u)=>{let f=y[u];f!==void 0&&(f.order=i)})}function ua(e,n,t){t.generation+=1;let{stamp:r,generation:o}=t,s=[e];for(r[e]=o;s.length>0;){let a=s.pop();if(a===n)return!0;for(let l of y[a].successors)r[l]!==o&&(r[l]=o,s.push(l))}return!1}function mo(){let e=-1;return T.forEach(n=>{n!==void 0&&n>e&&(e=n)}),e}function da(e){return e.map(n=>Math.abs(n))}function kn(e){T.forEach((n,t)=>{n!==void 0&&(T[t]=n+e)})}function bo(e,n,t){let r=new Map,o=[[e,t]],s=0;for(;s<o.length;){let[a,l]=o[s];s+=1;let i=T[a];i!==void 0&&i<l&&(!r.has(a)||r.get(a)<l)&&(r.set(a,l),y[a].successors.forEach(u=>{T[u]>i&&u!==n&&o.push([u,l+1])}),a!==e&&y[a].predecessors.forEach(u=>{T[u]>i&&u!==n&&o.push([u,l+1])}))}r.forEach((a,l)=>{T[l]=a})}function fa(){y.forEach(e=>{e.tracks=[]}),w.forEach(e=>{e.indexSequence.forEach(n=>{y[Math.abs(n)].tracks.push(e.id)})}),y.forEach(e=>{e.tracks!==void 0&&(e.degree=e.tracks.length)})}function pa(){let e=w[0],n=w.slice(1,w.length);if(v.showReads&&E.length>0&&(n=n.concat(E)),ha(n,e),v.showReads&&E.length>0){for(let t=0;t<E.length;t++)if(E[t]!==n[t+w.length-1])throw new Error("Read inequality")}}function ha(e,n){let t=new Map,r=n?new Set(n.sequence):null;for(let o of e){let s=o.sequence;for(let a=0;a<s.length;a+=1){let l=O(s[a]),i=oe(l);if(r&&!r.has(l)){let u=a>0?oe(s[a-1]).order:void 0,f=a<s.length-1?oe(s[a+1]).order:void 0;(u===void 0||u<i.order)&&(f===void 0||i.order<f)&&go(t,l,k(s[a])?1:-1),(u===void 0||u>i.order)&&(f===void 0||i.order>f)&&go(t,l,k(s[a])?-1:1)}}}for(let o of e)o.sequence.forEach((s,a)=>{let l=t.get(O(s));l!==void 0&&l>0&&(o.sequence[a]=Ao(s),o.indexSequence[a]=-o.indexSequence[a])});t.forEach((o,s)=>{if(o>0){let a=oe(s);a.seq=To(a.seq),a.switched=!0}})}function go(e,n,t){e.set(n,(e.get(n)??0)+t)}function ma(){let e=0,n=20,t=-1,r=y.slice();r.sort(fn);let o=ba();r.forEach(s=>{s.order>=0&&(s.order>t&&(t=s.order,e=n+10*o[s.order]),s.x=e,n=Math.max(n,e+40+s.pixelWidth))})}function ba(){let e=[],n=[],t=[],r=[];for(let o=0;o<=Q;o+=1)e.push(0),n.push(0),t.push(0);w.forEach(o=>{for(let s=1;s<o.path.length;s+=1){let a=o.path[s],l=o.path[s-1];if(a.order===l.order)a.isForward?e[a.order]=e[a.order]+1:n[a.order]=n[a.order]+1;else{let i=Math.abs((a.y??0)-(l.y??0));t[a.order]=Math.max(i/17.5,t[a.order])}}}),r.push(Math.max(0,e[0]-1));for(let o=1;o<=Q;o+=1)r.push(Math.max(Math.max(0,e[o]-1)+Math.max(0,n[o-1]-1),t[o]));return r}function ga(){let e,n,t,r,o,s,a=[],l=[];for(let i=0;i<=Q;i+=1){me[i]=[],l[i]=new Map,a[i]=[];for(let u=0;u<w.length;u+=1)a[i][u]=null}w.forEach((i,u)=>{n=Math.abs(i.indexSequence[0]),t=un(i.indexSequence[0]),r=y[n],i.path=[],i.path.push({order:r.order,lane:null,isForward:t,node:n}),X(r.order,n,u,0,a,l),e=1;for(let f=1;f<i.sequence.length;f+=1)if(o=r,s=t,n=Math.abs(i.indexSequence[f]),t=un(i.indexSequence[f]),r=y[n],r.order>o.order){s||(i.path.push({order:o.order,lane:null,isForward:!0,node:null}),X(o.order,null,u,e,a,l),e+=1);for(let c=o.order+1;c<r.order;c+=1)i.path.push({order:c,lane:null,isForward:!0,node:null}),X(c,null,u,e,a,l),e+=1;t?(i.path.push({order:r.order,lane:null,isForward:!0,node:n}),X(r.order,n,u,e,a,l),e+=1):(i.path.push({order:r.order,lane:null,isForward:!0,node:null}),X(r.order,null,u,e,a,l),e+=1,i.path.push({order:r.order,lane:null,isForward:!1,node:n}),X(r.order,n,u,e,a,l),e+=1)}else if(r.order<o.order){s&&(i.path.push({order:o.order,lane:null,isForward:!1,node:null}),X(o.order,null,u,e,a,l),e+=1);for(let c=o.order-1;c>r.order;c-=1)i.path.push({order:c,lane:null,isForward:!1,node:null}),X(c,null,u,e,a,l),e+=1;t?(i.path.push({order:r.order,lane:null,isForward:!1,node:null}),X(r.order,null,u,e,a,l),e+=1,i.path.push({order:r.order,lane:null,isForward:!0,node:n}),X(r.order,n,u,e,a,l),e+=1):(i.path.push({order:r.order,lane:null,isForward:!1,node:n}),X(r.order,n,u,e,a,l),e+=1)}else t!==s?(i.path.push({order:r.order,lane:null,isForward:t,node:n}),X(r.order,n,u,e,a,l),e+=1):(i.path.push({order:r.order,lane:null,isForward:!t,node:null}),X(r.order,null,u,e,a,l),e+=1,i.path.push({order:r.order,lane:null,isForward:t,node:n}),X(r.order,n,u,e,a,l),e+=1)});for(let i=0;i<=Q;i+=1)xa(me[i],i)}function X(e,n,t,r,o,s){let a={trackID:t,segmentID:r,compareToFromSame:o[e][t]??null},l=n===null?void 0:s[e].get(n);if(l===void 0){let i={type:"single",node:n,tracks:[a]};me[e].push(i),n!==null&&s[e].set(n,i)}else l.type="multiple",l.tracks.push(a);o[e][t]=a}function ya(e,n){let t;e.forEach(r=>{r.idealLane=0,r.tracks.forEach(o=>{if(o.segmentID===0)o.idealLane=o.trackID,o.idealY=null;else if(w[o.trackID].path[o.segmentID-1].order===n-1)o.idealLane=w[o.trackID].path[o.segmentID-1].lane??void 0,o.idealY=w[o.trackID].path[o.segmentID-1].y;else if(o.segmentID<w[o.trackID].path.length-1&&w[o.trackID].path[o.segmentID+1].order===n-1)o.idealLane=w[o.trackID].path[o.segmentID+1].lane??void 0,o.idealY=w[o.trackID].path[o.segmentID+1].y;else{for(t=o.segmentID-1;t>=0&&w[o.trackID].path[t].order!==n-1;)t-=1;t<0?(o.idealLane=o.trackID,o.idealY=null):(o.idealLane=w[o.trackID].path[t].lane??void 0,o.idealY=w[o.trackID].path[t].y)}r.idealLane+=o.idealLane}),r.idealLane/=r.tracks.length})}function xa(e,n){let t=0,r=new Set,o=20,s=!1,a=-1;ya(e,n),e.sort(xo),e.forEach(l=>{l.node!==null?(y[l.node].topLane=t,s&&(o-=10),y[l.node].y=o,y[l.node].contentHeight=0,s=!1):(s?o-=25:o>20&&(o-=10),s=!0),l.tracks.sort(xo),l.tracks.forEach(i=>{i.lane=t,i.trackID===a&&l.node===null&&s&&(o+=10),w[i.trackID].path[i.segmentID].lane=t,w[i.trackID].path[i.segmentID].y=o,i.idealY!=null&&r.add(i.idealY-o),t+=1,o+=w[i.trackID].width,l.node!==null&&(y[l.node].contentHeight+=w[i.trackID].width),a=i.trackID}),o+=25}),Ma(e,r)}function Ma(e,n){let t=0,r=Number.MAX_SAFE_INTEGER;n.forEach(o=>{yo(e,o)<r&&(r=yo(e,o),t=o)}),e.forEach(o=>{o.node!==null&&(y[o.node].y+=t),o.tracks.forEach(s=>{let a=w[s.trackID].path[s.segmentID];a.y=a.y+t})})}function wa(e,n){e.order<0||me[e.order]===void 0||(me[e.order].forEach(t=>{if(t.node!==null){let r=y[t.node];r!==e&&r.y>e.y&&(r.y+=n,t.tracks.forEach(o=>{let s=w[o.trackID].path[o.segmentID];s.y=s.y+n}))}else t.tracks.forEach(r=>{let o=w[r.trackID].path[r.segmentID];o.y>=e.y&&(o.y=o.y+n)})}),Ae[e.order].length>0&&Ae[e.order].forEach(t=>{y[t]!==e&&y[t].y>e.y&&(y[t].y+=n)}))}function yo(e,n){let t=0;return e.forEach(r=>{r.tracks.forEach(o=>{o.idealY!=null&&w[o.trackID].type!=="read"&&(t+=Math.abs(o.idealY-n-w[o.trackID].path[o.segmentID].y)*w[o.trackID].width)})}),t}function xo(e,n){return e.idealLane!==void 0?n.idealLane!==void 0?e.idealLane<n.idealLane?-1:e.idealLane>n.idealLane?1:0:-1:n.idealLane!==void 0?1:0}function fn(e,n){if(!e)return n?-1:0;if(!n)return 1;if(e.order!==void 0){if(n.order!==void 0){if(e.order<n.order)return-1;if(e.order>n.order)return 1;if(e.y!==void 0&&n.y!==void 0){if(e.y<n.y)return-1;if(e.y>n.y)return 1}return 0}return-1}return n.order!==void 0?1:0}function Sa(){let e,n,t={};Fn.forEach(r=>{let o=0;for(;o<w.length&&w[o].name!==r.track;)o+=1;o<w.length&&(e=0,w[o].path.forEach(s=>{s.node!==null&&(t={},n=e+y[s.node].sequenceLength-1,e>=r.start&&e<=r.end&&(t.start=0),e<r.start&&n>=r.start&&(t.start=r.start-e),n<=r.end&&n>=r.start&&(t.end=n-e,n<r.end&&(t.continue=!0)),n>r.end&&e<=r.end&&(t.end=r.end-e),t.start!==void 0&&(t.type=r.type,t.name=r.name,s.features===void 0&&(s.features=[]),s.features.push(t)),e=n+1)}))})}function Na(){let e=!0,n=4,t=v.trackWidth;w.forEach(r=>{r.freq!==void 0?r.width=Math.round((Math.log(r.freq)+1)*n):(r.width=t,r.type!==void 0&&r.type==="read"&&(r.width=n)),r.width!==n&&(e=!1)}),e&&w.forEach(r=>{r.freq!==void 0&&(r.width=t)})}function La(e){let n=e.path[0],t=y[n.node],r=e.firstNodeOffset??0;return Yn(t,n.isForward?r:t.sequenceLength-r)}function Ra(e){let n=e.path[e.path.length-1],t=y[n.node],r=e.finalNodeCoverLength??0;return Yn(t,n.isForward?r:t.sequenceLength-r)}function Co(e,n){return n>e.sequenceLength?null:Yn(e,n)}function Yn(e,n){let[t,r]=Ta(e);if(e.sequenceLength===0)return t;let o=Math.min(Math.max(n,0),e.sequenceLength);return t+o/e.sequenceLength*(r-t)}function Ea(){let e,n,t,r,o,s,a,l,i;for(let c=0;c<=Q;c+=1)Hn.push(0),Xn.push(0);let u=[],f=[];y.forEach(c=>{c.x!==void 0&&(u[c.order]=c.x,f[c.order]===void 0?f[c.order]=c.x+c.pixelWidth:f[c.order]=Math.max(f[c.order],c.x+c.pixelWidth))}),w.sort(Js),w.forEach(c=>{a="plain",o=Ee(c,a),s=ks(c);let d=a,h=()=>(d!==a&&(d=a,o=Ee(c,a)),o);t=c.path[0].y,c.type!=="read"?c.sequence[0].startsWith("-")?e=f[c.path[0].order]+20:e=u[c.path[0].order]-20:e=La(c);for(let p=0;p<c.path.length;p+=1)c.path[p].y===t?c.path[p].features!==void 0&&(i=p>0&&c.path[p-1].order===c.path[p].order,l=wo(c.path[p],u[c.path[p].order],f[c.path[p].order],a,c,e,t,o,i),a=l.highlight,e=l.xStart):(c.path[p-1].isForward?n=f[c.path[p-1].order]:n=u[c.path[p-1].order],n!==e&&(o=h(),G.rectangles.push({xStart:Math.min(e,n),yStart:t,xEnd:Math.max(e,n),yEnd:t+c.width-1,color:o,alpha:s,id:c.id,name:c.name,type:c.type})),c.path[p].order-1===c.path[p-1].order?(e=n,n=u[c.path[p].order],r=c.path[p].y,o=h(),G.curves.push({xStart:e,yStart:t,xEnd:n+1,yEnd:r,width:c.width,color:o,alpha:s,laneChange:Math.abs(c.path[p].lane-c.path[p-1].lane),id:c.id,name:c.name,type:c.type,nodeStart:c.path[p-1]?.node,nodeEnd:c.path[p]?.node}),e=n,t=r):c.path[p].order+1===c.path[p-1].order?(e=n,n=f[c.path[p].order],r=c.path[p].y,o=h(),G.curves.push({xStart:e+1,yStart:t,xEnd:n,yEnd:r,width:c.width,color:o,alpha:s,laneChange:Math.abs(c.path[p].lane-c.path[p-1].lane),id:c.id,name:c.name,type:c.type,nodeStart:c.path[p-1]?.node,nodeEnd:c.path[p]?.node}),e=n,t=r):c.path[p-1].isForward?(r=c.path[p].y,No(1,n,t,r,c.width,o,c.id,c.path[p].order,c.type,c.name),e=f[c.path[p].order],t=c.path[p].y):(r=c.path[p].y,No(-1,n,t,r,c.width,o,c.id,c.path[p].order,c.type,c.name),e=u[c.path[p].order],t=c.path[p].y),c.path[p].features!==void 0&&(i=c.path[p-1].order===c.path[p].order,l=wo(c.path[p],u[c.path[p].order],f[c.path[p].order],a,c,e,t,o,i),a=l.highlight,e=l.xStart));c.type!=="read"?c.path[c.path.length-1].isForward?n=f[c.path[c.path.length-1].order]+20:n=u[c.path[c.path.length-1].order]-20:n=Ra(c),G.rectangles.push({xStart:Math.min(e,n),yStart:t,xEnd:Math.max(e,n),yEnd:t+c.width-1,color:o,alpha:s,id:c.id,name:c.name,type:c.type})})}var Aa=1e9;function Mo(e,n){cn=new Map;let t=new Map,r=v.ignoreStrand;for(let d of e){let h=d.indexSequence;if(!(!h||h.length<2))for(let p=0;p<h.length-1;p+=1){let m=h[p],b=h[p+1],g=Math.abs(m),M=Math.abs(b),R=y[g],N=y[M];if(!R||!N)continue;let S;if(r){let L=`${m}>${b}`,A=`${-b}>${-m}`;S=L<A?L:A}else S=`${m}>${b}`;let x=t.get(S);x===void 0?t.set(S,{sSigned:m,dSigned:b,sName:m<0?`-${R.name}`:R.name,dName:b<0?`-${N.name}`:N.name,count:1,sourceTrackID:d.sourceTrackID}):x.count+=1}}let o=0;for(let d of t.values())d.count>o&&(o=d.count);let s=dn,a=60,l=d=>o<=1?s:s+(a-s)*Math.sqrt(d/o),i=[],u=0;for(let d of t.values()){let h=Aa+u,p=`${d.count.toLocaleString()} ${n}${d.count===1?"":"s"}: Node ${d.sName} \u2192 Node ${d.dName}`;cn.set(h,{count:d.count,label:p});let m=y[Math.abs(d.sSigned)]?.sequenceLength??0;i.push({id:h,sourceTrackID:d.sourceTrackID,type:"read",name:p,sequence:[d.sName,d.dName],indexSequence:[],path:[],width:l(d.count),firstNodeOffset:m,finalNodeCoverLength:0,sequenceNew:[{nodeName:d.sName,mismatches:[]},{nodeName:d.dName,mismatches:[]}],mapping_quality:0,is_secondary:!1,sample_name:null,read_group:null,score:d.count}),u+=1}let f="?",c=0;for(let d of y){if(!d)continue;let h=d.contentHeight;h>c&&(c=h,f=d.name)}return He&&console.log(`[coarsened] ${n}s_in=${e.length} edges_out=${i.length} tallestNodeBeforePlace=${f}(${c.toFixed(1)})`),i}function wo(e,n,t,r,o,s,a,l,i){let u=r,f,c,d,h;n-=8,t+=8;let p=y[e.node].sequenceLength;return e.features.sort((m,b)=>m.start-b.start),e.features.forEach(m=>{u!==m.type&&(f=Ee(o,u),e.isForward?(d=n+Math.round(m.start*(t-n+1)/p),i&&d<n+8&&(h=n+Math.round((m.end+1)*(t-n+1)/p)-1,c=Ee(o,m.type),G.featureRectangles.push({xStart:d,yStart:a,xEnd:h,yEnd:a+o.width-1,color:c,id:o.id,name:o.name,type:o.type})),d>s+1&&G.featureRectangles.push({xStart:s,yStart:a,xEnd:d-1,yEnd:a+o.width-1,color:f,id:o.id,name:o.name,type:o.type})):(d=t-Math.round(m.start*(t-n+1)/p),i&&d>t-8&&(h=t-Math.round((m.end+1)*(t-n+1)/p)-1,c=Ee(o,m.type),G.featureRectangles.push({xStart:h,yStart:a,xEnd:d,yEnd:a+o.width-1,color:c,id:o.id,name:o.name,type:o.type})),s>d+1&&G.featureRectangles.push({xStart:d+1,yStart:a,xEnd:s,yEnd:a+o.width-1,color:f,id:o.id,name:o.name,type:o.type})),s=d,u=m.type),(m.end<p-1||m.continue===void 0)&&(f=Ee(o,u),e.isForward?(h=n+Math.round((m.end+1)*(t-n+1)/p)-1,G.featureRectangles.push({xStart:s,yStart:a,xEnd:h,yEnd:a+o.width-1,color:f,id:o.id,name:o.name,type:o.type})):(h=t-Math.round((m.end+1)*(t-n+1)/p)-1,G.featureRectangles.push({xStart:h,yStart:a,xEnd:s,yEnd:a+o.width-1,color:f,id:o.id,name:o.name,type:o.type})),s=h+1,u="plain")}),{xStart:s,highlight:u}}var So=7;function No(e,n,t,r,o,s,a,l,i,u){let f=e===1?Xn:Hn,c=10*f[l],d=n+e*(c+5),h=So,p=Math.min(So,o),m=Math.min(t,r),b=Math.max(t,r),g=e===1?n:d-1,M=e===1?d:n,R=L=>{G.verticalRectangles.push({xStart:g,yStart:L,xEnd:M,yEnd:L+o-1,color:s,id:a,name:u,type:i})};R(t),G.verticalRectangles.push({xStart:e===1?d+h:d-h-p,yStart:m+o+h-1,xEnd:e===1?d+h+p-1:d-h-1,yEnd:b-h+1,color:s,id:a,name:u,type:i}),R(r);let N=d+e*h,S=d+e*(h+p),x=`M ${d} ${b}`;x+=` Q ${N} ${b} ${N} ${b-h}`,x+=` H ${S}`,x+=` Q ${S} ${b+o} ${d} ${b+o}`,x+=" Z ",G.corners.push({path:x,color:s,id:a,type:i}),x=`M ${d} ${m}`,x+=` Q ${S} ${m} ${S} ${m+o+h}`,x+=` H ${N}`,x+=` Q ${N} ${m+o} ${d} ${m+o}`,x+=" Z ",G.corners.push({path:x,color:s,id:a,type:i}),f[l]+=1}function Ta(e){let n=e.x-4,t=e.x+e.pixelWidth+4;return[n,t]}function Lo(){switch(y.forEach(e=>{e.sequenceLength===void 0&&(e.sequenceLength=e.seq.length)}),v.nodeWidthOption){case"compressed":y.forEach(e=>{e.width=1+Math.log(e.sequenceLength)/Math.log(2),e.pixelWidth=Math.round((e.width-1)*8.401)});break;case"small":y.forEach(e=>{e.width=e.sequenceLength/100,e.pixelWidth=Math.round((e.width-1)*8.401)});break;case"fixed":y.forEach(e=>{e.width=10,e.pixelWidth=Math.round(e.width*8.401)});break;case"normal":{let e=v.charWidth;y.forEach(n=>{n.width=n.sequenceLength;let t=n.seq?.length??1;n.pixelWidth=Math.round(e*t)});break}default:throw new Error(`${v.nodeWidthOption} not implemented`)}}function Ca(){let e,n,t=Array.from(y,()=>new Set),r=Array.from(y,()=>new Set),o;v.showReads&&E.length>0?o=w.concat(E):o=w,o.forEach(l=>{for(let i=0;i<l.sequence.length;i+=1)k(l.sequence[i])?(e=O(l.sequence[i]),i>0?(n=l.sequence[i-1],k(n)?r[I.get(e)].add(O(n)):(r[I.get(e)].add(n),r[I.get(e)].add(ue(n)))):l.type==="haplotype"&&r[I.get(e)].add("None"),i<l.sequence.length-1?(n=l.sequence[i+1],k(n)?t[I.get(e)].add(O(n)):(t[I.get(e)].add(n),t[I.get(e)].add(ue(n)))):l.type==="haplotype"&&t[I.get(e)].add("None")):(i>0?(e=l.sequence[i-1],t[I.get(l.sequence[i])].add(e),k(e)&&t[I.get(l.sequence[i])].add(O(e))):l.type==="haplotype"&&t[I.get(l.sequence[i])].add("None"),i<l.sequence.length-1?(e=l.sequence[i+1],r[I.get(l.sequence[i])].add(e),k(e)&&r[I.get(l.sequence[i])].add(O(e))):l.type==="haplotype"&&r[I.get(l.sequence[i])].add("None"))});let s=t.map(l=>Array.from(l)),a=r.map(l=>Array.from(l));if(v.showReads&&E.length>0){let l=y.slice();l.sort(fn);let i=new Map,u=new Map,f=d=>{let h=i.get(d);if(h===void 0)throw new Error(`No merge offset recorded for node ${d}`);return h},c=d=>{let h=u.get(d);if(h===void 0)throw new Error(`No merge origin recorded for node ${d}`);return h};l.forEach(d=>{let h=We(I.get(d.name),s,a);if(h){let p=f(h)+oe(h).sequenceLength;i.set(d.name,p),i.set(ue(d.name),p),u.set(d.name,c(h)),u.set(ue(d.name),c(h))}else i.set(d.name,0),i.set(ue(d.name),0),u.set(d.name,d.name),u.set(ue(d.name),d.name)}),E.forEach(d=>{let h=d.sequenceNew??[];d.firstNodeOffset=(d.firstNodeOffset??0)+f(d.sequence[0]),d.finalNodeCoverLength=(d.finalNodeCoverLength??0)+f(d.sequence[d.sequence.length-1]);for(let p=d.sequence.length-1;p>=0;p-=1){let m=O(d.sequence[p]),b=We(I.get(m),s,a);if(b&&qn(I.get(b),s,a))if(p>0){d.sequence.splice(p,1);let g=oe(b).sequenceLength;h[p].mismatches.forEach(M=>{M.pos+=g}),h[p-1].mismatches=h[p-1].mismatches.concat(h[p].mismatches),h.splice(p,1)}else{d.sequence[0]=c(d.sequence[0]);let g=h[0],M=f(g.nodeName);g.mismatches.forEach(R=>{R.pos+=M}),g.nodeName=c(g.nodeName)}}})}for(let l=0;l<y.length;l+=1)if(qn(l,s,a)&&!We(l,s,a)){let i=l;for(;qn(i,s,a);)i=I.get(O(a[i][0])),y[l].sequenceLength+=y[i].sequenceLength,y[l].seq+=y[i].seq}w.forEach(l=>{for(let i=l.sequence.length-1;i>=0;i-=1){e=O(l.sequence[i]);let u=I.get(e);We(u,s,a)&&l.sequence.splice(i,1)}});for(let l=y.length-1;l>=0;l-=1)We(l,s,a)&&y.splice(l,1)}function We(e,n,t){if(n[e].length!==1||n[e][0]==="None")return!1;let r=O(n[e][0]),o=I.get(r);return t[o].length!==1||t[o][0]==="None"?!1:r}function qn(e,n,t){if(t[e].length!==1||t[e][0]==="None")return!1;let r=O(t[e][0]),o=I.get(r);return!(n[o].length!==1||n[o][0]==="None")}function Ia(e){let n=v.focusReadNames?new Set(v.focusReadNames):null;return e.filter(t=>t.is_secondary!==!0&&(t.mapping_quality??0)>=v.mappingQualityCutoff&&(n===null||t.name!==void 0&&n.has(t.name)))}function _a(e,n){let t=new Map;for(let r of e){let o=n(r),s=t.get(o);s===void 0?t.set(o,[r]):s.push(r)}return t}function va(e,n){return e.xStart<n.xStart?1:e.xStart>n.xStart?-1:e.yStart>n.yStart?1:e.yStart<n.yStart?-1:0}function Io(e,n){let t=_a(e.filter(s=>s.type===n),s=>`${s.nodeStart},${s.nodeEnd}`);t.forEach(s=>{let a=.3,l=.4/s.length;s.length<=5&&(a=.5,l=0),s.sort(va),s.forEach(i=>{let u,f;i.yStart<i.yEnd?(u=i.xStart+(i.xEnd-i.xStart)*(1-a),a+=l,f=i.xStart+(i.xEnd-i.xStart)*(1-a)):(u=i.xStart+(i.xEnd-i.xStart)*a,a+=l,f=i.xStart+(i.xEnd-i.xStart)*a);let c=`M ${i.xStart} ${i.yStart}`;c+=` C ${u} ${i.yStart} ${u} ${i.yEnd} ${i.xEnd} ${i.yEnd}`,c+=` V ${i.yEnd+i.width}`,c+=` C ${f} ${i.yEnd+i.width} ${f} ${i.yStart+i.width} ${i.xStart} ${i.yStart+i.width}`,c+=" Z",i.path=c})});let r=Array.from(t.keys());r.sort(function(s,a){let l=s.split(",").map(Number),i=a.split(",").map(Number),u=l[0]??0,f=l[1]??0,c=i[0]??0,d=i[1]??0;return u<c?-1:u>c?1:f<d?-1:f>d?1:0});let o=[];for(let s of r){let a=t.get(s);a!==void 0&&o.push(...a)}return o}function _o(e){let n=`M ${e.x-9} ${e.y} Q ${e.x-9} ${e.y-9} ${e.x} ${e.y-9}`,t=e.x,r=e.y-9;return e.width>1&&(t+=e.pixelWidth,n+=` L ${t} ${r}`),n+=` Q ${t+9} ${r} ${t+9} ${r+9}`,t+=9,r+=9,e.contentHeight>0&&(r+=e.contentHeight-0,n+=` L ${t} ${r}`),n+=` Q ${t} ${r+9} ${t-9} ${r+9}`,t-=9,r+=9,e.width>1&&(t-=e.pixelWidth,n+=` L ${t} ${r}`),n+=` Q ${t-9} ${r} ${t-9} ${r-9}`,t-=9,r-=9,e.contentHeight>0&&(r-=e.contentHeight-0,n+=` L ${t} ${r}`),n}function Ba(e){let n=new Map;for(let[t,r]of e.pathVisits??[])for(let o of r){let s=`${o.path} ${t}`,a=n.get(s);a?a.push(o.strand):n.set(s,[o.strand])}return n}function Oa(e){let n=new Map(e.nodes.map(s=>[s.id,s])),t=Ba(e),r=(e.paths??[]).map((s,a)=>({path:s,index:a})),o=r.findIndex(s=>Ne(s.path.name).name===e.referencePath);return o>0&&r.unshift(...r.splice(o,1)),r.flatMap(({path:s,index:a})=>{let l=Ne(s.name).name,i=new Map,u=s.nodeIds.flatMap(f=>{let c=n.get(f);if(!c)return[];let d=i.get(c.name)??0;i.set(c.name,d+1);let h=t.get(`${l} ${c.name}`)?.[d]??vn(c);return[Bn(c,h)]});return u.length>0?[{id:a,name:s.name,sequence:u,sourceTrackID:0,type:"haplotype"}]:[]})}function Pa(e){return e.nodes.map(n=>({name:n.id,seq:"",sequenceLength:Math.max(1,n.length)}))}function Te(e){return(e.paths?.length??0)>0}var Ga=600,Da=3,$a=15;function qa(e){return Math.max(Da,Math.min($a,Math.floor(Ga/Math.max(1,e))))}function vo(e){let n=e.paths??[],t=Oa(e),r=e.reads?Qt(e,e.reads,n.length):[],o=r.length>0?pt:ft,s=n.map((l,i)=>o(i,n.length)),a=Eo(Pa(e),t,r,{nodeWidthOption:"compressed",trackWidth:qa(t.length),trackColor:l=>l.type==="read"?Zt(l):s[l.id]});return a?{layout:a,pathColors:s}:void 0}function zn(e){let n=[];return e.nodes.forEach(t=>{t.order>=0&&n.push(t)}),n}function Bo(e,n,t,r){let o=e.y+e.contentHeight/2-n;return[{x:t,y:o},{x:r,y:o}]}function Oo(e,n){return{minY:e.bounds.minY-n,maxY:e.bounds.maxY-n}}function Po(e){let n=Te(e)?vo(e):void 0;if(!n)return;let{layout:t,pathColors:r}=n,o=t.bounds.minY,s={};for(let i of zn(t))s[i.name]=Bo(i,o,i.x,i.x+i.pixelWidth);let{minY:a,maxY:l}=Oo(t,o);return{nodePositions:s,tubeMap:{layout:t,yOffset:o,pathColors:r},extent:{minX:t.bounds.minX,maxX:t.bounds.maxX,minY:a,maxY:l}}}function Fa(e,n){let t=new Map(e.nodes.map(l=>[l.id,l])),r=new Map;for(let l of zn(n)){let i=r.get(l.order);i?i.push(l):r.set(l.order,[l])}let o=[],s=-1/0;for(let l of[...r.keys()].sort((i,u)=>i-u)){let i=r.get(l),u=1/0,f=-1/0;for(let c of i){let d=t.get(c.name);d&&_(d)&&(u=Math.min(u,d.stable.start),f=Math.max(f,d.stable.start+c.sequenceLength))}u===1/0&&(u=f=s),u<s&&(u=s,f=Math.max(s,f)),s=Math.max(s,f),o.push({order:l,x0:Math.min(...i.map(c=>c.x)),x1:Math.max(...i.map(c=>c.x+c.pixelWidth)),bp0:u,bp1:f})}let a=o.find(l=>l.bp0!==-1/0);for(let l of o)l.bp0===-1/0&&(l.bp0=l.bp1=a?.bp0??0);return o}function Wa(e){return e.nodes.some(_)&&Te(e)}function Go(e){let n=Wa(e)?vo(e):void 0;if(!n)return;let{layout:t,pathColors:r}=n,o=Fa(e,t),s=new Map(o.map(i=>[i.order,i])),a=t.bounds.minY,l={};for(let i of zn(t)){let u=s.get(i.order);l[i.name]=Bo(i,a,u.bp0,u.bp1)}return{nodePositions:l,tubeMap:{layout:t,yOffset:a,columns:o,pathColors:r},referenceAxis:!0,pixelRows:!0,extent:Oo(t,a)}}function Do(e){return e.sample??Ke(e.name)}function Ha(e){return e.haplotype!==void 0&&e.sample!==void 0?`${e.sample}#${e.haplotype}`:Do(e)}function Xa(e,n,t,r=!0){if(!t||!r)return{ids:e,complete:!0};let o=-1,s=-1,a=-1/0,l=1/0;if(e.forEach((u,f)=>{let c=n.get(u);c&&(c.end<=t.start&&c.end>a&&(a=c.end,o=f),c.start>=t.end&&c.start<l&&(l=c.start,s=f))}),o<0||s<0)return{ids:e,complete:!1};let i=e.slice(Math.min(o,s)+1,Math.max(o,s));return{ids:o<s?i:i.reverse(),complete:!0}}function $o(e,n,t){let r=e.paths??[],o=r.find(m=>Ne(m.name).name===e.referencePath)??r[0];if(!o||r.length<2)return;let s=new Map(e.nodes.map(m=>[m.id,m.length])),a=e.anchorPaths?.find(m=>m.name===Ne(o.name).name)?.start??0,l=n&&n.end>n.start?n:void 0,i=new Map,u=a;for(let m of o.nodeIds){let b=s.get(m)??0;i.has(m)||i.set(m,{start:u,end:u+b}),u+=b}let f=l===void 0||[...i.values()].some(m=>m.end<=l.start)&&[...i.values()].some(m=>m.start>=l.end),c=m=>{let{ids:b,complete:g}=Xa(m.nodeIds,i,l,f),M=[],R=0,N=0,S=0;for(let x of b){let L=s.get(x)??0,A=i.get(x),C=M.at(-1),P=C?.referenceStart,D=A&&P!==void 0&&S>=0&&P+C.bp===A.start,$=A&&P!==void 0&&S<=0&&A.end===P;!A&&C&&!C.onReference?C.bp+=L:D||$?(C.bp+=L,D||(C.referenceStart=A.start,C.reversed=!0),S=D?1:-1):(M.push({start:R,bp:L,onReference:A!==void 0,referenceStart:A?.start}),S=0),R+=L,A||(N+=L)}return{name:m.name,label:Ha(m),sample:Do(m),haplotype:m.haplotype,bp:R,offReferenceBp:N,complete:g,runs:M}},d=m=>m.nodeIds.some(b=>i.has(b)),h=new Set(r.filter(m=>m!==o&&d(m)).map(m=>m.name));return{origin:l?l.start:a,unit:t,reference:c(o),rows:r.filter(m=>m!==o&&(d(m)||!h.has(m.name))).map(c).sort((m,b)=>Number(b.complete)-Number(m.complete)||b.bp-m.bp||m.label.localeCompare(b.label))}}var ka=1.3;function Ya(e){let n=e.reference.bp;for(let t of e.rows)n=Math.max(n,t.bp);return{maxX:e.origin+n*ka,maxY:e.rows.length*20}}function qo(e,n){let t=we(e),r=$o(e,n);if(t.length===0||!r)return;let o=Se(t),s=[{label:r.reference.label,y:0},...r.rows.map((a,l)=>({label:a.label,y:(l+1)*20}))];return{nodePositions:o,rowLabels:s,referenceAxis:!0,pixelRows:!0,extent:Ya(r)}}var za="Force-directed layout",be=e=>e.nodes.some(_),Fo=e=>be(e)&&(e.paths?.length??0)>1,Wo=e=>be(e)&&Te(e),Ho=e=>be(e)&&e.nodes.some(Qe),jn=[{value:"auto",label:"Anchored",description:"x is reference bp, one row per stable rank. Needs rGFA tags or a reference path.",run:to,cutMargins:!0,available:be,drawsLocally:be},{value:"samplerows",label:"Sample rows",description:"x is reference bp, one row per contributing assembly. Needs rGFA tags or a reference path.",run:io,cutMargins:!0,available:Ho,drawsLocally:Ho},{value:"walkrows",label:"Walk rows",description:"x is each walk's own bp: one bar per haplotype, sequence the reference also carries in blue and sequence it does not in purple, so a repeat expansion reads as bar length. Needs W or P lines.",run:qo,cutMargins:!1,available:Fo,drawsLocally:Fo},{value:"ordered",label:"Ordered",description:"x is reference order, not bp: every node gets room, bubbles read as lenses. Needs rGFA tags or a reference path.",run:ao,cutMargins:!1,available:be,drawsLocally:be},{value:"tubemap",label:"Tube map",description:"sequenceTubeMap: every path a coloured tube through boxed nodes, columns in node order and node widths log-scaled. Needs P or W lines.",run:Po,cutMargins:!1,available:Te,drawsLocally:Te},{value:"tubemapref",label:"Tube map on reference",description:"The tube map with each column at the reference bp it covers, so the tubes line up with the tracks around them. Needs P or W lines and a reference path.",run:Go,cutMargins:!0,available:Wo,drawsLocally:Wo},{value:"force",label:za,description:"OGDF FMMM, via the external Bandage engine.",run:()=>{},cutMargins:!1,available:()=>!0,drawsLocally:()=>!1}],Kc=jn.map(e=>e.value);function Qc(e,n){return!ja(e).drawsLocally(n)}function ja(e){return jn.find(n=>n.value===e)??jn[0]}var Y=ie(pn(),1);var F=ie(pn(),1);function ge(e,n){return(0,F.packAbgr)((0,F.abgrRed)(e),(0,F.abgrGreen)(e),(0,F.abgrBlue)(e),Math.round((0,F.abgrAlpha)(e)*n))}function eu(e,n){return(0,F.packAbgr)(Math.min(255,Math.round((0,F.abgrRed)(e)*n)),Math.min(255,Math.round((0,F.abgrGreen)(e)*n)),Math.min(255,Math.round((0,F.abgrBlue)(e)*n)),(0,F.abgrAlpha)(e))}var ko=new WeakMap;function Yo(e,n,t,r,o){let{scaleX:s,scaleY:a}=t,l=ko.get(e);if(l?.graph===n&&l.scaleX===s&&l.scaleY===a&&l.deletions===r&&l.version===o)return l.curves;let i=new Map;for(let u=0;u<n.edges.length;u++){let f=n.edges[u],c=e[f.from],d=e[f.to];if(c?.length&&d?.length){let h=r?.get(u);i.set(u,De(c,d,Tt(f),0,0,t,h?Rn(e,h):[]))}}return ko.set(e,{graph:n,scaleX:s,scaleY:a,deletions:r,version:o,curves:i}),i}var Va=(0,Y.packAbgr)(119,119,119,217),zo=(0,Y.packAbgr)(136,136,136,217),Ua=(0,Y.packAbgr)(24,24,28,240),Ka=2.2,Qa=(0,Y.packAbgr)(30,30,36,245),Za=1.8,Ce=.18,Ja=11,ei=1.2,ni=1e-6,ti=7,oi=3,ri=3,si=1.5,ai=.5,ii=.45,Vn=12,li=1e-6,de=new Float64Array(64),ye=new Float64Array(64),xe=new Float64Array(64);function ci(e,n,t,r){let o=e.length*Vn+1;de.length<o&&(de=new Float64Array(o),ye=new Float64Array(o),xe=new Float64Array(o));let s=e[e.length-1];de[0]=s.x1,ye[0]=s.y1,xe[0]=0;let a=1,l=0;for(let i=e.length-1;i>=0&&l<r;i--){let u=e[i];for(let f=1;f<=Vn&&l<r;f++){let c=1-(f/Vn)**2,d=1-c,h=d*d*d,p=3*d*d*c,m=3*d*c*c,b=c*c*c,g=h*u.x0+p*u.cx0+m*u.cx1+b*u.x1,M=h*u.y0+p*u.cy0+m*u.cy1+b*u.y1;l+=Math.hypot(g-de[a-1],(M-ye[a-1])*t)*n,de[a]=g,ye[a]=M,xe[a]=l,a++}}return a}function jo(e,n){let t=1;for(;t<e-1&&xe[t]<n;)t++;let r=xe[t-1],o=xe[t],s=o>r?Math.min(1,Math.max(0,(n-r)/(o-r))):0;return{x:de[t-1]+(de[t]-de[t-1])*s,y:ye[t-1]+(ye[t]-ye[t-1])*s}}function ui({curves:e,edgeHalfWidth:n,nodeHalfWidth:t,scale:r,yToX:o,color:s,fromSegments:a,toSegments:l}){let i=ti+oi*n,u=ci(e,r,o,t+i),f=xe[u-1],c=Math.min(t,f),d=f<li?di(a,r,o):f-c,h=Math.min(1,d/i);if(h<ai)return;let p=jo(u,c),m=jo(u,c+i*h),b=p.x-m.x,g=(p.y-m.y)*o,M=e[e.length-1];return{x:p.x,y:p.y,angle:Math.hypot(b,g)>0?Math.atan2(g,b):Ni(M,o,fi(M,l)),length:i*h,halfWidth:(ri+si*n)*h,color:s}}function di(e,n,t){let r=0;for(let o=1;o<e.length;o++)r+=Math.hypot(e[o].x-e[o-1].x,(e[o].y-e[o-1].y)*t);return r*n}function fi(e,n){if(n.length<2)return;let t=n[0],r=n[n.length-1];return Math.hypot(e.x1-t.x,e.y1-t.y)<=Math.hypot(e.x1-r.x,e.y1-r.y)?{from:t,to:n[1]}:{from:r,to:n[n.length-2]}}function pi(e,n=1){let t=n===1?e:e.map(o=>({x:o.x,y:o.y*n})),r=[];for(let o=0;o<t.length;o++){let s=0,a=0;if(o===0){let l=t[1].x-t[0].x,i=t[1].y-t[0].y,u=Math.hypot(l,i);u>0&&(s=-i/u,a=l/u)}else if(o===t.length-1){let l=t[o].x-t[o-1].x,i=t[o].y-t[o-1].y,u=Math.hypot(l,i);u>0&&(s=-i/u,a=l/u)}else{let l=t[o].x-t[o-1].x,i=t[o].y-t[o-1].y,u=Math.hypot(l,i),f=t[o+1].x-t[o].x,c=t[o+1].y-t[o].y,d=Math.hypot(f,c);if(u>0&&d>0){let h=-i/u,p=l/u,m=-c/d,b=f/d;s=(h+m)/2,a=(p+b)/2;let g=h*s+p*a;g>.1&&(s/=g,a/=g)}else u>0?(s=-i/u,a=l/u):d>0&&(s=-c/d,a=f/d)}r.push({nx:s,ny:a})}return r}function hi(e,n,t,r=1){return e.map((o,s)=>({x:o.x+n[s].nx*t,y:o.y+n[s].ny*t/r}))}function hn(e,n,t,r){return(0,Y.packAbgr)(Math.round(e*255),Math.round(n*255),Math.round(t*255),Math.round(r*255))}function mn(e,n,t){let r=(1-Math.abs(2*t-1))*n,o=r*(1-Math.abs(e/60%2-1)),s=t-r/2,a=0,l=0,i=0;return e<60?(a=r,l=o):e<120?(a=o,l=r):e<180?(l=r,i=o):e<240?(l=o,i=r):e<300?(a=o,i=r):(a=r,i=o),[a+s,l+s,i+s]}var mi=[[68,1,84],[59,82,139],[33,145,140],[94,201,98],[253,231,37]],bi=[[237,137,44],[158,42,122]],gi=[[220,50,50],[50,120,220]];function Un(e,n){let t=Math.max(0,Math.min(1,n))*(e.length-1),r=Math.min(e.length-2,Math.floor(t)),o=t-r,s=e[r],a=e[r+1];return(0,Y.packAbgr)(Math.round(s[0]+(a[0]-s[0])*o),Math.round(s[1]+(a[1]-s[1])*o),Math.round(s[2]+(a[2]-s[2])*o),255)}function yi(e,n){let[t,r,o]=mn(ut(e),yn,xn);return hn(t,r,o,n)}var xi=300,Mi=.7,wi=.5,Ko=[60,65,72],Si=(0,Y.packAbgr)(...Ko,255),du=`rgb(${Ko.join(", ")})`;function fu(e,n){let t=Jt(e);if(n&&n.end>n.start)return{start:n.start,span:n.end-n.start,midpoints:t};let r=1/0,o=-1/0;for(let s of t.values())r=Math.min(r,s),o=Math.max(o,s);return{start:r===1/0?0:r,span:o>r?o-r:1,midpoints:t}}function Kn(e){let n=1/0,t=-1/0,r=1/0,o=-1/0,s=0;for(let a of e.nodes)a.stable!==void 0&&a.stable.rank>s&&(s=a.stable.rank),a.depth<n&&(n=a.depth),a.depth>t&&(t=a.depth),a.length<r&&(r=a.length),a.length>o&&(o=a.length);return{minDepth:n,maxDepth:t,minLength:r,maxLength:o,maxRank:s,nodeCount:e.nodes.length}}function Qn(e,n,t,r){switch(t){case"random":return yi(e.id,1);case"depth":return Un(mi,r.maxDepth>r.minDepth?(e.depth-r.minDepth)/(r.maxDepth-r.minDepth):.5);case"node-length":return Un(gi,r.maxLength>r.minLength?(e.length-r.minLength)/(r.maxLength-r.minLength):.5);case"stable-rank":return e.stable===void 0?(0,Y.packAbgr)(160,160,160,255):e.stable.rank===0?(0,Y.packAbgr)(52,152,219,255):Un(bi,r.maxRank>1?(e.stable.rank-1)/(r.maxRank-1):0);case"reference-position":{let o=r.referenceRamp,s=o?.midpoints.get(e.id);if(o===void 0||s===void 0)return(0,Y.packAbgr)(160,160,160,255);if(e.stable!==void 0&&e.stable.rank>0)return Si;let a=Math.max(0,Math.min(1,(s-o.start)/o.span)),[l,i,u]=mn(a*xi,Mi,wi);return hn(l,i,u,1)}case"grey":return(0,Y.packAbgr)(160,160,160,255);case"rainbow":{let o=r.nodeCount>1?n/r.nodeCount*360:0,[s,a,l]=mn(o,.75,.5);return hn(s,a,l,1)}default:return(0,Y.packAbgr)(52,152,219,255)}}function Ni(e,n=1,t){let r=[[e.cx1,e.cy1],[e.cx0,e.cy0],[e.x0,e.y0]];for(let[o,s]of r){let a=e.x1-o,l=(e.y1-s)*n;if(Math.hypot(a,l)>0)return Math.atan2(l,a)}return t?Math.atan2((t.to.y-t.from.y)*n,t.to.x-t.from.x):0}function Li(e,n){let t=!1;for(let r=0;r<e.length-1&&!t;r++){let o=e[r],s=e[r+1];t=Math.max(o.x,s.x)>=n.minX&&Math.min(o.x,s.x)<=n.maxX&&Math.max(o.y,s.y)>=n.minY&&Math.min(o.y,s.y)<=n.maxY}return t}function Ri(e,n){for(let t of e){let r=Math.min(t.x0,t.cx0,t.cx1,t.x1),o=Math.max(t.x0,t.cx0,t.cx1,t.x1),s=Math.min(t.y0,t.cy0,t.cy1,t.y1),a=Math.max(t.y0,t.cy0,t.cy1,t.y1);if(o>=n.minX&&r<=n.maxX&&a>=n.minY&&s<=n.maxY)return!0}return!1}var Vo=new WeakMap;function Ei(e){let n=Vo.get(e);if(!n){let t=[],r=new Map,o=e.paths??[];for(let s=0;s<o.length;s++){let[a,l,i]=mn(dt(s,o.length),yn,xn);t.push(hn(a,l,i,.85)),r.has(o[s].name)||r.set(o[s].name,s)}n={pathColorByIndex:t,pathIndexByName:r},Vo.set(e,n)}return n}var Uo=new WeakMap,Ai=new Map;function Ti(e){let n=Uo.get(e);if(!n){n=new Map;for(let[t,r]of(e.paths??[]).entries())for(let o of new Set(r.nodeIds)){let s=n.get(o);s?s.push(t):n.set(o,[t])}Uo.set(e,n)}return n}function pu(e){let{nodePositions:n,graph:t,nodeById:r,colorScheme:o,contigThickness:s,connectorThickness:a,drawPaths:l,nodeWidth:i="uniform",highlight:u,axis:f,linearLayout:c,viewportBounds:d,referenceRamp:h,deletions:p,hiddenEdges:m,version:b=0}=e,g=i==="depth"?ht(t):0,M=f.scaleX,R=Je(f),N=Yo(n,t,f,p,b),S=[],x=new Map,L=[],A=new Map,C=[],P=new Map,D={...Kn(t),referenceRamp:h},$=new Map;if(o==="rainbow")for(let B=0;B<t.nodes.length;B++)$.set(t.nodes[B].id,B);let W=t.paths?.length??0,{pathColorByIndex:H,pathIndexByName:U}=Ei(t),j=B=>H[U.get(B)??-1]??zo,V=l&&W>1?Ti(t):Ai,Xe=M>(c?1:ii);for(let B=0;B<t.edges.length;B++){let q=t.edges[B],K=n[q.from],Z=n[q.to],J=N.get(B);if(!K?.length||!Z?.length||!J||m?.has(B))continue;let re=q.pathIds?.length??0,se=p?.get(B)!==void 0,ae=u?.edgeIndexes.has(B)??!1,ee=a/2*(se?Ka:1)*(ae?Za:1),fe=l&&re>0,ve=se?Ua:Va,Me=u?ae?Qa:ge(ve,Ce):ve,ke=r.get(q.to),gn=(ke?Mn(ke,s,i,g):s)/2,Be=(ze,je,ne,Pe)=>{let pe=ze===0&&je===0?J:Et(J,ze,je);if(d&&!Ri(pe,d))return;if(se)for(let xr of Rt(pe,Ja/M))C.push({curves:xr,thickness:ee,color:ne});else C.push({curves:pe,thickness:ee,color:ne});let rt=Xe&&Pe!==void 0?ui({curves:pe,edgeHalfWidth:ee,nodeHalfWidth:gn,scale:M,yToX:R,color:Pe,fromSegments:K,toSegments:Z}):void 0;rt&&L.push(rt)},Ye=C.length,Oe=L.length;if(!fe)Be(0,0,Me,Me);else{let ze=At(K,Z,re,f),je=Math.floor((re-1)/2);for(let ne=0;ne<re;ne++){let Pe=ze[ne],pe=j(q.pathIds[ne]);Be(Pe.x,Pe.y,u&&q.pathIds[ne]!==u.name?ge(pe,Ce):pe,ne===je?Me:void 0)}}let tt=C.length-Ye;tt>0&&P.set(B,{start:Ye,count:tt});let ot=L.length-Oe;ot>0&&A.set(B,{start:Oe,count:ot})}for(let[B,q]of Object.entries(n)){let K=r.get(B);if(!K||d&&!Li(q,d))continue;let Z=Qn(K,$.get(B)??0,o,D),J=u!==void 0&&!u.nodeIds.has(B),re=J?ge(Z,Ce):Z,_e=Mn(K,s,i,g),se=_e/2,ae=S.length,ee=V.get(B),fe=_e/W,ve=q.length>=2;if(ve&&ee?.length&&fe>=ei){let ke=pi(q,R),gn=1/Math.max(M,ni);for(let Be of ee){let Ye=(Be-(W-1)/2)*fe*gn,Oe=H[Be]??zo;S.push({points:hi(q,ke,Ye,R),thickness:fe/2,color:J?ge(Oe,Ce):Oe})}}else ve&&S.push({points:q,thickness:se,color:re});let Me=S.length-ae;Me>0&&x.set(B,{start:ae,count:Me})}return{nodeStrokes:S,nodeStrokeRuns:x,arrows:L,arrowRuns:A,edgeCurves:C,edgeCurveRuns:P}}var Jo=ie(pn(),1);var Ie=10,Ci=Ie+8,Zn=5;function yu(e,{colorScheme:n,referenceRamp:t,walkNodes:r}){let o={...Kn(e),referenceRamp:t},s=[];return e.nodes.forEach((a,l)=>{if(_(a)){let i=Qn(a,l,n,o);s.push({node:a.id,bp0:a.stable.start,bp1:a.stable.start+a.length,color:(0,Jo.abgrToCssRgba)(r&&!r.has(a.id)?ge(i,Ce):i)})}}),s.sort((a,l)=>a.bp0-l.bp0)}function Jn(e,n,t){let r=n*e.scale+e.translateX,o=t*e.scale+e.translateX;return r<o?[r,o]:[o,r]}var Qo=2;function Ii(e,n,t){let r=-n.translateX/n.scale,o=(t-n.translateX)/n.scale,s=Math.min(r,o),a=Math.max(r,o),l=0,i=0;for(let u of e)l+=Math.max(0,Math.min(u.bp1,s)-u.bp0),i+=Math.max(0,u.bp1-Math.max(u.bp0,a));return n.scale>=0?{left:Math.round(l),right:Math.round(i)}:{left:Math.round(i),right:Math.round(l)}}function xu(e,n,t,r){if(r<0||r>Ie+Qo)return;let o,s=Qo;for(let a of e){let[l,i]=Jn(n,a.bp0,a.bp1),u=Math.max(0,l-t,t-i);u<=s&&(s=u,o=a.node)}return o}function Mu(e,n){let t=e?.[Math.floor(e.length/2)];return t?n(t):void 0}function Zo(e,n,t){let r=Ie/2,o=n+t*(Zn+1);e.moveTo(n+t,r),e.lineTo(o,r-Zn/1.4),e.lineTo(o,r+Zn/1.4),e.closePath()}function wu(e,n,t,{width:r,lit:o,darkMode:s}){let a=s?"#ffffff":"#18181c";e.fillStyle=s?"#1f1f1f":"#ffffff",e.fillRect(0,0,r,Ci);for(let i of n){let[u,f]=Jn(t,i.bp0,i.bp1);f>=0&&u<=r&&(e.fillStyle=i.color,e.fillRect(u,0,Math.max(f-u,1),Ie))}let l=Ii(n,t,r);if((l.left>0||l.right>0)&&(e.beginPath(),l.left>0&&Zo(e,0,1),l.right>0&&Zo(e,r,-1),e.fillStyle=a,e.fill(),e.strokeStyle=s?"#1f1f1f":"#ffffff",e.lineWidth=1,e.stroke()),o){let[i,u]=Jn(t,o.start,o.end),f=Math.max(u-i,3),c=(i+u)/2-f/2;e.strokeStyle=a,e.lineWidth=1.5,e.strokeRect(c,.75,f,Ie-1.5),o.anchor&&(e.beginPath(),e.moveTo(c+f/2,Ie),e.lineTo(o.anchor.x,o.anchor.y),e.lineWidth=1,e.stroke())}}function Ru(e,n){let t=new Map(e.nodes.map(o=>[o.id,o])),r=new Map;n.nodes.forEach(o=>{let s=t.get(o.name);if(o.order>=0&&s&&_(s)){let{refName:a,start:l}=s.stable;(r.get(a)??r.set(a,[]).get(a)).push({bp0:l,bp1:l+o.sequenceLength,x0:o.x,x1:o.x+o.pixelWidth})}});for(let o of r.values())o.sort((s,a)=>s.bp0-a.bp0);return r}function et(e,n,t=!1){let r=0,o=e.length-1;for(;r<o;){let l=r+o>>1,i=e[l].bp1;(t?i<n:i<=n)?r=l+1:o=l}let s=e[r];if(n<=s.bp0)return s.x0;let a=s.bp1-s.bp0;return a>0?s.x0+(Math.min(n,s.bp1)-s.bp0)/a*(s.x1-s.x0):s.x0}function nt(e,n,t){let r=et(e,n),o=et(e,t,!0);return{x0:Math.min(r,o),x1:Math.max(r,o)}}function Eu(e){let n;for(let t of e.values())(!n||t.length>n.length)&&(n=t);return n}var er=5,_i=10,vi=110,Bi=3;function Oi(e){return e?"#b0b0b8":"#55555c"}function Pi(e){let n=10**Math.floor(Math.log10(e)),t=e/n;return(t<=1?1:t<=2?2:t<=5?5:10)*n}function nr(e){let n=[];for(let t of e)n.every(r=>Math.abs(t.at-r.at)>=(t.width+r.width)/2+_i)&&n.push(t);return n.map(({text:t,at:r})=>({text:t,at:r}))}function Gi(e,n,t,r){let o=e[0],s=e.at(-1);if(!o||!s)return;let a=n(o.x0),l=n(s.x1),i=s.bp1-o.bp0;if(l<=a||i<=0)return;let u=Pi(i/(l-a)*vi),f=[],c=[];for(let d=Math.ceil(o.bp0/u)*u;d<=s.bp1;d+=u){let h=Math.round(n(et(e,d)))+.5;if(h>=0&&h<=t){f.push(h);let p=d.toLocaleString("en-US");c.push({text:p,at:h,width:r(p)})}}return{spans:[{x0:a,x1:l}],ticks:f,labels:nr(c)}}function Di(e,n,t,r){let o=e.map(p=>({box:p,s0:n(p.x0),s1:n(p.x1)})).filter(({s0:p,s1:m})=>m>=0&&p<=t),s=o[0],a=o.at(-1);if(!s||!a)return;let l=(p,m)=>({text:p,at:m,width:r(p)}),i=(p,m)=>l(p.toLocaleString("en-US"),m),u=[],f=[],c=[],d=[];for(let{box:p,s0:m,s1:b}of o){u.push({x0:m,x1:b}),f.push(Math.round(m)+.5),b-m>=Bi&&f.push(Math.round(b)+.5);let g=l(sn(p.bp1-p.bp0),(m+b)/2);g.width<=b-m&&c.push(g),p!==s.box&&d.push(i(p.bp0,m))}let h=[i(s.box.bp0,s.s0),i(a.box.bp1,a.s1)];return{spans:u,ticks:f,labels:nr([...h,...c,...d])}}function $i(e,n,t,r,o=!1){return(o?Gi:Di)(e,n,t,r)}function Au(e,n,t,r,o=!1){let{x:s,width:a,darkMode:l}=t;e.font="10px sans-serif";let i=$i(n,s,a,c=>e.measureText(c).width,o);if(!i)return;let u=Oi(l),f=r+.5;e.strokeStyle=u,e.fillStyle=u,e.lineWidth=1,e.textAlign="center",e.textBaseline="top",e.beginPath();for(let{x0:c,x1:d}of i.spans)e.moveTo(Math.max(0,c),f),e.lineTo(Math.min(a,d),f);for(let c of i.ticks)e.moveTo(c,r),e.lineTo(c,r+er);e.stroke();for(let{text:c,at:d}of i.labels)e.fillText(c,d,r+er+2)}function Iu(e,n){let t=new Map(e.nodes.map(o=>[o.id,o])),r=[];return n.nodes.forEach(o=>{let s=t.get(o.name);o.order>=0&&s&&_(s)&&r.push({node:o.name,bp0:s.stable.start,bp1:s.stable.start+o.sequenceLength,x0:o.x,x1:o.x+o.pixelWidth})}),r}function _u(e,n,t,r){let o=[];for(let s of e){let a=t(s.x0),l=t(s.x1);l>=0&&a<=r&&o.push({node:s.node,top0:n(s.bp0),top1:n(s.bp1),bottom0:a,bottom1:l})}return o}var qi=3;function vu(e,n,t,r){if(n<=0||r<0||r>n)return;let o=r/n,s,a=qi;for(let l of e){let i=l.top0+(l.bottom0-l.top0)*o,u=l.top1+(l.bottom1-l.top1)*o,f=Math.max(0,Math.min(i,u)-t,t-Math.max(i,u));f<a&&(a=f,s=l.node)}return s}function tr(e,n,t){e.moveTo(n.top0,0),e.lineTo(n.top1,0),e.lineTo(n.bottom1,t),e.lineTo(n.bottom0,t),e.closePath()}function or(e,n,t){e.moveTo(n.top0,0),e.lineTo(n.bottom0,t),e.moveTo(n.top1,0),e.lineTo(n.bottom1,t)}function Bu(e,n,t,{highlightNode:r,darkMode:o}){if(t<=0||n.length===0)return;let s=o?"255,255,255":"0,0,0";e.beginPath();for(let l of n)tr(e,l,t);e.fillStyle=`rgba(${s},0.07)`,e.fill(),e.beginPath();for(let l of n)or(e,l,t);e.strokeStyle=`rgba(${s},0.35)`,e.lineWidth=.5,e.stroke();let a=n.find(l=>l.node===r);a&&(e.beginPath(),tr(e,a,t),e.fillStyle="rgba(255,192,203,0.5)",e.fill(),e.beginPath(),or(e,a,t),e.strokeStyle="#ff0000",e.lineWidth=1,e.stroke())}function Fi(e,n){let t=new Map(e.nodes.map(o=>[o.id,o])),r=[];return n.nodes.forEach((o,s)=>{let a=t.get(o.name);if(o.order>=0&&a&&_(a)){let l=a.stable.start;r.push({bp0:l,bp1:l+o.sequenceLength,index:s,node:o})}}),r.sort((o,s)=>o.bp0-s.bp0)}function Wi(e,n,t){let r=0,o=e.length-1;for(;r<o;){let s=r+o+1>>1;e[s].bp0<=n?r=s:o=s-1}return[e[r],e[r-1]].find(s=>s&&t.has(s.index)&&s.bp0<=n&&n<=s.bp1)}function Gu(e,n,t){let r=Fi(e,n),o=[];for(let s of n.tracks){let a=s.name?t.get(s.name):void 0;if(a?.length){let l=new Map;for(let i of s.path)i.node!==null&&i.y!==void 0&&!l.has(i.node)&&l.set(i.node,i.y);for(let i of a){let u=Wi(r,i.start,l);if(u){let{node:f,bp0:c,bp1:d}=u,h=p=>f.x+(Math.min(p,d)-c)/(d-c)*f.pixelWidth;o.push({x0:h(i.start),x1:h(i.end),y:l.get(u.index),height:s.width})}}}}return o}function Hi(e){return e?"rgba(255,255,255,0.9)":"rgba(0,0,0,0.8)"}function Du(e,n,{x:t,y:r,width:o,darkMode:s}){if(n.length!==0){e.fillStyle=Hi(s),e.beginPath();for(let a of n){let l=t(a.x0),i=Math.max(t(a.x1),l+1);if(i>=0&&l<=o){let u=r(a.y);e.rect(l,u,i-l,r(a.y+a.height)-u)}}e.fill()}}function Xi(e,n){let t=[],r=e.sequenceNew??[];return r.forEach((o,s)=>{let a=n.nodeMap.get(O(o.nodeName)),l=a===void 0?void 0:n.nodes[a],i=s;for(;i<e.path.length&&e.path[i].node!==a;)i+=1;let u=e.path[i]?.y;if(l&&u!==void 0){let f={readId:e.id,y:u,height:e.width,nodeY:l.y},c=d=>Co(l,d);for(let d of o.mismatches){let h=c(d.pos);if(h!==null){if(d.type==="insertion"){let p=s===0&&d.pos===e.firstNodeOffset||s===r.length-1&&d.pos===e.finalNodeCoverLength;t.push({...f,kind:"insertion",x:h,seq:d.seq,softClip:p})}else if(d.type==="deletion"&&d.length!==void 0){let p=c(d.pos+d.length);p!==null&&t.push({...f,kind:"deletion",x0:h,x1:p})}else if(d.type==="substitution"&&d.seq!==void 0){let p=c(d.pos+d.seq.length);p!==null&&t.push({...f,kind:"substitution",x0:h,x1:p,seq:d.seq})}}}}}),t}function rr(e){return e.reads.flatMap(n=>Xi(n,e))}function ur(e){let n=e.trim().split(/[\s,]+/),t=[],r=0,o=()=>Number(n[r++]);for(;r<n.length;){let s=n[r++];s==="M"||s==="L"?t.push({op:s,x:o(),y:o()}):s==="H"?t.push({op:s,x:o()}):s==="V"?t.push({op:s,y:o()}):s==="Q"?t.push({op:s,x1:o(),y1:o(),x:o(),y:o()}):s==="C"?t.push({op:s,x1:o(),y1:o(),x2:o(),y2:o(),x:o(),y:o()}):s==="Z"&&t.push({op:s})}return t}function dr(e){let n=1/0,t=-1/0;for(let r of e)"x"in r&&(n=Math.min(n,r.x),t=Math.max(t,r.x)),"x1"in r&&(n=Math.min(n,r.x1),t=Math.max(t,r.x1));return{x0:n,x1:t}}function sr(e,n,t){let r=ur(e);return{commands:r,color:n,alpha:t,...dr(r)}}function ar(e,n){let{shapes:t}=e,r=[...t.rectangles,...t.verticalRectangles].filter(a=>a.type===n).map(a=>({color:a.color,alpha:a.alpha??1,x0:a.xStart,x1:a.xEnd+1,y0:a.yStart,y1:a.yEnd+1})),o=Io(t.curves,n).map(a=>sr(a.path,a.color,a.alpha??1)),s=t.corners.filter(a=>a.type===n).map(a=>sr(a.path,a.color,1));return{rects:r,shapes:[...o,...s]}}function Xu(e){let n=[];return e.nodes.forEach(t=>{if(t.order>=0){let r=ur(_o(t));n.push({name:t.name,commands:r,...dr(r)})}}),{bounds:e.bounds,layers:[ar(e,"haplotype"),ar(e,"read")],nodes:n,mismatches:rr(e).filter(t=>t.kind!=="insertion"||!t.softClip)}}function ir(e,n,t,r){let o=0,s=0;for(let a of n)switch(a.op){case"M":o=a.x,s=a.y,e.moveTo(t(o),r(s));break;case"L":o=a.x,s=a.y,e.lineTo(t(o),r(s));break;case"H":o=a.x,e.lineTo(t(o),r(s));break;case"V":s=a.y,e.lineTo(t(o),r(s));break;case"Q":o=a.x,s=a.y,e.quadraticCurveTo(t(a.x1),r(a.y1),t(o),r(s));break;case"C":o=a.x,s=a.y,e.bezierCurveTo(t(a.x1),r(a.y1),t(a.x2),r(a.y2),t(o),r(s));break;case"Z":e.closePath();break}}function lr(e,n,t,r){let o=new Map;for(let s of n)if(t(s)){let a=`${s.color}|${s.alpha}`,l=o.get(a);l?l.push(s):o.set(a,[s])}for(let s of o.values()){e.fillStyle=s[0].color,e.globalAlpha=s[0].alpha,e.beginPath();for(let a of s)r(a);e.fill()}e.globalAlpha=1}var ki=12,Yi=.15;function zi(e){return Math.max(Yi,Math.min(1,e/ki))}function ku(e,n,t){let{x:r,y:o,width:s}=t,a=u=>r(u.x1)>=0&&r(u.x0)<=s;for(let u of n.layers)lr(e,u.rects,a,f=>{let c=r(f.x0),d=o(f.y0);e.rect(c,d,r(f.x1)-c,o(f.y1)-d)}),lr(e,u.shapes,a,f=>{ir(e,f.commands,r,o)});let l=t.darkMode?"#d0d0d0":"#000000",i=t.darkMode?"rgba(40,40,40,0.4)":"rgba(255,255,255,0.4)";e.lineWidth=2*Math.max(.25,Math.min(1,t.yScale));for(let u of n.nodes)if(a(u)){let f=u.name===t.highlightNode;e.beginPath(),ir(e,u.commands,r,o),e.fillStyle=f?"rgba(255,192,203,0.5)":i,e.strokeStyle=f?"#ff0000":l,e.fill(),e.globalAlpha=f?1:zi(r(u.x1)-r(u.x0)),e.stroke(),e.globalAlpha=1}Ui(e,n.mismatches,t)}var fr=12,ji=6;function Vi(e){return fr*e>=ji}function cr(e,{x:n,width:t}){let[r,o]=e.kind==="insertion"?[e.x,e.x]:[e.x0,e.x1];return n(o)>=0&&n(r)<=t}function Ui(e,n,t){let{x:r,y:o,yScale:s,darkMode:a}=t;if(!Vi(s)||n.length===0)return;let l=fr*s;e.fillStyle="grey",e.beginPath();for(let i of n)if(i.kind==="deletion"&&cr(i,t)){let u=r(i.x0),f=o(i.y);e.rect(u,f,r(i.x1)-u,o(i.y+i.height)-f)}e.fill(),e.font=`${l}px monospace`,e.textBaseline="alphabetic",e.fillStyle=a?"#ffffff":"#000000";for(let i of n)cr(i,t)&&(i.kind==="substitution"?e.fillText(i.seq,r(i.x0)+1,o(i.y+i.height)):i.kind==="insertion"&&e.fillText("*",r(i.x)-3,o(i.y+i.height)))}function ju(e,n){let t=[];for(let r of n){let o=e.get(r.refName),s=o?.[0],a=o?.at(-1);if(!o||!s||!a)continue;let l=Math.max(r.start,s.bp0),i=Math.min(r.end,a.bp1);i<=l||t.push({gene:r,...nt(o,l,i),exons:r.exons.filter(u=>u.end>l&&u.start<i).map(u=>nt(o,Math.max(u.start,l),Math.min(u.end,i)))})}return t}var pr=22,Ki=4,bn=8,Qi="11px sans-serif",Zi=8;function Vu(e){let n=e.flatMap(o=>[{x:o.x0,step:1},{x:o.x1,step:-1}]).sort((o,s)=>o.x-s.x||o.step-s.step),t=0,r=0;for(let{step:o}of n)t+=o,r=Math.max(r,t);return Math.min(Ki,r)}function Uu(e,n,t,r,o=0){let{x:s,width:a,darkMode:l}=t,i=Math.max(1,Math.floor((r-o)/pr)),u=[],f=l?"#e0e0e6":"#1c1c22";e.font=Qi,e.textAlign="center",e.textBaseline="alphabetic",e.fillStyle=f,e.strokeStyle=f,e.lineWidth=1;let c=n.map(d=>({g:d,s0:s(d.x0),s1:s(d.x1)})).filter(({s0:d,s1:h})=>h>=0&&d<=a).sort((d,h)=>d.s0-h.s0);for(let{g:d,s0:h,s1:p}of c){let m=d.gene.strand>0?" \u2192":d.gene.strand<0?"\u2190 ":"",b=d.gene.strand<0?`${m}${d.gene.name}`:`${d.gene.name}${m}`,g=(Math.max(h,0)+Math.min(p,a))/2,M=e.measureText(b).width/2,R=Math.min(h,g-M),N=Math.max(p,g+M),S=u.findIndex(L=>L+Zi<=R);S===-1&&(S=Math.min(u.length,i-1)),u[S]=Math.max(u[S]??-1/0,N);let x=r-S*pr-bn/2;e.beginPath(),e.moveTo(h,x+.5),e.lineTo(p,x+.5),e.stroke();for(let L of d.exons){let A=s(L.x0);e.fillRect(A,x-bn/2,Math.max(1,s(L.x1)-A),bn)}e.fillText(b,g,x-bn/2-2)}}function el(e,n){return Math.abs(n-e)<=.1*Math.max(e,n)}function nl(e,n){let t=e.length,r=n.length,o=Array.from({length:t+1},()=>new Array(r+1).fill(0)),s=Array.from({length:t+1},()=>new Array(r+1).fill(!1));for(let l=t-1;l>=0;l--)for(let i=r-1;i>=0;i--){let u=Math.abs(e[l]-n[i])+o[l+1][i+1],f=t-l>r-i?o[l+1][i]:1/0,c=r-i>t-l?o[l][i+1]:1/0;o[l][i]=Math.min(u,f,c),s[l][i]=o[l][i]===u}let a=new Array(t).fill(void 0);for(let l=0,i=0;l<t&&i<r;)s[l][i]?a[l++]=i++:t-l>r-i?l++:i++;return a}function mr(e,n){return e.map((t,r)=>({item:t,i:r})).sort((t,r)=>n(t.item)-n(r.item))}var gr=e=>e.spanningReads!==0;function br(e,n){return{...n,agrees:e.complete&&gr(n)?el(e.bp,n.bp):void 0}}function tl(e,n){let t=new Map,r=new Set;e.forEach((l,i)=>{let u=n.find(f=>f.haplotype!==void 0&&f.haplotype===l.haplotype);u&&(r.add(u),t.set(i,br(l,u)))});let o=n.filter(l=>l.haplotype===void 0),s=mr(e,l=>l.bp),a=[[s.filter(l=>l.item.complete),o.filter(gr)],[s,o]];for(let[l,i]of a){let u=l.filter(c=>!t.has(c.i)),f=mr(i.filter(c=>!r.has(c)),c=>c.bp);nl(u.map(c=>c.item.bp),f.map(c=>c.item.bp)).forEach((c,d)=>{if(c!==void 0){let h=u[d],p=f[c].item;r.add(p),t.set(h.i,br(h.item,p))}})}return t}function Zu(e,n){if(!n)return e;let t=new Map,r=new Map;for(let o of e)r.set(o.sample,[...r.get(o.sample)??[],o]);for(let[o,s]of r){let a=n[o];if(a)for(let[l,i]of tl(s,a))t.set(s[l],i)}return e.map(o=>({...o,call:t.get(o)}))}var yr=ie(Ue(),1),nd=(0,te.createSvgIcon)((0,yr.jsx)("path",{d:"M19.14 12.94c.04-.3.06-.61.06-.94 0-.32-.02-.64-.07-.94l2.03-1.58c.18-.14.23-.41.12-.61l-1.92-3.32c-.12-.22-.37-.29-.59-.22l-2.39.96c-.5-.38-1.03-.7-1.62-.94l-.36-2.54c-.04-.24-.24-.41-.48-.41h-3.84c-.24 0-.43.17-.47.41l-.36 2.54c-.59.24-1.13.57-1.62.94l-2.39-.96c-.22-.08-.47 0-.59.22L2.74 8.87c-.12.21-.08.47.12.61l2.03 1.58c-.05.3-.09.63-.09.94s.02.64.07.94l-2.03 1.58c-.18.14-.23.41-.12.61l1.92 3.32c.12.22.37.29.59.22l2.39-.96c.5.38 1.03.7 1.62.94l.36 2.54c.05.24.24.41.48.41h3.84c.24 0 .44-.17.47-.41l.36-2.54c.59-.24 1.13-.56 1.62-.94l2.39.96c.22.08.47 0 .59-.22l1.92-3.32c.12-.22.07-.47-.12-.61zM12 15.6c-1.98 0-3.6-1.62-3.6-3.6s1.62-3.6 3.6-3.6 3.6 1.62 3.6 3.6-1.62 3.6-3.6 3.6"}),"Settings");export{te as a,gs as b,pn as c,he as d,Kl as e,Ql as f,Zl as g,Jl as h,il as i,Lt as j,Et as k,Je as l,At as m,$s as n,dl as o,hl as p,ml as q,bl as r,Dt as s,sn as t,as as u,Ft as v,El as w,$o as x,Ya as y,qe as z,Cs as A,vn as B,Rs as C,Es as D,Ba as E,jn as F,Kc as G,Qc as H,ja as I,eu as J,Yo as K,xi as L,du as M,fu as N,pu as O,Ci as P,yu as Q,Ii as R,xu as S,Mu as T,wu as U,Ru as V,Eu as W,Oi as X,Au as Y,Iu as Z,_u as _,vu as $,Bu as aa,Gu as ba,Hi as ca,Du as da,Xu as ea,ku as fa,Vi as ga,cr as ha,ju as ia,pr as ja,Vu as ka,Uu as la,Ji as ma,Zu as na,Il as oa,_l as pa,ps as qa,vl as ra,Bl as sa,In as ta,Yl as ua,zl as va,nd as wa};
|
|
2
|
-
//# sourceMappingURL=chunk-BI3M6FZD.js.map
|