a75cfac22803c312074299aa0dcf286ca336ce89 braney Wed Aug 27 13:28:08 2025 -0700 ongoing work on quickLift hgc page plus a fix for genepreds that don't successfully quickLift diff --git src/hg/hgc/hgc.c src/hg/hgc/hgc.c index 25374b41a83..bd49bb6affe 100644 --- src/hg/hgc/hgc.c +++ src/hg/hgc/hgc.c @@ -3916,63 +3916,30 @@ atoi(words[0])); } } if (subChain->qStrand == '-') reverseIntRange(&qStart, &qEnd, subChain->qSize); qSeq = otherChromSeq(otherTbf, otherDb, subChain->qName, qStart, qEnd); if (subChain->qStrand == '-') reverseComplement(qSeq->dna, qSeq->size); subChain->score = chainCalcScoreSubChain(subChain, scoreScheme, gapCalc, qSeq, tSeq); subSetScore = subChain->score; } chainFree(&toFree); -char *liftDb = cloneString(trackDbSetting(tdb, "quickLiftDb")); - -if (liftDb != NULL) - { - int seqStart = cartInt(cart, "l"); - int seqEnd = cartInt(cart, "r"); - char *chromName = cartString(cart, "c"); - char *quickLiftFile = cloneString(trackDbSetting(tdb, "quickLiftUrl")); - struct quickLiftRegions *hr, *regions = quickLiftGetRegions(database, liftDb, quickLiftFile, chromName, seqStart, seqEnd); - - printf("Variants in Window:
"); - printf("\n"); - printf(""); - for(hr = regions; hr; hr = hr->next) - { - if (hr->type == QUICKTYPE_NOTHING) - continue; - - char position[128]; - char *ourPos, *otherPos; - snprintf(position, 128, "%s:%ld-%ld", hr->chrom, hr->chromStart, hr->chromEnd); - ourPos = cloneString(addCommasToPos(database, position)); - snprintf(position, 128, "%s:%ld-%ld", hr->oChrom, hr->oChromStart, hr->oChromEnd); - otherPos = cloneString(addCommasToPos(database, position)); - printf("
TypePosOther PosOther CountOther BasesCountBases
%s%s%s%d%.*s%d%.*s", quickTypeStrings[hr->type], ourPos,otherPos, hr->otherBaseCount, hr->otherBaseCount, hr->otherBases, hr->baseCount, hr->baseCount, hr->bases); - - printf(""); - hgcAnchorSomewhereExt("htcChainAli", item, tdb->track, chain->tName, hr->chromStart - 10, hr->chromEnd + 10, tdb->track); - printf("DNA alignment"); - } - printf("
"); - } - printf("%s position: %s:%d-%d" " size: %d
\n", trackHubSkipHubName(thisOrg), hgTracksName(), cartSidUrlString(cart), database, chain->tName, chain->tStart+1, chain->tEnd, chain->tName, chain->tStart+1, chain->tEnd, chain->tEnd-chain->tStart); printf("Strand: %c
\n", chain->qStrand); qChainRangePlusStrand(chain, &qs, &qe); if (sameWord(otherDb, "seq")) { printf("%s position: %s:%d-%d size: %d
\n", otherOrg, chain->qName, qs, qe, chain->qEnd - chain->qStart); } else { /* prints link to other db browser only if db exists and is active */