e2abb3d38e8bf2296e4be847deff63409bc82d1c angie Thu Dec 7 12:55:15 2017 -0800 In vpPep, don't distinguish between spansUtrCds and cantPredict. (btw spansUtrCds was only used for 5' UTR/CDS boundary) This should not affect output of anything because they were treated the same anyway. diff --git src/hg/lib/hgHgvs.c src/hg/lib/hgHgvs.c index e3f24a6..0557301 100644 --- src/hg/lib/hgHgvs.c +++ src/hg/lib/hgHgvs.c @@ -2791,31 +2791,31 @@ if (addParens) dyStringAppendC(dy, '('); int refLen = vpPep->end - vpPep->start; // When predicting frameshift/extension, the length of ref may be different from refLen int refExtLen = vpPep->ref ? strlen(vpPep->ref) : refLen; int altLen = vpPep->alt ? strlen(vpPep->alt) : 0; char refStartAbbr[4]; if (vpPep->ref) aaToAbbr(vpPep->ref[0], refStartAbbr, sizeof(refStartAbbr)); else // If ref is null then we should be writing just '=' or '?' but prevent garbage just in case: safecpy(refStartAbbr, sizeof(refStartAbbr), "?"); // protSeq may or may not end with X, so treat protSeq->size accordingly boolean hitsStopCodon = (vpPep->end > protSeq->size || ((protSeq->dna[protSeq->size-1] == 'X') && vpPep->end == protSeq->size)); -if (vpPep->cantPredict || vpPep->spansUtrCds || isStartLoss(vpPep)) +if (vpPep->cantPredict || isStartLoss(vpPep)) dyStringAppend(dy, "?"); else if (vpPep->frameshift) { dyStringPrintf(dy, "%s%d", refStartAbbr, vpPep->start+1); if (altLen == 1) dyStringAppend(dy, "Ter"); else { char altStartAbbr[4]; aaToAbbr(vpPep->alt[0], altStartAbbr, sizeof(altStartAbbr)); // For stop-loss extension, make it "ext*" if (hitsStopCodon && altLen > refExtLen) dyStringPrintf(dy, "%sext*%d", altStartAbbr, altLen - refExtLen); else dyStringPrintf(dy, "%sfsTer%d", altStartAbbr, altLen);