--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/x/static/js/6391.aabba4a1.chunk.js	Sat Oct 05 23:58:05 2024 +0000
@@ -0,0 +1,2 @@
+(globalThis.webpackChunk_jbrowse_web=globalThis.webpackChunk_jbrowse_web||[]).push([[6391],{78369:(e,t,r)=>{"use strict";r.r(t),r.d(t,{default:()=>u});var a=r(75099),s=r(46377),i=r(99546),n=r(99834),o=r(66885),c=r(82088),d=r(44728);class h{constructor(e,t){this.record=e,this._store=t}_get_name(){return this.record.readName}_get_start(){return this.record.alignmentStart-1}_get_end(){return this.record.alignmentStart+(this.record.lengthOnRef??1)-1}_get_cram_read_features(){return this.record.readFeatures}_get_type(){return"match"}_get_score(){return this.record.mappingQuality}_get_flags(){return this.record.flags}_get_strand(){return this.record.isReverseComplemented()?-1:1}_read_group_id(){return this._store.samHeader.readGroups?.[this.record.readGroupId]}_get_qual(){return(this.record.qualityScores||[]).join(" ")}qualRaw(){return this.record.qualityScores}_get_refName(){return this._store.refIdToName(this.record.sequenceId)}_get_is_paired(){return!!this.record.mate}_get_pair_orientation(){return this.record.isPaired()?this.record.getPairOrientation():void 0}_get_template_length(){return this.record.templateLength||this.record.templateSize}_get_next_ref(){return this.record.mate?this._store.refIdToName(this.record.mate.sequenceId):void 0}_get_next_segment_position(){return this.record.mate?`${this._store.refIdToName(this.record.mate.sequenceId)}:${this.record.mate.alignmentStart}`:void 0}_get_next_pos(){return this.record.mate?.alignmentStart}_get_tags(){const e=this._read_group_id(),{tags:t}=this.record;return void 0!==e?{...t,RG:e}:t}_get_seq(){return this.record.getReadBases()}_get_CIGAR(){return function(e,t,r,a){let s="",i="",n="M",o=0;if(!a)return"";const c=a.seq,d=a.start;let h=t,u=0,g=0;if(void 0!==e)for(const{code:t,refPos:r,sub:a,data:l}of e)if(u=r-h,s+=c.slice(h-d,r-d),h=r,g>0&&u&&(i+=`${g}I`,g=0),o&&"M"!==n&&(i+=`${o}${n}`,o=0),u&&(n="M",o+=u),"b"===t){const e=l.split(","),t=String.fromCharCode(...e);s+=t,h+=t.length,o+=t.length}else"B"===t||"X"===t?(s+=a,h++,o++):"D"===t||"N"===t?(h+=l,o&&(i+=`${o}${n}`),i+=l+t,o=0):"I"===t||"S"===t?(s+=l,o&&(i+=`${o}${n}`),i+=l.length+t,o=0):"i"===t?(o&&(i+=`${o}${n}`),g++,s+=l,o=0):"P"===t?(o&&(i+=`${o}${n}`),i+=`${l}P`):"H"===t&&(o&&(i+=`${o}${n}`),i+=`${l}H`,o=0);else u=r-s.length;return s.length!==r&&(u=r-s.length,s+=c.slice(h-d,h-d+u),o&&"M"!==n&&(i+=`${o}${n}`,o=0),n="M",o+=u),u&&g>0&&(i+=`${g}I`),o&&(i+=`${o}${n}`),i}(this.record.readFeatures,this.record.alignmentStart,this.record.readLength,this.record._refRegion)}tags(){return Object.getOwnPropertyNames(h.prototype).filter((e=>e.startsWith("_get_")&&"_get_mismatches"!==e&&"_get_cram_read_features"!==e)).map((e=>e.replace("_get_","")))}id(){return`${this._store.id}-${this.record.uniqueId}`}get(e){const t=`_get_${e}`;if(this[t])return this[t]()}parent(){}children(){}set(){}pairedFeature(){return!1}_get_clipPos(){const e=this.get("mismatches");if(e.length){const t=-1===this.get("strand")?e.at(-1):e[0],{type:r,cliplen:a}=t;if("softclip"===r||"hardclip"===r)return a}return 0}toJSON(){return{...Object.fromEntries(this.tags().map((e=>[e,this.get(e)])).filter((e=>void 0!==e[1]))),uniqueId:this.id()}}_get_mismatches(){const e=this.record.readFeatures,t=this.qualRaw();return function(e,t,r){if(!e)return[];const a=new Array(e.length);let s=0,i=0,n=0,o=0,c=t;for(const{refPos:d,code:h,pos:u,data:g,sub:l,ref:m}of e)if(o=n-c,c=n,o&&i>0&&(a[s++]={start:n,type:"insertion",base:`${i}`,length:0},i=0),n=d-1-t,"X"===h)a[s++]={start:n,length:1,base:l,qual:r?.[u-1],altbase:m?.toUpperCase(),type:"mismatch"};else if("I"===h)a[s++]={start:n,type:"insertion",base:`${g.length}`,length:0};else if("N"===h)a[s++]={type:"skip",length:g,start:n,base:"N"};else if("S"===h){const e=g.length;a[s++]={start:n,type:"softclip",base:`S${e}`,cliplen:e,length:1}}else if("P"===h);else if("H"===h){const e=g;a[s++]={start:n,type:"hardclip",base:`H${e}`,cliplen:e,length:1}}else"D"===h?a[s++]={type:"deletion",length:g,start:n,base:"*"}:"b"===h||"q"===h||"B"===h||"i"===h&&i++;return o&&i>0&&(a[s++]={start:n,type:"insertion",base:`${i}`,length:0},i=0),a.slice(0,s)}(e,this.get("start"),t)}}class u extends s.BaseFeatureDataAdapter{samHeader={};seqIdToOriginalRefName=[];async configurePre(){const e=this.getConf("cramLocation"),t=this.getConf("craiLocation"),r=this.pluginManager,s=new a.bQ({cramFilehandle:(0,n.openLocation)(e,r),index:new a.Wb({filehandle:(0,n.openLocation)(t,r)}),seqFetch:(...e)=>this.seqFetch(...e),checkSequenceMD5:!1});if(!this.getSubAdapter)throw new Error("Error getting subadapter");const i=this.getConf("sequenceAdapter");if(!i)throw new Error("no sequenceAdapter supplied to CramAdapter config");return{cram:s,sequenceAdapter:(await this.getSubAdapter(i)).dataAdapter}}async configure(){return this.configureP||(this.configureP=this.configurePre().catch((e=>{throw this.configureP=void 0,e}))),this.configureP}async getHeader(e){const{cram:t}=await this.configure();return t.cram.getHeaderText()}async seqFetch(e,t,r){t-=1;const{sequenceAdapter:a}=await this.configure(),s=this.refIdToOriginalName(e)||this.refIdToName(e);if(!s)throw new Error("unknown");const n=(await(0,d._)(a.getFeatures({refName:s,start:t,end:r,assemblyName:""}).pipe((0,c.$)()))).sort(((e,t)=>e.get("start")-t.get("start"))).map((e=>{const a=e.get("start"),s=e.get("end"),i=Math.max(t-a,0),n=Math.min(r-a,s-a)-i;return(e.get("seq")||e.get("residues")).slice(i,i+n)})).join(""),o=r-t;if(n.length!==o)throw new Error(`fetching ${s}:${(0,i.toLocale)(t-1)}-${(0,i.toLocale)(r)} returned ${(0,i.toLocale)(n.length)} bases, should have returned ${(0,i.toLocale)(o)}`);return n}async setupPre(e){const{statusCallback:t=()=>{}}=e||{};return(0,i.updateStatus)("Downloading index",t,(async()=>{const e=await this.configure(),{cram:t}=e,r=await t.cram.getSamHeader(),a=[],s={};r.filter((e=>"SQ"===e.tag)).forEach(((e,t)=>{const r=e.data.find((e=>"SN"===e.tag));if(r){const e=r.value;s[e]=t,a[t]=e}}));const i=r.filter((e=>"RG"===e.tag)).map((e=>e.data.find((e=>"ID"===e.tag))?.value)),n={idToName:a,nameToId:s,readGroups:i};return this.samHeader=n,{samHeader:n,...e}}))}async setup(e){return this.setupP||(this.setupP=this.setupPre(e).catch((e=>{throw this.setupP=void 0,e}))),this.setupP}async getRefNames(e){const{samHeader:t}=await this.setup(e);if(!t.idToName)throw new Error("CRAM file has no header lines");return t.idToName}refNameToId(e){return this.samHeader.nameToId?this.samHeader.nameToId[e]:this.seqIdToRefName?this.seqIdToRefName.indexOf(e):void 0}refIdToName(e){return this.samHeader.idToName?.[e]||this.seqIdToRefName?.[e]}refIdToOriginalName(e){return this.seqIdToOriginalRefName[e]}getFeatures(e,t){const{signal:r,filterBy:a,statusCallback:s=()=>{}}=t||{},{refName:n,start:c,end:d,originalRefName:h}=e;return(0,o.ObservableCreate)((async e=>{const{cram:o,samHeader:u}=await this.setup(t),g=this.refNameToId(n);if(void 0===g)return console.warn("Unknown refName",n),void e.complete();h&&(this.seqIdToOriginalRefName[g]=h);const l=await(0,i.updateStatus)("Downloading alignments",s,(()=>o.getRecordsForRange(g,c,d)));(0,i.checkAbortSignal)(r),await(0,i.updateStatus)("Processing alignments",s,(()=>{const{flagInclude:t=0,flagExclude:r=0,tagFilter:s,readName:i}=a||{};for(const a of l){const n=a.flags;if((n&t)===t||n&r){if(s){const e="RG"===s.tag?u.readGroups?.[a.readGroupId]:a.tags[s.tag],t=s.value;if("*"===t?void 0!==e:`${e}`!=`${t}`)continue}i&&a.readName!==i||e.next(this.cramRecordToFeature(a))}}e.complete()}))}),r)}freeResources(){}cramRecordToFeature(e){return new h(e,this)}async getMultiRegionFeatureDensityStats(e,t){return{bytes:await this.bytesForRegions(e,t),fetchSizeLimit:this.getConf("fetchSizeLimit")}}async bytesForRegions(e,t){const{cram:r}=await this.configure();return(await Promise.all(e.map((e=>{const{refName:t,start:a,end:s}=e,i=this.refNameToId(t);return void 0!==i?r.index.getEntriesForRange(i,a,s):[{sliceBytes:0}]})))).flat().reduce(((e,t)=>e+t.sliceBytes),0)}}},81239:()=>{}}]);
+//# sourceMappingURL=6391.aabba4a1.chunk.js.map
\ No newline at end of file