jbrowse-plugin-msaview 2.5.2 → 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, treeFileLocation, connectedViewId, connectedFeature, displayName, colorSchemeName, colWidth, rowHeight, treeAreaWidth, treeWidth, drawNodeBubbles, labelsAlignRight, showBranchLen, querySeqName, highlightColumns, }) => {
-        if (!data && !msaFileLocation) {
+    ({ 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,6 +23,8 @@ export default function LaunchMsaViewExtensionPointF(pluginManager) {
                 msaData: data?.msa,
                 treeData: data?.tree,
                 msaUrl: msaFileLocation?.uri,
+                msaIndexedLocation,
+                msaName,
                 treeUrl: treeFileLocation?.uri,
                 querySeqName,
             },

package/dist/MsaViewPanel/afterCreateAutoruns.js

@@ -1,5 +1,6 @@
 import { getSession } from '@jbrowse/core/util';
 import { doLaunchBlast } from './doLaunchBlast';
+import { fetchIndexedMsa } from './fetchIndexedMsa';
 import { genomeToMSA } from './genomeToMSA';
 import { loadProteinDomains } from './loadProteinDomains';
 import { cleanupOldData, generateDataStoreId, retrieveMsaData, storeMsaData, } from './msaDataStore';
@@ -121,7 +122,7 @@ export function processInit(self) {
         void (async () => {
             try {
                 self.setError(undefined);
-                const { msaData, msaUrl, treeData, treeUrl, querySeqName } = init;
+                const { msaData, msaUrl, msaIndexedLocation, msaName, treeData, treeUrl, querySeqName, } = init;
                 if (msaUrl) {
                     const id = getUniprotIdFromAlphaFoldUrl(msaUrl);
                     if (id) {
@@ -143,6 +144,18 @@ export function processInit(self) {
                     const data = await response.text();
                     self.setMSA(data);
                 }
+                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) {
                     self.setTree(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,6 +1,10 @@
 export interface MsaViewInitState {
     msaData?: string;
     msaUrl?: string;
+    msaIndexedLocation?: {
+        uri: string;
+    };
+    msaName?: string;
     treeData?: string;
     treeUrl?: string;
     querySeqName?: string;