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("
| Type | Pos | Other Pos | Other Count | Other Bases | Count | Bases | "); - 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("|
| %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(" |