d8a949135724f97cfc92b52111ac223fece3dbf5 markd Mon Jun 29 14:53:25 2020 -0700 got entire tree compliing with geneFind changes diff --git src/hg/lib/coordConv.c src/hg/lib/coordConv.c index cd404f9..8407da3 100644 --- src/hg/lib/coordConv.c +++ src/hg/lib/coordConv.c @@ -255,34 +255,34 @@ if(seq == NULL || db == NULL) errAbort("coordConv::doDnaAlignment() - dnaSeq and/or db can't be NULL."); if(strlen(seq->dna) != seq->size) errAbort("coordConv::doDnaAlignment() - there seems to be something fishy about %s: the size doesn't equal the length", seq->name); /* if there are too many n's it can cause the blat server to hang */ if(strstr(seq->dna, "nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn") ) return NULL; makeTempName(&pslTn,"ccR", ".psl"); f = mustOpen(pslTn.forCgi, "w"); gvo = gfOutputPsl(920, FALSE, FALSE, f, FALSE, FALSE); gfOutputHead(gvo, f); /* align to genome, both strands */ conn = gfConnect(blatHost, port); -gfAlignStrand(&conn, nibDir, seq, FALSE, 20, tFileCache, gvo); +gfAlignStrand(&conn, nibDir, seq, FALSE, 20, tFileCache, gvo, NULL); reverseComplement(seq->dna, seq->size); conn = gfConnect(blatHost, port); -gfAlignStrand(&conn, nibDir, seq, TRUE, 20 , tFileCache, gvo); +gfAlignStrand(&conn, nibDir, seq, TRUE, 20 , tFileCache, gvo, NULL); gfOutputQuery(gvo, f); carefulClose(&f); pslList = pslLoadAll(pslTn.forCgi); remove(pslTn.forCgi); gfOutputFree(&gvo); return pslList; } void getAlignmentsForSeqs(struct coordConvRep *ccr, char *blatHost, char *port, char *nibDir) /* Do alignments for the the dnaSeqs in a coordConvRep */ { struct hash *tFileCache = gfFileCacheNew(); if(!(ccr->midSeq && ccr->upSeq && ccr-> downSeq)) errAbort("coordConv::getAlignmentsForSeqs() - can't have any NULL dnaSeqs."); ccr->midPsl = doDnaAlignment(ccr->midSeq, ccr->to->version, blatHost, port, nibDir, tFileCache);