03ffd59683ac24c0eb8ac90b8f5589f71570a87f angie Tue Feb 27 14:16:28 2018 -0800 HGVS p. term: if a single amino acid is duplicated, show the single position instead of range. diff --git src/hg/lib/hgHgvs.c src/hg/lib/hgHgvs.c index 7f70293..c295038 100644 --- src/hg/lib/hgHgvs.c +++ src/hg/lib/hgHgvs.c @@ -2844,23 +2844,26 @@ } else if (refLen == 0) { // Insertion; expand to two-AA range around insertion point rangeStart--; aaToAbbr(pSeq[rangeStart], refStartAbbr, sizeof(refStartAbbr)); rangeEnd++; } hitsStopCodon = (rangeEnd > protSeq->size || ((protSeq->dna[protSeq->size-1] == 'X') && rangeEnd == protSeq->size)); char refLastAbbr[4]; if (hitsStopCodon) aaToAbbr('X', refLastAbbr, sizeof(refLastAbbr)); else aaToAbbr(pSeq[rangeEnd-1], refLastAbbr, sizeof(refLastAbbr)); + if (dupLen == 1) + dyStringPrintf(dy, "%s%d", refStartAbbr, rangeStart+1); + else dyStringPrintf(dy, "%s%d_%s%d", refStartAbbr, rangeStart+1, refLastAbbr, rangeEnd); } hgvsAppendChangesFromPepRefAlt(dy, vpPep->ref, vpPep->alt, dupLen); } if (addParens) dyStringAppendC(dy, ')'); return dyStringCannibalize(&dy); }