2b50ebde25b8ad6c5a43797ae0ceb85ae66fd0bc markd Sun Jun 11 13:58:19 2023 -0700 fixed a couple of calls to pslTransMap to add new alignments diff --git src/hg/lib/hgHgvs.c src/hg/lib/hgHgvs.c index 1a27a47..23a2428 100644 --- src/hg/lib/hgHgvs.c +++ src/hg/lib/hgHgvs.c @@ -1756,31 +1756,31 @@ } } // Not contained in a deletion from reference genome (txAli target) -- return NULL. return NULL; } static struct psl *mapPsl(struct hgvsVariant *hgvs, char *acc, struct psl *txAli, struct genbankCds *cds, int *retUpstream, int *retDownstream) /* Use pslTransMap to map hgvs onto the genome. */ { // variantPsl contains the anchor if a non-cdsStart anchor is used because // the actual position might be outside the bounds of the transcript sequence (intron/UTR) struct psl *variantPsl = pslFromHgvsNuc(hgvs, acc, txAli->qSize, txAli->qEnd, cds, retUpstream, retDownstream); -struct psl *mappedToGenome = pslTransMap(pslTransMapNoOpts, variantPsl, txAli); +struct psl *mappedToGenome = pslTransMap(pslTransMapNoOpts, variantPsl, pslTypeUnspecified, txAli, pslTypeUnspecified); // If there is a deletion in the genome / insertion in the transcript then pslTransMap cannot // map those bases to the genome. In that case take a harder look and return the deletion // coords. if (mappedToGenome == NULL) mappedToGenome = mapToDeletion(variantPsl, txAli); pslFree(&variantPsl); return mappedToGenome; } static char *pslTableForAcc(char *db, char *acc) /* Based on acc (and whether db has NCBI RefSeq alignments), pick a PSL (or genePred for GENCODE) * table where acc should be found (table may or may not exist). Don't free the returned string. */ { char *pslTable = NULL; if (startsWith("LRG_", acc))