jbrowse-plugin-msaview 2.6.0 → 2.6.1

package/dist/LaunchMsaViewExtensionPoint/index.js

@@ -1,8 +1,8 @@
 export default function LaunchMsaViewExtensionPointF(pluginManager) {
     pluginManager.addToExtensionPoint('LaunchView-MsaView', 
     // @ts-expect-error
-    ({ session, data, msaFileLocation, msaTabixLocation, msaIndexLocation, msaId, treeFileLocation, connectedViewId, connectedFeature, displayName, colorSchemeName, colWidth, rowHeight, treeAreaWidth, treeWidth, drawNodeBubbles, labelsAlignRight, showBranchLen, querySeqName, highlightColumns, }) => {
-        if (!data && !msaFileLocation && !msaTabixLocation) {
+    ({ session, data, msaFileLocation, msaIndexedLocation, msaName, treeFileLocation, connectedViewId, connectedFeature, displayName, colorSchemeName, colWidth, rowHeight, treeAreaWidth, treeWidth, drawNodeBubbles, labelsAlignRight, showBranchLen, querySeqName, highlightColumns, }) => {
+        if (!data && !msaFileLocation && !msaIndexedLocation) {
             throw new Error('No MSA data or file location provided when launching MSA view');
         }
         session.addView('MsaView', {
@@ -23,9 +23,8 @@ export default function LaunchMsaViewExtensionPointF(pluginManager) {
                 msaData: data?.msa,
                 treeData: data?.tree,
                 msaUrl: msaFileLocation?.uri,
-                msaTabixLocation,
-                msaIndexLocation,
-                msaId,
+                msaIndexedLocation,
+                msaName,
                 treeUrl: treeFileLocation?.uri,
                 querySeqName,
             },

package/dist/MsaViewPanel/afterCreateAutoruns.js

@@ -1,6 +1,6 @@
 import { getSession } from '@jbrowse/core/util';
 import { doLaunchBlast } from './doLaunchBlast';
-import { fetchTabixMsa } from './fetchTabixMsa';
+import { fetchIndexedMsa } from './fetchIndexedMsa';
 import { genomeToMSA } from './genomeToMSA';
 import { loadProteinDomains } from './loadProteinDomains';
 import { cleanupOldData, generateDataStoreId, retrieveMsaData, storeMsaData, } from './msaDataStore';
@@ -122,7 +122,7 @@ export function processInit(self) {
         void (async () => {
             try {
                 self.setError(undefined);
-                const { msaData, msaUrl, msaTabixLocation, msaIndexLocation, msaId, treeData, treeUrl, querySeqName, } = init;
+                const { msaData, msaUrl, msaIndexedLocation, msaName, treeData, treeUrl, querySeqName, } = init;
                 if (msaUrl) {
                     const id = getUniprotIdFromAlphaFoldUrl(msaUrl);
                     if (id) {
@@ -144,23 +144,16 @@ export function processInit(self) {
                     const data = await response.text();
                     self.setMSA(data);
                 }
-                else if (msaTabixLocation) {
-                    const feature = self.connectedFeature;
-                    if (feature) {
-                        const fasta = await fetchTabixMsa({
-                            location: msaTabixLocation,
-                            indexLocation: msaIndexLocation,
-                            msaId: msaId ?? String(feature.name),
-                            refName: String(feature.refName),
-                            start: Number(feature.start),
-                            end: Number(feature.end),
-                        });
-                        if (fasta) {
-                            self.setMSA(fasta);
-                        }
-                        else {
-                            throw new Error(`No alignment for ${msaId ?? String(feature.name)} in ${msaTabixLocation.uri}`);
-                        }
+                else if (msaIndexedLocation && msaName) {
+                    const fasta = await fetchIndexedMsa({
+                        location: msaIndexedLocation,
+                        name: msaName,
+                    });
+                    if (fasta) {
+                        self.setMSA(fasta);
+                    }
+                    else {
+                        throw new Error(`No alignment named ${msaName} in ${msaIndexedLocation.uri}`);
                     }
                 }
                 if (treeData) {

package/dist/MsaViewPanel/fetchIndexedMsa.d.ts

@@ -0,0 +1,6 @@
+export declare function fetchIndexedMsa({ location, name, }: {
+    location: {
+        uri: string;
+    };
+    name: string;
+}): Promise<string | undefined>;

package/dist/MsaViewPanel/fetchIndexedMsa.js

@@ -0,0 +1,36 @@
+import { BgzfFilehandle } from '@gmod/bgzf-filehandle';
+import { openLocation } from '@jbrowse/core/util/io';
+// Pull one transcript's whole multiple-alignment out of a single bgzip file
+// keyed by NAME. Given the `.fa.gz` uri, the bgzip index (`.gzi`) and name index
+// (`.idx`, a TSV `name<TAB>offset<TAB>length` of each block's uncompressed byte
+// offset + length) are found by suffix. We fetch the small `.idx` once, look up
+// the name, and random-read just that block — already valid FASTA
+// (`>hg38\nSEQ\n>panTro4\nSEQ\n...`). One genome-scale alignment serves any gene
+// with no per-gene files or coordinates. See react-msaview's gene-explorer.
+export async function fetchIndexedMsa({ location, name, }) {
+    const open = (uri) => openLocation({ uri, locationType: 'UriLocation' });
+    const idxText = await open(`${location.uri}.idx`).readFile('utf8');
+    const entry = lookup(idxText, name);
+    if (entry) {
+        const fh = new BgzfFilehandle({
+            filehandle: open(location.uri),
+            gziFilehandle: open(`${location.uri}.gzi`),
+        });
+        const bytes = await fh.read(entry.length, entry.offset);
+        return new TextDecoder().decode(bytes).trim();
+    }
+    return undefined;
+}
+// ids are matched versionless, so version drift between the alignment build and
+// the live annotation never breaks the lookup
+const versionless = (s) => s.replace(/\.\d+$/, '');
+function lookup(idxText, name) {
+    const want = versionless(name);
+    for (const line of idxText.split('\n')) {
+        const [id, offset, length] = line.split('\t');
+        if (id && versionless(id) === want) {
+            return { offset: Number(offset), length: Number(length) };
+        }
+    }
+    return undefined;
+}

package/dist/MsaViewPanel/types.d.ts

@@ -1,13 +1,10 @@
 export interface MsaViewInitState {
     msaData?: string;
     msaUrl?: string;
-    msaTabixLocation?: {
+    msaIndexedLocation?: {
         uri: string;
     };
-    msaIndexLocation?: {
-        uri: string;
-    };
-    msaId?: string;
+    msaName?: string;
     treeData?: string;
     treeUrl?: string;
     querySeqName?: string;