bd8427c06c3c63aa073850539abcb6ea284d43b5 braney Fri Oct 7 10:52:30 2016 -0700 fix a problem in mafGene with translated 1bp exons. diff --git src/hg/lib/mafGene.c src/hg/lib/mafGene.c index 7c2383e..7a545b9 100644 --- src/hg/lib/mafGene.c +++ src/hg/lib/mafGene.c @@ -267,31 +267,34 @@ aaSeq *outSeq; int exonCount = 0; struct exonInfo *gi = giList; for(; gi; gi = gi->next) { if (gi->exonSize > 1) exonCount++; } for(gi = giList; gi; gi = gi->next, exonNum++) { struct speciesInfo *siTemp = si; if (gi->exonSize == 1) + { + siTemp->curPosString = siTemp->curPosString->next; continue; + } for(; siTemp ; siTemp = siTemp->next) { char *ptr = exonBuffer; switch(gi->frame) { case 0: /* just copy the sequence over */ memcpy(ptr, &siTemp->nucSequence[gi->exonStart], gi->exonSize); ptr += gi->exonSize; break; case 1: /* we need to grab one nuc from the end * of the previous exon */ *ptr++ = siTemp->nucSequence[gi->exonStart - 1];