jbrowse-plugin-graphgenomeviewer 4.0.7 → 4.0.8

package/dist/chunks/chunk-VXUVWSY5.js

@@ -0,0 +1,2 @@
+import{c as Yn}from"./chunk-ZNGDWZ6Z.js";import{A as jn,E as Vn,F as Kn,a as Te,b as G,c as Ke,d as ce,e as ue,f as Wn,h as Ue,q as xe,u as Xn,v as zn,w as dn,x as fn}from"./chunk-NBDRTXH2.js";import{a as Ve,b as Fn,c as Hn}from"./chunk-F4C3BEPV.js";import{a as je}from"./chunk-XEJ47DLD.js";import{a as ze,b as me}from"./chunk-FY3BTKY4.js";var kn=ze((Ga,Un)=>{Un.exports=JBrowseExports["@mui/material/SvgIcon"]});var Nn=ze((ai,Nt)=>{Nt.exports=JBrowseExports["@jbrowse/core/util/colorBits"]});var is=ze((Yi,Pt)=>{Pt.exports=JBrowseExports.mobx});var ls=ze((Xi,Dt)=>{Dt.exports=JBrowseExports["@jbrowse/render-core/canvas2dUtils"]});var ee=me(kn(),1);function qa(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 Qn(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 mn(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 ot(e){let n=e[Math.floor(e.length/2)];return n?mn(n,.5):void 0}function rt(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 Fo(e,n,t){let r=mn(e,n),o=mn(e,t),s=(t-n)/3,a=Zn(e,n),l=Zn(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 Zn(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 st(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([Fo(r,a/s,(a+1)/s)])}return t}function at(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 Jn(e,n){let t=e.length-1;return n==="end"?{at:e[t],inward:e[t-1]}:{at:e[0],inward:e[1]}}var Ho=.5,Wo=80;function et(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 nt(e,n){return Math.max(0,n.min-e,e-n.max)}function tt(e,n,t){let r=nt(e[0].x,n),o=nt(e[e.length-1].x,n);return r===o?t:r<o?"start":"end"}function Yo(e,n){return{from:tt(e,et(n),"end"),to:tt(n,et(e),"start")}}var Xo=.75,zo=1.15,jo=.35,Vo=60;function Ko(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 g=m.x-e,b=m.y-n,w=g*-i+b*l;Math.abs(w)>Math.abs(u)&&(u=w);let E=g*l+b*i;f=Math.min(f,E),c=Math.max(c,E)}let d=o.length?c-f:0,h=Math.max(Math.abs(u)*zo/Xo,d*jo),p=s===void 0?h:Math.min(h,s);return u<0?-p:p}function ke({scaleX:e,scaleY:n}){return n/e}var Uo=3;function pn(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 it(e,n,t,r,o=Uo){let s=ke(r),a=pn(e[e.length-1],n[0],s)??pn(e[0],e.at(-1),s)??pn(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 hn(e,n){return e.map(t=>({x:t.x,y:t.y*n}))}function ko(e,n){for(let t of e)t.y0/=n,t.cy0/=n,t.cy1/=n,t.y1/=n;return e}function ve(e,n,t,r,o,s,a=[],l){let i=s.scaleX,u=ke(s),f=l??(s.pixelRows?Vo:void 0);if(u!==1)return ko(ve(hn(e,u),hn(n,u),t,r,o*u,{scaleX:i,scaleY:i},hn(a,u),f===void 0?void 0:f*u),u);let c=t?{from:"end",to:"start"}:Yo(e,n),d=Jn(e,c.from),h=Jn(n,c.to),p=d.at,m=h.at,g=p.x+r,b=p.y+o,w=m.x+r,E=m.y+o;if(t){let L=1,N=0;if(d.inward){let U=p.x-d.inward.x,V=p.y-d.inward.y,K=Math.hypot(U,V);K>0&&(L=U/K,N=V/K)}let R=Math.min(50,50*i),x=-N,A=L,v=(p.x+m.x)/2+r+x*R,_=(p.y+m.y)/2+o+A*R,j=g+L*R,W=b+N*R,D=w-L*R,q=E-N*R;return[{x0:g,y0:b,cx0:j,cy0:W,cx1:j+x*R,cy1:W+A*R,x1:v,y1:_},{x0:v,y0:_,cx0:D+x*R,cy0:q+A*R,cx1:D,cy1:q,x1:w,y1:E}]}else{let L=d.inward??{x:p.x-(m.x-p.x),y:p.y-(m.y-p.y)},N=h.inward??{x:m.x-(p.x-m.x),y:m.y-(p.y-m.y)},R=Math.hypot(w-g,E-b),x=(oe,Z,re,Se)=>{let se=Z.x-oe.x,ae=Z.y-oe.y,ie=Math.hypot(se,ae),he=ie===0?0:(re*se+Se*ae)/ie;return Math.max(0,Math.min(R*Ho,Wo,he*.5))},A=m.x-p.x,v=m.y-p.y,[_,j]=Qn(L.x,L.y,p.x,p.y,x(L,p,A,v)),[W,D]=Qn(N.x,N.y,m.x,m.y,x(N,m,-A,-v)),q=a.length?Ko(g,b,w,E,a,f):0,U=R===0?1:R,V=(w-g)/U,K=(E-b)/U,qe=-K*q,O=V*q,F=V*Math.abs(q),Q=K*Math.abs(q);return[{x0:g,y0:b,cx0:_+r+qe-F,cy0:j+o+O-Q,cx1:W+r+qe+F,cy1:D+o+O+Q,x1:w,y1:E}]}}var Qo=1;function gn(e,n){return n.flatMap(t=>e[t]??[])}function lt(e,n,t){let r=e[n.from],o=e[n.to];return r?.length&&o?.length?ve(r,o,n.from===n.to,0,0,t,gn(e,n.bypassed)):void 0}function Zo(e){let n=new Map;for(let t of e.nodes)if(G(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 Jo(e,n,t){let r=[];if(e)for(let o=Wn(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 er(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 Ya(e){let n=new Map(e.nodes.map(o=>[o.id,o])),t=Zo(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&&G(a)&&G(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>=Qo&&er(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:Jo(t.get(p),c,d)})}}}return r}function ja(e){let n=new Map,t=new Map,r=(d,h,p)=>{let m=d.get(h);m?m.push(p):d.set(h,[p])};for(let d of e.edges)r(n,d.from,d.to),r(t,d.to,d.from);let o=(d,h)=>d!==h&&n.get(d)?.length===1&&t.get(h)?.length===1,s=new Map,a=new Map;for(let d of e.nodes){let h=t.get(d.id);if(s.has(d.id)||h?.length===1&&o(h[0],d.id))continue;let p=[d.id];for(;;){let m=n.get(p.at(-1))?.[0];if(m===void 0||s.has(m)||!o(p.at(-1),m))break;p.push(m)}for(let m of p)s.set(m,d.id);p.length>1&&a.set(d.id,p)}for(let d of e.nodes)s.has(d.id)||s.set(d.id,d.id);let l=new Map(e.nodes.map(d=>[d.id,d])),i=[];for(let d of e.nodes){if(s.get(d.id)!==d.id)continue;let h=a.get(d.id);if(!h){i.push(d);continue}let p=h.map(b=>l.get(b)),m=p.filter(G),g=m.length?m.reduce((b,w)=>b.stable.start<=w.stable.start?b:w):p[0];i.push({id:d.id,name:d.name,length:p.reduce((b,w)=>b+w.length,0),depth:d.depth,stable:g.stable})}let u=new Set;for(let d of a.values())for(let h=1;h<d.length;h++)u.add(`${d[h-1]}>${d[h]}`);let f=new Set,c=[];for(let d of e.edges){let h=`${s.get(d.from)}>${s.get(d.to)}`;u.has(`${d.from}>${d.to}`)||f.has(h)||(f.add(h),c.push({from:s.get(d.from),to:s.get(d.to)}))}return{graph:{...e,nodes:i,edges:c,paths:void 0,pathVisits:void 0},runs:a}}function Va(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=yn(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]=ut(a,l,c,m),c=m})}return r}function Ka(e,n,t){if(e.length===0)return[];let r=yn(e),o=r.at(-1);return ut(e,r,n*o,t*o)}function ct(e){let n=yn(e);return bn(e,n,n.at(-1)/2)}function yn(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 bn(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 ut(e,n,t,r){let o=[bn(e,n,t)];for(let s=0;s<e.length;s++)n[s]>t&&n[s]<r&&o.push({...e[s]});return o.push(bn(e,n,r)),o}function xn(e){return e.length*6.2+8}function nr(e,n){return e.x0<n.x1&&e.x1>n.x0&&e.y0<n.y1&&e.y1>n.y0}function dt(e,n){return e.x1<=0||e.x0>=n.width||e.y1<=0||e.y0>=n.height}function tr(e,n){return e.x0>=0&&e.x1<=n.width&&e.y0>=0&&e.y1<=n.height}function Mn(e,n=[]){let t=[...n];return r=>dt(r,e)||t.some(o=>nr(o,r))?!1:(t.push(r),!0)}function Oe(e,n,t){let r=[];for(let s of e){let a=xn(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(dt(f,n))break;if(tr(f,n)&&t(f)){r.push({...s,y:u,w:a});break}}}return r}var or=26,rr=4e4,sr=8,ar=8,ir=5.7,lr=6,cr=13,Be=3,Qe=2;function en(e){return(e.length*ir+lr+Be)/2}var de=(cr+Be)/2,ur=en("1 bp");function Ze(e,n,t){return{x0:n-e,x1:n+e,y0:t-de,y1:t+de}}function mt(e,n){return{x0:6,x1:6+e.length*6.2+8+Be,y0:n-8-Be/2,y1:n+8+Be/2}}function Je(e,n,t){return t<n?e:Math.min(Math.max(e,n),t)}function gt(e){return e>=1e6?`${+(e/1e6).toFixed(1)} Mb`:e>=1e3?`${+(e/1e3).toFixed(1)} kb`:`${e} bp`}function ft(e){return`${gt(e)} del`}function dr(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 fr(e,n){return Math.abs(n.x)*(e.x1-e.x0)/2+Math.abs(n.y)*(e.y1-e.y0)/2}function pr(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 hr(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 mr(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 gr(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 ht=new WeakMap;function br(e,n,t,r){let{scaleX:o,scaleY:s}=t,a=ht.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=lt(e,i,t),f=u?ot(u):void 0;if(u&&f){let{minX:c,minY:d,maxX:h,maxY:p}=rt(u);l.push({deletion:i,curves:u,apex:f,extent:Math.max((h-c)*o,(p-d)*s)})}}return ht.set(e,{deletions:n,scaleX:o,scaleY:s,version:r,arcs:l}),l}function bt({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:g,curves:b,apex:w,extent:E}of br(e,t,o,u))if(E>=or){let L=ft(g.bp),N=en(L),R=Ze(N,0,0),x=w.x*f+s,A=w.y*c+a,v=E>=N*2/Qe,_=dr(b,w,f,c),j=fr(R,_)+ar,W=v?x:Je(x+_.x*j,N,l-N),D=v?A:Je(A+_.y*j,de,i-de),q=Math.hypot(W-x,D-A),U;if(q>0){let V={x:(W-x)/q,y:(D-A)/q},K=pr(R,V);U={arcX:x,arcY:A,labelX:W-V.x*K,labelY:D-V.y*K}}d.push({label:{key:`del:${g.edgeIndex}`,text:L,x:W,y:D,kind:"deletion",leader:U},box:Ze(N,W,D)})}let p=new Set,m=[...r??[]].sort((g,b)=>b.bp-g.bp);for(let g of m){let b=mr(g.nodeIds,e);if(!b)continue;let w=Math.max(b.minX*f+s,0),E=Math.min(b.maxX*f+s,l),L=Math.max(b.minY*c+a,0),N=Math.min(b.maxY*c+a,i);if(E<w||N<L)continue;let R=Math.max(E-w,N-L),x=ft(g.bp),A=en(x);if(R>=A*2/Qe){for(let j of g.nodeIds)p.add(j);let v=Je((w+E)/2,A,l-A),_=Je((L+N)/2,de,i-de);d.push({label:{key:`alleledel:${g.nodeIds.join(",")}`,text:x,x:v,y:_,kind:"deletion"},box:Ze(A,v,_)})}}for(let[g,b]of gr(e,n)){if(p.has(g))continue;let{x:w,y:E}=ct(b),L=E*c+a;if(L+de<=0||L-de>=i)continue;let N=hr(b,f,c);if(N<ur*2/Qe)continue;let R=gt(n.get(g)),x=en(R);if(N>=x*2/Qe){let A=w*f+s;h.push({label:{key:`node:${g}`,text:R,x:A,y:L,kind:"node"},box:Ze(x,A,L)})}}return{deletions:d,nodes:h}}function yt(e,n){return Math.max(sr,Math.floor(e*n/rr))}function wn(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 xr=3.4,xt=6,Mr=8,wr=18,Lr=60,Nr=8,Mt=15;function Ln(e,n){return[...e].sort((t,r)=>n(r)-n(t))}function ri(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:A})=>mt(x,A*a+o));e.legendSize.width>0&&u.push({x0:n-xt-e.legendSize.width,x1:n,y0:0,y1:xt+e.legendSize.height}),e.poppedFrom&&u.push({x0:0,x1:60+e.poppedFrom.label.length*6.2*1.2,y0:0,y1:44});let f=Mn(l,u),c=e.drawnRowLabels.length>0,d=c?Math.min(...e.drawnRowLabels.map(x=>x.y)):1/0,h=e.contigThickness*xr,p=Ln(e.bubbleHalos,x=>x.members),m=Oe(p.map(x=>{let{x:A,y:v}=i({x:x.labelAt.x,y:Math.min(x.labelAt.y,d)}),_=v-h/2-6;return{item:x,x:A,y:c?Math.max(_,Mt):_,text:x.label}}),l,f),g=Math.max(2,Math.min(Nr,Math.floor((o-Lr)/19))),b=Oe(Ln(e.bubbleGlyphs,x=>x.bubble.end-x.bubble.start).flatMap(x=>{let A=(x.bubble.start+x.bubble.end)/2*s+r,v=xn(x.label)/2;return A<0||A>n?[]:[{item:{glyph:x,glyphX:A},x:Math.min(Math.max(A,v),n-v),y:Mt,text:x.label,stack:g-1}]}),l,f),w=Oe(Ln(e.genePins,x=>x.gene.end-x.gene.start).map(x=>{let{x:A,y:v}=i(x.at);return{item:x,x:A,y:v+e.contigThickness+wr,text:x.covered<.98?`${x.gene.name} \u2026`:x.gene.name}}),l,f),E=e.nodePositions&&e.labelsNodeSizes?bt({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:[]},L=wn(E.deletions,f),N=Oe(p.flatMap(x=>x.routes.map(A=>{let{x:v,y:_}=i(A.at);return{item:{halo:x,route:A},x:v,y:_+4,text:A.text,stack:Mr}})),l,f),R=wn(E.nodes,f,yt(n,t));return{bubbles:m,glyphs:b,genes:w,routes:N,sizes:[...L,...R]}}var Y=me(Nn(),1);var P=me(Nn(),1);function _e(e,n){return(0,P.packAbgr)((0,P.abgrRed)(e),(0,P.abgrGreen)(e),(0,P.abgrBlue)(e),Math.round((0,P.abgrAlpha)(e)*n))}function ii(e,n){return(0,P.packAbgr)(Math.min(255,Math.round((0,P.abgrRed)(e)*n)),Math.min(255,Math.round((0,P.abgrGreen)(e)*n)),Math.min(255,Math.round((0,P.abgrBlue)(e)*n)),(0,P.abgrAlpha)(e))}var Sr=4;function Ge(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 Er({nodeId:e,nodeById:n,neighbors:t}){let r=new Set([e]),o=[e],s=1/0,a=-1/0;for(let l=0;l<Sr&&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&&G(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 Rr({nodeId:e,nodeById:n,neighbors:t}){let r=n.get(e);return r===void 0?void 0:G(r)?{start:r.stable.start,end:r.stable.start+r.length}:Er({nodeId:e,nodeById:n,neighbors:t})}function St(e){let n=new Map(e.nodes.map(o=>[o.id,o])),t=Ge(e),r=new Map;for(let o of e.nodes){let s=Rr({nodeId:o.id,nodeById:n,neighbors:t});s&&r.set(o.id,(s.start+s.end)/2)}return r}var Et=new WeakMap;function Rt(e,n,t,r,o){let{scaleX:s,scaleY:a}=t,l=Et.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,ve(c,d,f.from===f.to,0,0,t,h?gn(e,h):[]))}}return Et.set(e,{graph:n,scaleX:s,scaleY:a,deletions:r,version:o,curves:i}),i}var Ar=(0,Y.packAbgr)(119,119,119,217),At=(0,Y.packAbgr)(136,136,136,217),Ir=(0,Y.packAbgr)(24,24,28,240),Cr=2.2,Tr=(0,Y.packAbgr)(30,30,36,245),vr=1.8,nn=.18,Or=11,Br=1.2,_r=1e-6,Gr=12,$r=.45;function Pr(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,g=f/d;s=(h+m)/2,a=(p+g)/2;let b=h*s+p*a;b>.1&&(s/=b,a/=b)}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 Dr(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 tn(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 on(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 qr=[[68,1,84],[59,82,139],[33,145,140],[94,201,98],[253,231,37]],Fr=[[237,137,44],[158,42,122]],Hr=[[220,50,50],[50,120,220]];function Sn(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 Wr(e,n){let[t,r,o]=on(Xn(e),dn,fn);return tn(t,r,o,n)}var Yr=300,Xr=.7,zr=.5,jr=(0,Y.packAbgr)(60,65,72,255);function Mi(e,n){let t=St(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 Vr(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 Kr(e,n,t,r){switch(t){case"random":return Wr(e.id,1);case"depth":return Sn(qr,r.maxDepth>r.minDepth?(e.depth-r.minDepth)/(r.maxDepth-r.minDepth):.5);case"node-length":return Sn(Hr,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):Sn(Fr,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 jr;let a=Math.max(0,Math.min(1,(s-o.start)/o.span)),[l,i,u]=on(a*Yr,Xr,zr);return tn(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]=on(o,.75,.5);return tn(s,a,l,1)}default:return(0,Y.packAbgr)(52,152,219,255)}}function Ur(e,n=1){let t=e.x1-e.cx1,r=(e.y1-e.cy1)*n;return Math.hypot(t,r)>0?Math.atan2(r,t):Math.atan2((e.y1-e.y0)*n,e.x1-e.x0)}function kr(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 Qr(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 It=new WeakMap;function Zr(e){let n=It.get(e);if(!n){let t=[],r=new Map,o=e.paths??[];for(let s=0;s<o.length;s++){let[a,l,i]=on(zn(s,o.length),dn,fn);t.push(tn(a,l,i,.85)),r.has(o[s].name)||r.set(o[s].name,s)}n={pathColorByIndex:t,pathIndexByName:r},It.set(e,n)}return n}var Ct=new WeakMap,Jr=new Map;function es(e){let n=Ct.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])}Ct.set(e,n)}return n}function wi(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:g=0}=e,b=i==="depth"?Vn(t):0,w=f.scaleX,E=ke(f),L=Rt(n,t,f,p,g),N=[],R=new Map,x=[],A=new Map,v=[],_=new Map,j={...Vr(t),referenceRamp:h},W=new Map;if(o==="rainbow")for(let O=0;O<t.nodes.length;O++)W.set(t.nodes[O].id,O);let D=t.paths?.length??0,{pathColorByIndex:q,pathIndexByName:U}=Zr(t),V=O=>q[U.get(O)??-1]??At,K=l&&D>1?es(t):Jr,qe=w>(c?1:$r);for(let O=0;O<t.edges.length;O++){let F=t.edges[O],Q=n[F.from],oe=n[F.to],Z=L.get(O);if(!Q?.length||!oe?.length||!Z||m?.has(O))continue;let re=F.pathIds?.length??0,se=p?.get(O)!==void 0,ae=u?.edgeIndexes.has(O)??!1,ie=a/2*(se?Cr:1)*(ae?vr:1),he=l&&re>0,Ee=se?Ir:Ar,ye=u?ae?Tr:_e(Ee,nn):Ee,Fe=(Ye,Xe,J,Ie)=>{let le=Ye===0&&Xe===0?Z:at(Z,Ye,Xe);if(!(d&&!Qr(le,d))){if(se)for(let Ce of st(le,Or/w))v.push({curves:Ce,thickness:ie,color:J});else v.push({curves:le,thickness:ie,color:J});if(qe&&Ie!==void 0){let Ce=le[le.length-1];x.push({x:Ce.x1,y:Ce.y1,angle:Ur(Ce,E),size:Gr,color:Ie})}}},He=v.length,Re=x.length;if(!he)Fe(0,0,ye,ye);else{let Ye=it(Q,oe,re,f),Xe=Math.floor((re-1)/2);for(let J=0;J<re;J++){let Ie=Ye[J],le=V(F.pathIds[J]);Fe(Ie.x,Ie.y,u&&F.pathIds[J]!==u.name?_e(le,nn):le,J===Xe?ye:void 0)}}let We=v.length-He;We>0&&_.set(O,{start:He,count:We});let Ae=x.length-Re;Ae>0&&A.set(O,{start:Re,count:Ae})}for(let[O,F]of Object.entries(n)){let Q=r.get(O);if(!Q||d&&!kr(F,d))continue;let oe=Kr(Q,W.get(O)??0,o,j),Z=u!==void 0&&!u.nodeIds.has(O),re=Z?_e(oe,nn):oe,Se=Kn(Q,s,i,b),se=Se/2,ae=N.length,ie=K.get(O),he=Se/D,Ee=F.length>=2;if(Ee&&ie?.length&&he>=Br){let Fe=Pr(F,E),He=1/Math.max(w,_r);for(let Re of ie){let We=(Re-(D-1)/2)*he*He,Ae=q[Re]??At;N.push({points:Dr(F,Fe,We,E),thickness:he/2,color:Z?_e(Ae,nn):Ae})}}else Ee&&N.push({points:F,thickness:se,color:re});let ye=N.length-ae;ye>0&&R.set(O,{start:ae,count:ye})}return{nodeStrokes:N,nodeStrokeRuns:R,arrows:x,arrowRuns:A,edgeCurves:v,edgeCurveRuns:_}}function Ei(e){return Te(e)?Tt(e):void 0}function Tt(e){return{sample:Ve(e.stable.refName),haplotype:Fn(e.stable.refName),refName:Hn(e.stable.refName),start:e.stable.start,end:e.stable.start+e.length}}function ns(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 ts(e){return`${e.haplotype??e.sample}\0${e.refName}`}function Ri(e,{maxGap:n=1/0}={}){let t=new Map;for(let o of e.nodes)if(Te(o)){let s=Tt(o),a=ts(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=ns(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 os(e,n){return(n.haplotype===void 0?void 0:e(n.haplotype))??e(n.sample)}function Ai(e,n){return e.flatMap(t=>{let r=os(n,t);return r===void 0?[]:[{...t,sample:r}]})}function Ii(e){return`${e.refName}:${e.start+1}-${e.end}`}function En(e){return`${e.refName}:${(e.start+1).toLocaleString()}-${e.end.toLocaleString()}`}var vt=me(je(),1),Ot=(0,ee.createSvgIcon)((0,vt.jsx)("path",{d:"M9.01 14H2v2h7.01v3L13 15l-3.99-4zm5.98-1v-3H22V8h-7.01V5L11 9z"}),"CompareArrows");var Bt=me(je(),1),_t=(0,ee.createSvgIcon)((0,Bt.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 Gt=me(je(),1),Rn=(0,ee.createSvgIcon)((0,Gt.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 rs="Linear genome view",ss="Linear synteny view",as="No synteny track in this session aligns these assemblies";function $t({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 Fi({contributors:e,syntenyTracks:n,onShowLinear:t,onShowSynteny:r}){let o=n[0]?.coverage??e.length;return[...$t({label:rs,icon:Rn,entries:e,entryLabel:s=>`${s.sample} ${En(s)}`,onSelect:s=>()=>{t({location:s,assembly:s.sample})}}),...e.length>=2?$t({label:`${ss} (${o} assemblies)`,icon:Ot,entries:n,entryLabel:s=>s.coverage<o?`${s.name} (${s.coverage} of ${o})`:s.name,nameTheOnlyEntry:!1,disabledHelpText:as,onSelect:s=>()=>{r(s.trackId)}}):[]]}function Hi({own:e,reference:n,highlight:t,onShowLinear:r,onHighlight:o}){let s=[e?{label:`Open in ${e.assembly} \u2014 ${En(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:_t,onClick:()=>{o(t)}}]:[],...s.map(({label:a,target:l})=>({label:a,icon:Rn,onClick:()=>{r(l)}}))]}var cs=10;function us(e,n){let t=e.end-e.start,r=e.longestAlleleLength-e.shortestAlleleLength;return n.find(o=>r<=cs*(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 ds=40;function ge(e){return e>=1e4?`${(e/1e3).toFixed(0)} kb`:e>=1e3?`${(e/1e3).toFixed(1)} kb`:`${e} bp`}function fs(e,n){let t=ge(e),r=ge(n),o=r.slice(r.indexOf(" "));return t.endsWith(o)?`${t.slice(0,-o.length)}\u2013${r}`:`${t}\u2013${r}`}function qt(e){return e>=Yn?"\u22652.1B routes":e>=1e5?`${(e/1e3).toFixed(0)}k routes`:`${e} routes`}function ji(e,n=[]){let t=ps(e,n);return e.partial?{...t,label:`${t.label}, partial`}:t}function ps(e,n){let t=e.end-e.start,{shortestAlleleLength:r,longestAlleleLength:o,pathCount:s,segmentCount:a,inversion:l}=e,i=fs(r,o);if(a>=ds)return{kind:"superbubble",label:`${i} superbubble, ${a} segments, ${qt(s)}${l?", inv":""}`};if(l)return{kind:"inversion",label:`${ge(o)} inv`};let u=o!==r?us(e,n):void 0;if(u){let c=u.name?` (${u.name})`:"";return{kind:"repeat",label:`${i} repeat array${c}, ${qt(s)}`}}let f=s>2?`, ${s} alleles`:"";return r===o&&r===t?t<=1?{kind:"snp",label:"SNP"}:{kind:"substitution",label:`${ge(t)} sub`}:t===0||r===t?{kind:"insertion",label:`\u2264${ge(o-t)} ins${f}`}:o===t?{kind:"deletion",label:`${ge(t-r)} del${f}`}:{kind:"complex",label:`${ge(t)} ref \u2192 ${i}, ${s} alleles`}}var Vi={snp:"#5b6b7a",substitution:"#5b6b7a",insertion:"#2f8fd6",deletion:"#c94040",inversion:"#e07b00",repeat:"#8e3fbf",superbubble:"#8e3fbf",complex:"#8e3fbf"},Ki={snp:"SNP",substitution:"substitution",insertion:"insertion",deletion:"deletion",inversion:"inversion",repeat:"repeat array",superbubble:"superbubble",complex:"site"};function Ui(e){return e.segments.split(",").filter(Boolean)}function rn(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 hs(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 ms(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 gs(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=gs(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 Ht(e){let n=new Map(e.nodes.filter(Te).map(s=>[s.id,s])),t=Ge(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}=hs(s.id,n,t,r),i=ms(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 ys=1,xs=.015;function sn({graph:e,span:n,rowY:t,positions:r}){let o=n*xs,s=Ht(e),a=new Map(e.nodes.map(f=>[f.id,f])),l=[];for(let f of s){let c=f.refSpan-f.pathLength;c>=ys&&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 g=a.get(p);if(g){let b=f.nodeOffsets[m],w=d/f.nodeIds.length,E=N=>f.start+N*h,L=t(g);r[p]=h?[{x:E(b),y:L},{x:E(b+g.length),y:L}]:[{x:f.start+m*w,y:L},{x:f.start+(m+1)*w,y:L}]}})}let i=Ge(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 Ms(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 Wt(e,n){let t=ce(e);if(t.length===0)return;let r=Ue(t,n),o=Ms(e),s=ue(t),a=sn({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 ws=12,An=20;function Yt(e){return 10+8*Math.log2(1+e.length)}function fe(e,n){return e<n?-1:e>n?1:0}function zt(e,n){return fe(e.at,n.at)||fe(e.side,n.side)||fe(e.hops,n.hops)}var Xt={at:1/0,side:0,hops:0};function Ls(e){let n=new Map(e.nodes.map(r=>[r.id,r])),t=Es(e);return e.nodes.map(r=>r.id).sort((r,o)=>zt(t.get(r)??Xt,t.get(o)??Xt)||fe(n.get(r).name,n.get(o).name)||fe(r,o))}function Ns(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 b of f.nodeIds){let w=r.get(b);w!==void 0&&(h!==void 0&&w!==h&&(w>h?c++:d++),h=w)}let p=d>c?[...f.nodeIds].reverse():f.nodeIds,m,g=new Set;for(let b of p)if(o.has(b)){if(m===void 0)for(let w of g)i(b,w);m=b}else t.has(b)&&(m===void 0?g.add(b):(l(m,b),m=b))}let u=[];for(let f=a;f!==void 0;f=o.get(f))u.push(f);return u}function Ss(e,n){let t=[...n.entries()].sort(([a,l],[i,u])=>zt(l,u)||fe(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 Es(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(G).sort((s,a)=>fe(s.stable.start,a.stable.start)||fe(s.name,a.name)),o=new Map;return Ns(e,r).forEach((s,a)=>{o.set(s,{at:a,side:0,hops:0})}),Ss(t,o),o}function Rs(e){let n=new Map(e.nodes.map(c=>[c.id,c])),t=Ls(e),r=new Map(t.map((c,d)=>[c,d])),o=c=>G(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 b of h)p=Math.min(p,u.get(b));let m=u.get(d),g=Math.max(m,p-1);u.set(d,Math.floor((m+g)/2))}let f=[];for(let c of t)(f[u.get(c)]??=[]).push(c);return{layers:f,layerOf:u,preds:s}}function jt(e){if(!e.nodes.some(G))return;let n=new Map(e.nodes.map(c=>[c.id,c])),t=c=>G(n.get(c)),{layers:r,layerOf:o,preds:s}=Rs(e),a=r.map(c=>{let d=0;for(let h of c)d=Math.max(d,Yt(n.get(h)));return d}),l=[],i=0;for(let c=0;c<r.length;c++)l[c]=i,i+=a[c]+ws;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(b=>u.has(b)),g=0;for(let b of m)g+=u.get(b)/m.length;d.push({id:p,ideal:g===0?An:g})}d.sort((p,m)=>p.ideal-m.ideal);let h=new Set([0]);for(let{id:p,ideal:m}of d){let g=Math.round(m/An);g===0&&(g=m>=0?1:-1);for(let b=1;h.has(g);b++)h.has(g+b)?h.has(g-b)||(g-=b):g+=b;h.add(g),u.set(p,g*An)}}let f={};for(let c of e.nodes){let d=o.get(c.id),h=Yt(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 As(e,n){let t=new Map;for(let o of e.nodes)if(Ke(o)){let{sample:s}=rn(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 Vt(e,n,t=[]){let r=ce(e),o=As(e,t);if(r.length===0||o.length===0)return;let s=Ue(r,n),a=new Map(o.map((d,h)=>[d,h+1])),l=ue(r),i=o.length+1,u=sn({graph:e,span:s,rowY:d=>{let h=d.stable;return((h?a.get(rn(h.refName).sample):void 0)??i)*20},positions:l}),c=[{label:rn(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 In(){return{rectangles:[],curves:[],corners:[],verticalRectangles:[],featureRectangles:[]}}var De=!1,Kt=["#1f77b4","#ff7f0e","#2ca02c","#d62728","#9467bd","#8c564b","#e377c2","#7f7f7f"];function uo(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=>Kt[n.id%Kt.length]),trackAlpha:e.trackAlpha??(()=>1)}}var T=uo({}),Bn=null,y=[],M=[],S=[],C=new Map,Le=[],I=[],be=[],Gn=[],$n=[],k=-1,$=In(),_n={minX:0,maxX:0,minY:0,maxY:0},$e,an=new Map;function fo(e,n,t=[],r={}){if(T=uo(r),Bn=r.bed??null,$=In(),be=[],Gn=[],$n=[],Le=[],_n={minX:0,maxX:0,minY:0,maxY:0},$e=void 0,an=new Map,e.length===0||n.length===0)return;let o=Cn(e);y=[],o.forEach((i,u)=>{y[u+1]=i}),M=Cn(n);let s=t.length>0;S=Cn(ha(t));for(let i=M.length-1;i>=0;i-=1){let u=M[i];if(u.type===void 0&&(u.type="haplotype"),u.hidden===!0){M.splice(i,1);continue}u.indexOfFirstBase!==void 0&&($e=u.name)}if(M.length===0)return;Cs(0),C=Qt(),Me(M),T.showReads&&S.length>0&&Me(S),co(),T.mergeNodesFlag&&(Zt(),vn(),T.showReads&&S.length>0&&Tn(),pa(),C=Qt(),co(),Me(M),T.showReads&&S.length>0&&Me(S)),Zt(),ks(),De&&console.log(`${y.length} nodes.`),vn(),k=eo(),Qs(),vn(),k=eo();let a=T.showReads&&T.coarsenedReadView&&!s;if(a){let i=$e===void 0?-1:M.findIndex(c=>c.name===$e),u=i===-1?0:i,f=M.filter((c,d)=>d!==u);if(f.length>0){let c=so(f,"haplotype");c.length>0&&(M=[M[u]],S=c)}}ia(),na(),T.showExonsFlag&&Bn!==null&&aa();let l=T.showReads&&T.coarsenedReadView&&S.length>0&&!a;return T.showReads&&S.length>0?(Ts(),Tn(),Me(S),l&&(S=so(S,"read"),Tn(),Me(S)),Os(),M=M.concat(S)):y.forEach(i=>{i&&(i.incomingReads=[],i.outgoingReads=[],i.internalReads=[])}),Js(),ua(),De&&(console.log("Tracks:"),console.log(M),console.log("Nodes:"),console.log(y),console.log("Lane assignment:"),console.log(be)),Xs(),{nodes:y,tracks:M,reads:S,nodeMap:C,shapes:$,bounds:_n,maxOrder:k,trackForRuler:$e,coarsenedEdgeMeta:an}}function we(e,n="plain"){return T.trackColor(e,n)}function Is(e){return T.trackAlpha(e)}function Cn(e){return structuredClone(e)}function X(e){let n=String(e);return n.length>=1&&n.startsWith("-")}function B(e){return X(e)?String(e).substring(1):e}function pe(e){return X(e)?e:`-${e}`}function po(e){return X(e)?B(e):pe(e)}function ln(e){return(e||1/e)>=0}function Cs(e){let n=new Set,t=M[e].sequence;for(let r=0;r<t.length;r+=1){let o=t[r];if(X(o)){let s=B(o),a=t.indexOf(s);(a===-1||a>r)&&n.add(s)}}for(let r of M){t=r.sequence;for(let o=0;o<t.length;o+=1){let s=t[o];X(s)?n.has(B(s))&&(t[o]=B(s)):n.has(s)&&(t[o]=pe(s))}}y.forEach(r=>{r&&n.has(r.name)&&(r.seq=r.seq.split("").reverse().join(""))})}function Ts(){Le=[];for(let n=0;n<=k;n+=1)Le[n]=[];let e=new Map;y.forEach(n=>{n.y!==void 0&&Ut(e,n.order,n.y+n.contentHeight)}),M.forEach(n=>{n.type==="haplotype"&&n.path.forEach(t=>{Ut(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,Le[n.order].push(t))})}function Ut(e,n,t){let r=e.get(n);r===void 0?e.set(n,t):e.set(n,Math.max(r,t))}var cn=7;function vs(){y.forEach(e=>{e&&(e.incomingReads=[],e.outgoingReads=[],e.internalReads=[])}),S.forEach((e,n)=>{if((e.width===void 0||e.width===0)&&(e.width=cn),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 Os(){Fs(),vs();let e=y.slice();e.sort(un);let n=new Map;for(let a=0;a<S.length;a++){let l=S[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?cn:0;e.forEach(a=>{if(a)for(let l of t)Bs(n.get(l),a,r)});let o=qs(),s=[];S.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],g=c[p],b=!!(g&&m&&g.isForward!==void 0&&m.isForward!==void 0&&(m.order>g.order||g.order===m.order&&g.isForward===m.isForward));h.betweenCycleReverseTraversal=b,s.push({readIndex:l,pathIndex:p,previousY:m?.y,previousNode:m?.node})}})}),s.sort(Gs),s.forEach(a=>{let l=S[a.readIndex],i=l.path[a.pathIndex];i.y=o[i.order],o[i.order]+=l.width}),De&&(console.log("Reads:"),console.log(S))}function Bs(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,g]=p,b=S[m].path,w=g+1,E;for(let L=1;L<=g;L++){let N=b[g-L]?.y;if(N!==void 0){w=L,E=N;break}}return{entry:p,decisionStep:w,y:E}});i.sort(_s),o=i.map(p=>p.entry);let u=l,f=new Map;o.forEach(p=>{let m=S[p[0]];m.path[p[1]].y=u,kt(f,m,p[1],u,n),u+=m.width});let c=u;s.sort($s);let d=new Map;u=l,s.forEach(p=>{let m=S[p[0]],g=m.firstNodeOffset??0;m.path[p[1]].y=u,d.set(u,g);let b=f.get(u);b===void 0||b+1<=g?(u+=m.width,c=Math.max(c,u)):(f.set(u,0),o.forEach(w=>{let E=S[w[0]],L=E.path[w[1]];L.y!==void 0&&L.y>=u&&(L.y+=m.width,kt(f,E,w[1],L.y,n))}),u+=m.width,c+=m.width)}),a.sort(Ds),a.forEach(p=>{let m=S[p],g=m.firstNodeOffset??0,b=m.finalNodeCoverLength??0;for(u=l;g<(f.get(u)??-1/0)+2||b>(d.get(u)??1/0)-3;)u+=m.width||cn;m.path[0].y=u,f.set(u,b),c=Math.max(c,u)});let h=c-n.y-n.contentHeight;n.contentHeight+=h,sa(n,h)}function _s(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 kt(e,n,t,r,o){t===n.path.length-1?e.set(r,n.finalNodeCoverLength??0):e.set(r,o.sequenceLength)}function Gs(e,n){let t=S[e.readIndex],r=S[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 $s([e,n],[t,r]){let o=S[e],s=S[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?un(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 Ps(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 Ds(e,n){let t=S[e],r=S[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 qs(){let e=[];for(let n=0;n<=k;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))}),M.forEach(n=>{n.path.forEach(t=>{e[t.order]=Math.max(e[t.order],(t.y??0)+n.width)})}),e}function Fs(){S.forEach(e=>{let n=Math.abs(e.indexSequence[0]),t=ln(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=ln(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 Tn(){S.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]=B(e.sequence[a]);let t=e.sequenceNew??[];t.reverse(),e.sequenceNew=t;for(let a of t)a.nodeName=B(a.nodeName),Hs(a.mismatches,ne(a.nodeName).sequenceLength);let r=e.firstNodeOffset??0,o=ne(e.sequence[0]).sequenceLength;e.firstNodeOffset=o-(e.finalNodeCoverLength??0);let s=ne(e.sequence[e.sequence.length-1]).sequenceLength;e.finalNodeCoverLength=s-r}})}function Hs(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=ho(t.seq).split("").reverse().join(""))})}var Ws={A:"T",T:"A",C:"G",G:"C"};function ho(e){let n="";for(let t=e.length-1;t>=0;t-=1)n+=Ws[e.charAt(t)]??"N";return n}function Me(e){e.forEach(n=>{n.indexSequence=[],n.sequence.forEach(t=>{let r=t,o=C.get(B(r)),s=y[o]?.switched??!1;s&&(r=po(r));let a=C.get(B(r));if(a===0)throw new Error("Node "+B(r)+" has prohibited index 0");X(r)!==s?n.indexSequence.push(-a):n.indexSequence.push(a)})})}var Ys=20;function Xs(){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+Ys+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)))}),M.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)})}),_n=e}function ne(e){let n=C.get(B(e)),t=n===void 0?void 0:y[n];if(t===void 0)throw new Error(`Unknown node ${e}`);return t}function Qt(){return C=new Map,y.forEach((e,n)=>{e&&C.set(e.name,n)}),C}function Zt(){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)}};M.forEach(t),T.showReads&&S.length>0&&S.forEach(t),y.forEach((r,o)=>{r&&(r.successors=Array.from(e[o]),r.predecessors=Array.from(n[o]))})}function Jt(e){let n=0,t=0,r=0;e.forEach(o=>{let s=Math.abs(o);if(o<0){I[s]===void 0&&(I[s]=t);let a=I[s];a<r&&(r=a),n=a,t=a-1}else{I[s]===void 0&&(I[s]=n);let a=I[s];n=a+1,t=a}}),r<0&&Pn(-r)}function zs(e){let n=0,t=0;for(;n<e.length&&I[Math.abs(e[n])]===void 0;)n+=1;if(n>=e.length)return null;let r=e[n],o,s;r>=0?(o=I[r]-1,s=-1):(o=I[-r]+1,s=1);for(let a=n-1;a>=0;a-=1){let l=Math.abs(e[a]);I[l]===void 0&&(I[l]=o,t=Math.min(t,o),o+=s)}return t<0&&Pn(-t),n}var js=-1;function Vs(e){let n=[];for(let t of e)n.push(t??js);return n}function vn(){let e,n,t,r,o=0,s=T.showReads&&S.length>0?M.concat(S):M,a={stamp:new Int32Array(y.length),generation:0};I=new Array(y.length).fill(void 0),y.forEach(i=>{i&&(i.order=void 0)}),Jt(M[0].indexSequence);for(let i=1;i<s.length;i+=1){if(De&&console.log(`generating order for track ${i+1}`),t=zs(s[i].indexSequence),t===null){s[i].type==="haplotype"?Jt(s[i].indexSequence):(s.splice(i,1),S.splice(i-M.length,1),i-=1);continue}for(e=Us(s[i].indexSequence);t<e.length;){for(r=t,t+=1;t<e.length&&I[e[t]]===void 0;)t+=1;if(t<e.length){n=I[e[r]]+1;for(let u=r+1;u<t;u+=1)I[e[u]]=n,n+=1;if(I[e[t]]>I[e[r]])I[e[t]]<n&&no(e[t],e[t-1],n);else if(s[i].indexSequence[t]>=0&&!Ks(e[t],e[r],a))no(e[t],e[t-1],n);else if(s[i].indexSequence[r]<0||y[e[r+1]].degree<2&&I[e[t]]<I[e[r]]){n=I[e[r]]-1;for(let u=r+1;u<t;u+=1)I[e[u]]=n,o=Math.min(o,n),n-=1}}else if(s[i].indexSequence[r]>=0){n=I[e[r]]+1;for(let u=r+1;u<e.length;u+=1){let f=e[u];I[f]===void 0&&(I[f]=n,n+=1)}}else{n=I[e[r]]-1;for(let u=r+1;u<e.length;u+=1){let f=e[u];I[f]===void 0&&(I[f]=n,o=Math.min(o,n),n-=1)}}}}o<0&&Pn(-o);let l=Vs(I);I=l,l.forEach((i,u)=>{let f=y[u];f!==void 0&&(f.order=i)})}function Ks(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 eo(){let e=-1;return I.forEach(n=>{n!==void 0&&n>e&&(e=n)}),e}function Us(e){return e.map(n=>Math.abs(n))}function Pn(e){I.forEach((n,t)=>{n!==void 0&&(I[t]=n+e)})}function no(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=I[a];i!==void 0&&i<l&&(!r.has(a)||r.get(a)<l)&&(r.set(a,l),y[a].successors.forEach(u=>{I[u]>i&&u!==n&&o.push([u,l+1])}),a!==e&&y[a].predecessors.forEach(u=>{I[u]>i&&u!==n&&o.push([u,l+1])}))}r.forEach((a,l)=>{I[l]=a})}function ks(){y.forEach(e=>{e.tracks=[]}),M.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 Qs(){let e=M[0],n=M.slice(1,M.length);if(T.showReads&&S.length>0&&(n=n.concat(S)),Zs(n,e),T.showReads&&S.length>0){for(let t=0;t<S.length;t++)if(S[t]!==n[t+M.length-1])throw new Error("Read inequality")}}function Zs(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=B(s[a]),i=ne(l);if(r&&!r.has(l)){let u=a>0?ne(s[a-1]).order:void 0,f=a<s.length-1?ne(s[a+1]).order:void 0;(u===void 0||u<i.order)&&(f===void 0||i.order<f)&&to(t,l,X(s[a])?1:-1),(u===void 0||u>i.order)&&(f===void 0||i.order>f)&&to(t,l,X(s[a])?-1:1)}}}for(let o of e)o.sequence.forEach((s,a)=>{let l=t.get(B(s));l!==void 0&&l>0&&(o.sequence[a]=po(s),o.indexSequence[a]=-o.indexSequence[a])});t.forEach((o,s)=>{if(o>0){let a=ne(s);a.seq=ho(a.seq),a.switched=!0}})}function to(e,n,t){e.set(n,(e.get(n)??0)+t)}function Js(){let e=0,n=20,t=-1,r=y.slice();r.sort(un);let o=ea();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 ea(){let e=[],n=[],t=[],r=[];for(let o=0;o<=k;o+=1)e.push(0),n.push(0),t.push(0);M.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<=k;o+=1)r.push(Math.max(Math.max(0,e[o]-1)+Math.max(0,n[o-1]-1),t[o]));return r}function na(){let e,n,t,r,o,s,a=[],l=[];for(let i=0;i<=k;i+=1){be[i]=[],l[i]=new Map,a[i]=[];for(let u=0;u<M.length;u+=1)a[i][u]=null}M.forEach((i,u)=>{n=Math.abs(i.indexSequence[0]),t=ln(i.indexSequence[0]),r=y[n],i.path=[],i.path.push({order:r.order,lane:null,isForward:t,node:n}),H(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=ln(i.indexSequence[f]),r=y[n],r.order>o.order){s||(i.path.push({order:o.order,lane:null,isForward:!0,node:null}),H(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}),H(c,null,u,e,a,l),e+=1;t?(i.path.push({order:r.order,lane:null,isForward:!0,node:n}),H(r.order,n,u,e,a,l),e+=1):(i.path.push({order:r.order,lane:null,isForward:!0,node:null}),H(r.order,null,u,e,a,l),e+=1,i.path.push({order:r.order,lane:null,isForward:!1,node:n}),H(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}),H(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}),H(c,null,u,e,a,l),e+=1;t?(i.path.push({order:r.order,lane:null,isForward:!1,node:null}),H(r.order,null,u,e,a,l),e+=1,i.path.push({order:r.order,lane:null,isForward:!0,node:n}),H(r.order,n,u,e,a,l),e+=1):(i.path.push({order:r.order,lane:null,isForward:!1,node:n}),H(r.order,n,u,e,a,l),e+=1)}else t!==s?(i.path.push({order:r.order,lane:null,isForward:t,node:n}),H(r.order,n,u,e,a,l),e+=1):(i.path.push({order:r.order,lane:null,isForward:!t,node:null}),H(r.order,null,u,e,a,l),e+=1,i.path.push({order:r.order,lane:null,isForward:t,node:n}),H(r.order,n,u,e,a,l),e+=1)});for(let i=0;i<=k;i+=1)oa(be[i],i)}function H(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]};be[e].push(i),n!==null&&s[e].set(n,i)}else l.type="multiple",l.tracks.push(a);o[e][t]=a}function ta(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(M[o.trackID].path[o.segmentID-1].order===n-1)o.idealLane=M[o.trackID].path[o.segmentID-1].lane??void 0,o.idealY=M[o.trackID].path[o.segmentID-1].y;else if(o.segmentID<M[o.trackID].path.length-1&&M[o.trackID].path[o.segmentID+1].order===n-1)o.idealLane=M[o.trackID].path[o.segmentID+1].lane??void 0,o.idealY=M[o.trackID].path[o.segmentID+1].y;else{for(t=o.segmentID-1;t>=0&&M[o.trackID].path[t].order!==n-1;)t-=1;t<0?(o.idealLane=o.trackID,o.idealY=null):(o.idealLane=M[o.trackID].path[t].lane??void 0,o.idealY=M[o.trackID].path[t].y)}r.idealLane+=o.idealLane}),r.idealLane/=r.tracks.length})}function oa(e,n){let t=0,r=new Set,o=20,s=!1,a=-1;ta(e,n),e.sort(ro),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(ro),l.tracks.forEach(i=>{i.lane=t,i.trackID===a&&l.node===null&&s&&(o+=10),M[i.trackID].path[i.segmentID].lane=t,M[i.trackID].path[i.segmentID].y=o,i.idealY!=null&&r.add(i.idealY-o),t+=1,o+=M[i.trackID].width,l.node!==null&&(y[l.node].contentHeight+=M[i.trackID].width),a=i.trackID}),o+=25}),ra(e,r)}function ra(e,n){let t=0,r=Number.MAX_SAFE_INTEGER;n.forEach(o=>{oo(e,o)<r&&(r=oo(e,o),t=o)}),e.forEach(o=>{o.node!==null&&(y[o.node].y+=t),o.tracks.forEach(s=>{let a=M[s.trackID].path[s.segmentID];a.y=a.y+t})})}function sa(e,n){e.order<0||be[e.order]===void 0||(be[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=M[o.trackID].path[o.segmentID];s.y=s.y+n}))}else t.tracks.forEach(r=>{let o=M[r.trackID].path[r.segmentID];o.y>=e.y&&(o.y=o.y+n)})}),Le[e.order].length>0&&Le[e.order].forEach(t=>{y[t]!==e&&y[t].y>e.y&&(y[t].y+=n)}))}function oo(e,n){let t=0;return e.forEach(r=>{r.tracks.forEach(o=>{o.idealY!=null&&M[o.trackID].type!=="read"&&(t+=Math.abs(o.idealY-n-M[o.trackID].path[o.segmentID].y)*M[o.trackID].width)})}),t}function ro(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 un(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 aa(){let e,n,t={};Bn.forEach(r=>{let o=0;for(;o<M.length&&M[o].name!==r.track;)o+=1;o<M.length&&(e=0,M[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 ia(){let e=!0,n=4,t=T.trackWidth;M.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&&M.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 mo(t,n.isForward?r:t.sequenceLength-r)}function ca(e){let n=e.path[e.path.length-1],t=y[n.node],r=e.finalNodeCoverLength??0;return mo(t,n.isForward?r:t.sequenceLength-r)}function mo(e,n){let[t,r]=fa(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 ua(){let e,n,t,r,o,s,a,l,i;for(let c=0;c<=k;c+=1)Gn.push(0),$n.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))}),M.sort(Ps),M.forEach(c=>{a="plain",o=we(c,a),s=Is(c);let d=a,h=()=>(d!==a&&(d=a,o=we(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=ao(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(),$.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(),$.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(),$.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,lo(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,lo(-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=ao(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=ca(c),$.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 da=1e9;function so(e,n){an=new Map;let t=new Map,r=T.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],g=h[p+1],b=Math.abs(m),w=Math.abs(g),E=y[b],L=y[w];if(!E||!L)continue;let N;if(r){let x=`${m}>${g}`,A=`${-g}>${-m}`;N=x<A?x:A}else N=`${m}>${g}`;let R=t.get(N);R===void 0?t.set(N,{sSigned:m,dSigned:g,sName:m<0?`-${E.name}`:E.name,dName:g<0?`-${L.name}`:L.name,count:1,sourceTrackID:d.sourceTrackID}):R.count+=1}}let o=0;for(let d of t.values())d.count>o&&(o=d.count);let s=cn,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=da+u,p=`${d.count.toLocaleString()} ${n}${d.count===1?"":"s"}: Node ${d.sName} \u2192 Node ${d.dName}`;an.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 De&&console.log(`[coarsened] ${n}s_in=${e.length} edges_out=${i.length} tallestNodeBeforePlace=${f}(${c.toFixed(1)})`),i}function ao(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,g)=>m.start-g.start),e.features.forEach(m=>{u!==m.type&&(f=we(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=we(o,m.type),$.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&&$.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=we(o,m.type),$.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&&$.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=we(o,u),e.isForward?(h=n+Math.round((m.end+1)*(t-n+1)/p)-1,$.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,$.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 io=7;function lo(e,n,t,r,o,s,a,l,i,u){let f=e===1?$n:Gn,c=10*f[l],d=n+e*(c+5),h=io,p=Math.min(io,o),m=Math.min(t,r),g=Math.max(t,r),b=e===1?n:d-1,w=e===1?d:n,E=x=>{$.verticalRectangles.push({xStart:b,yStart:x,xEnd:w,yEnd:x+o-1,color:s,id:a,name:u,type:i})};E(t),$.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:g-h+1,color:s,id:a,name:u,type:i}),E(r);let L=d+e*h,N=d+e*(h+p),R=`M ${d} ${g}`;R+=` Q ${L} ${g} ${L} ${g-h}`,R+=` H ${N}`,R+=` Q ${N} ${g+o} ${d} ${g+o}`,R+=" Z ",$.corners.push({path:R,color:s,id:a,type:i}),R=`M ${d} ${m}`,R+=` Q ${N} ${m} ${N} ${m+o+h}`,R+=` H ${L}`,R+=` Q ${L} ${m+o} ${d} ${m+o}`,R+=" Z ",$.corners.push({path:R,color:s,id:a,type:i}),f[l]+=1}function fa(e){let n=e.x-4,t=e.x+e.pixelWidth+4;return[n,t]}function co(){switch(y.forEach(e=>{e.sequenceLength===void 0&&(e.sequenceLength=e.seq.length)}),T.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=T.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(`${T.nodeWidthOption} not implemented`)}}function pa(){let e,n,t=Array.from(y,()=>new Set),r=Array.from(y,()=>new Set),o;T.showReads&&S.length>0?o=M.concat(S):o=M,o.forEach(l=>{for(let i=0;i<l.sequence.length;i+=1)X(l.sequence[i])?(e=B(l.sequence[i]),i>0?(n=l.sequence[i-1],X(n)?r[C.get(e)].add(B(n)):(r[C.get(e)].add(n),r[C.get(e)].add(pe(n)))):l.type==="haplotype"&&r[C.get(e)].add("None"),i<l.sequence.length-1?(n=l.sequence[i+1],X(n)?t[C.get(e)].add(B(n)):(t[C.get(e)].add(n),t[C.get(e)].add(pe(n)))):l.type==="haplotype"&&t[C.get(e)].add("None")):(i>0?(e=l.sequence[i-1],t[C.get(l.sequence[i])].add(e),X(e)&&t[C.get(l.sequence[i])].add(B(e))):l.type==="haplotype"&&t[C.get(l.sequence[i])].add("None"),i<l.sequence.length-1?(e=l.sequence[i+1],r[C.get(l.sequence[i])].add(e),X(e)&&r[C.get(l.sequence[i])].add(B(e))):l.type==="haplotype"&&r[C.get(l.sequence[i])].add("None"))});let s=t.map(l=>Array.from(l)),a=r.map(l=>Array.from(l));if(T.showReads&&S.length>0){let l=y.slice();l.sort(un);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=Pe(C.get(d.name),s,a);if(h){let p=f(h)+ne(h).sequenceLength;i.set(d.name,p),i.set(pe(d.name),p),u.set(d.name,c(h)),u.set(pe(d.name),c(h))}else i.set(d.name,0),i.set(pe(d.name),0),u.set(d.name,d.name),u.set(pe(d.name),d.name)}),S.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=B(d.sequence[p]),g=Pe(C.get(m),s,a);if(g&&On(C.get(g),s,a))if(p>0){d.sequence.splice(p,1);let b=ne(g).sequenceLength;h[p].mismatches.forEach(w=>{w.pos+=b}),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 b=h[0],w=f(b.nodeName);b.mismatches.forEach(E=>{E.pos+=w}),b.nodeName=c(b.nodeName)}}})}for(let l=0;l<y.length;l+=1)if(On(l,s,a)&&!Pe(l,s,a)){let i=l;for(;On(i,s,a);)i=C.get(B(a[i][0])),y[l].sequenceLength+=y[i].sequenceLength,y[l].seq+=y[i].seq}M.forEach(l=>{for(let i=l.sequence.length-1;i>=0;i-=1){e=B(l.sequence[i]);let u=C.get(e);Pe(u,s,a)&&l.sequence.splice(i,1)}});for(let l=y.length-1;l>=0;l-=1)Pe(l,s,a)&&y.splice(l,1)}function Pe(e,n,t){if(n[e].length!==1||n[e][0]==="None")return!1;let r=B(n[e][0]),o=C.get(r);return t[o].length!==1||t[o][0]==="None"?!1:r}function On(e,n,t){if(t[e].length!==1||t[e][0]==="None")return!1;let r=B(t[e][0]),o=C.get(r);return!(n[o].length!==1||n[o][0]==="None")}function ha(e){let n=T.focusReadNames?new Set(T.focusReadNames):null;return e.filter(t=>t.is_secondary!==!0&&(t.mapping_quality??0)>=T.mappingQualityCutoff&&(n===null||t.name!==void 0&&n.has(t.name)))}function ma(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 ga(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 go(e,n){let t=ma(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(ga),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 bo(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 yo(e){return e.id.endsWith("-")?"-":"+"}function ya(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=>xe(s.path.name).name===e.referencePath);return o>0&&r.unshift(...r.splice(o,1)),r.flatMap(({path:s,index:a})=>{let l=xe(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;return i.set(c.name,d+1),[(t.get(`${l}	${c.name}`)?.[d]??yo(c))===yo(c)?f:`-${f}`]});return u.length>0?[{id:a,name:s.name,sequence:u,sourceTrackID:0,type:"haplotype"}]:[]})}function xa(e){return e.nodes.map(n=>({name:n.id,seq:"",sequenceLength:Math.max(1,n.length)}))}function Ne(e){return(e.paths?.length??0)>0}var Ma=600,wa=3,La=15;function Na(e){return Math.max(wa,Math.min(La,Math.floor(Ma/Math.max(1,e))))}function xo(e){let n=e.paths??[],t=ya(e);return fo(xa(e),t,[],{nodeWidthOption:"compressed",trackWidth:Na(t.length),trackColor:r=>jn(r.id,n.length)})}function Dn(e){let n=[];return e.nodes.forEach(t=>{t.order>=0&&n.push(t)}),n}function Mo(e,n,t,r){let o=e.y+e.contentHeight/2-n;return[{x:t,y:o},{x:r,y:o}]}function wo(e,n){return{minY:e.bounds.minY-n,maxY:e.bounds.maxY-n}}function Lo(e){let n=Ne(e)?xo(e):void 0;if(!n)return;let t=n.bounds.minY,r={};for(let a of Dn(n))r[a.name]=Mo(a,t,a.x,a.x+a.pixelWidth);let{minY:o,maxY:s}=wo(n,t);return{nodePositions:r,tubeMap:{layout:n,yOffset:t},extent:{minX:n.bounds.minX,maxX:n.bounds.maxX,minY:o,maxY:s}}}function Sa(e,n){let t=new Map(e.nodes.map(l=>[l.id,l])),r=new Map;for(let l of Dn(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&&G(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 Ea(e){return e.nodes.some(G)&&Ne(e)}function No(e){let n=Ea(e)?xo(e):void 0;if(!n)return;let t=Sa(e,n),r=new Map(t.map(a=>[a.order,a])),o=n.bounds.minY,s={};for(let a of Dn(n)){let l=r.get(a.order);s[a.name]=Mo(a,o,l.bp0,l.bp1)}return{nodePositions:s,tubeMap:{layout:n,yOffset:o,columns:t},referenceAxis:!0,pixelRows:!0,extent:wo(n,o)}}function So(e){let n=ce(e);if(n.length!==0)return{nodePositions:ue(n),referenceAxis:!0,pixelRows:!0}}function Eo(e){return e.sample??Ve(e.name)}function Ra(e){return e.haplotype!==void 0&&e.sample!==void 0?`${e.sample}#${e.haplotype}`:Eo(e)}function Aa(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 Ro(e,n,t){let r=e.paths??[],o=r.find(p=>xe(p.name).name===e.referencePath)??r[0];if(!o||r.length<2)return;let s=new Map(e.nodes.map(p=>[p.id,p.length])),a=new Set(o.nodeIds),l=e.anchorPaths?.find(p=>p.name===xe(o.name).name)?.start??0,i=n&&n.end>n.start?n:void 0,u=new Map,f=l;for(let p of o.nodeIds){let m=s.get(p)??0;u.has(p)||u.set(p,{start:f,end:f+m}),f+=m}let c=i===void 0||[...u.values()].some(p=>p.end<=i.start)&&[...u.values()].some(p=>p.start>=i.end),d=p=>{let{ids:m,complete:g}=Aa(p.nodeIds,u,i,c),b=[],w=0,E=0;for(let L of m){let N=s.get(L)??0,R=a.has(L),x=b.at(-1);x?.onReference===R?x.bp+=N:b.push({start:w,bp:N,onReference:R}),w+=N,R||(E+=N)}return{name:p.name,label:Ra(p),sample:Eo(p),bp:w,offReferenceBp:E,complete:g,runs:b}};return{origin:i?i.start:l,unit:t,reference:d(o),rows:r.filter(p=>p!==o).map(d).sort((p,m)=>Number(m.complete)-Number(p.complete)||m.bp-p.bp||p.label.localeCompare(m.label))}}var Ia=1.3;function Ca(e){let n=e.reference.bp;for(let t of e.rows)n=Math.max(n,t.bp);return{maxX:e.origin+n*Ia,maxY:e.rows.length*20}}function Ao(e,n){let t=ce(e),r=Ro(e,n);if(t.length===0||!r)return;let o=ue(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:Ca(r)}}var Ta="Force-directed layout",te=e=>e.nodes.some(G),Io=e=>te(e)&&(e.paths?.length??0)>1,Co=e=>te(e)&&Ne(e),To=e=>te(e)&&e.nodes.some(Ke),qn=[{value:"auto",label:"Anchored",description:"x is reference bp, one row per stable rank. Needs rGFA tags or a reference path.",run:Wt,cutMargins:!0,available:te,drawsLocally:te},{value:"samplerows",label:"Sample rows",description:"x is reference bp, one row per contributing assembly. Needs rGFA tags or a reference path.",run:Vt,cutMargins:!0,available:To,drawsLocally:To},{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:Ao,cutMargins:!1,available:Io,drawsLocally:Io},{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:jt,cutMargins:!1,available:te,drawsLocally:te},{value:"variants",label:"Variant map",description:"The reference as a line with one typed glyph per bubble, from the bubble index or from the graph itself: SNP, indel, deletion, inversion, repeat array. Click a glyph to open the graph inside it.",run:So,cutMargins:!0,available:te,drawsLocally:te},{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:Lo,cutMargins:!1,available:Ne,drawsLocally:Ne},{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:No,cutMargins:!0,available:Co,drawsLocally:Co},{value:"force",label:Ta,description:"OGDF FMMM, via the external Bandage engine.",run:()=>{},cutMargins:!1,available:()=>!0,drawsLocally:()=>!1}],Xl=qn.map(e=>e.value);function zl(e,n){return!va(e).drawsLocally(n)}function va(e){return qn.find(n=>n.value===e)??qn[0]}function Oa(e,n){return Math.abs(n-e)<=.1*Math.max(e,n)}function Ba(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 vo(e,n){return e.map((t,r)=>({item:t,i:r})).sort((t,r)=>n(t.item)-n(r.item))}var Oo=e=>e.spanningReads!==0;function _a(e,n){let t=new Map,r=new Set,o=vo(e,a=>a.bp),s=[[o.filter(a=>a.item.complete),n.filter(Oo)],[o,n]];for(let[a,l]of s){let i=a.filter(f=>!t.has(f.i)),u=vo(l.filter(f=>!r.has(f)),f=>f.bp);Ba(i.map(f=>f.item.bp),u.map(f=>f.item.bp)).forEach((f,c)=>{if(f!==void 0){let d=i[c],h=u[f].item;r.add(h),t.set(d.i,{...h,agrees:d.item.complete&&Oo(h)?Oa(d.item.bp,h.bp):void 0})}})}return t}function Vl(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 _a(s,a))t.set(s[l],i)}return e.map(o=>({...o,call:t.get(o)}))}function Po(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 Do(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 Bo(e,n,t){let r=Po(e);return{commands:r,color:n,alpha:t,...Do(r)}}function _o(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=go(t.curves,n).map(a=>Bo(a.path,a.color,a.alpha??1)),s=t.corners.filter(a=>a.type===n).map(a=>Bo(a.path,a.color,1));return{rects:r,shapes:[...o,...s]}}function kl(e){let n=[];return e.nodes.forEach(t=>{if(t.order>=0){let r=Po(bo(t));n.push({name:t.name,commands:r,...Do(r)})}}),{layers:[_o(e,"haplotype"),_o(e,"read")],nodes:n}}function Go(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 $o(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}function Ql(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)$o(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)}),$o(e,u.shapes,a,f=>{Go(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(),Go(e,u.commands,r,o),e.fillStyle=f?"rgba(255,192,203,0.5)":i,e.strokeStyle=f?"#ff0000":l,e.fill(),e.stroke()}}var qo=me(je(),1),ec=(0,ee.createSvgIcon)((0,qo.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{ee as a,is as b,Nn as c,ls as d,ge as e,ji as f,Vi as g,Ki as h,Ui as i,qa as j,rt as k,at as l,ke as m,it as n,Ls as o,Ya as p,ja as q,Va as r,Ka as s,gt as t,xr as u,xt as v,ri as w,Ro as x,Ca as y,Ge as z,Rr as A,qn as B,Xl as C,zl as D,va as E,ii as F,Rt as G,Yr as H,Mi as I,wi as J,Vl as K,kl as L,Ql as M,Ei as N,Ri as O,os as P,Ai as Q,Ii as R,En as S,Fi as T,Hi as U,ec as V};
+//# sourceMappingURL=chunk-VXUVWSY5.js.map