54406b80d5d435970989acff7b22dd6146c6b411 braney Sat Jan 22 15:24:59 2022 -0800 adding chrom alias support to big files diff --git src/hg/hgc/clinvarSubLolly.c src/hg/hgc/clinvarSubLolly.c index a0b78b2..1294a96 100644 --- src/hg/hgc/clinvarSubLolly.c +++ src/hg/hgc/clinvarSubLolly.c @@ -1,47 +1,48 @@ #include "common.h" #include "wiggle.h" #include "cart.h" #include "hgc.h" #include "hCommon.h" #include "hgColors.h" #include "bigBed.h" #include "hui.h" #include "subText.h" #include "trackHub.h" #include "clinvarSubLolly.h" +#include "chromAlias.h" char *statusByScore[] = { "Other", "Benign", "Likely Benign", "Uncertain significance", "Likely pathogenic", "Pathogenic", }; void printSubmissions(struct trackDb *tdb, char *chrom, int start, int end, unsigned wantScore, int numSubs, boolean not) /* Print out all the submissions at this position that match the * desired score (ie. clinical status). * If "not" is true, then print out all submissions that do NOT match score. * We're grabbing these records out of a second bigBed that's referenced by * "xrefDataUrl" in the trackDb. */ { char *xrefDataUrl = hReplaceGbdb(trackDbSetting(tdb, "xrefDataUrl")); -struct bbiFile *bbi = bigBedFileOpen(xrefDataUrl); +struct bbiFile *bbi = bigBedFileOpenAlias(xrefDataUrl, chromAliasGetHash(database)); struct asObject *as = bigBedAsOrDefault(bbi); struct lm *lm = lmInit(0); struct bigBedInterval *bbList = bigBedIntervalQuery(bbi, chrom, start, end, 0, lm); int count; if (not) count = slCount(bbList) - numSubs; else count = numSubs; // nothing to print if (count == 0) return; if (!not) @@ -68,31 +69,31 @@ } void doClinvarSubLolly(struct trackDb *tdb, char *item) /* Put up page for clinvarSubLolly track. */ { genericHeader(tdb, NULL); struct sqlConnection *conn = NULL; if (!trackHubDatabase(database)) conn = hAllocConnTrack(database, tdb); int start = cartInt(cart, "o"); int end = cartInt(cart, "t"); char *chrom = cartString(cart, "c"); char *fileName = bbiNameFromSettingOrTable(tdb, conn, tdb->table); -struct bbiFile *bbi = bigBedFileOpen(hReplaceGbdb(fileName)); +struct bbiFile *bbi = bigBedFileOpenAlias(hReplaceGbdb(fileName), chromAliasGetHash(database)); struct lm *lm = lmInit(0); struct bigBedInterval *bbList = bigBedIntervalQuery(bbi, chrom, start, end, 0, lm); struct bigBedInterval *bb; char *fields[bbi->fieldCount]; for (bb = bbList; bb != NULL; bb = bb->next) { if (!(bb->start == start && bb->end == end)) continue; // our names are unique char *name = cloneFirstWordByDelimiterNoSkip(bb->rest, '\t'); boolean match = (isEmpty(name) && isEmpty(item)) || sameOk(name, item); if (!match) continue;