36f8f6fb024b20cc523cdf9ebde7491eca84fd7c markd Sun Dec 6 20:33:20 2020 -0800 multiple request per connect works except hgBlat diff --git src/hg/lib/coordConv.c src/hg/lib/coordConv.c index f7e1675..7764dea 100644 --- src/hg/lib/coordConv.c +++ src/hg/lib/coordConv.c @@ -253,35 +253,35 @@ 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 */ -struct gfConnection *conn = gfConnect(blatHost, port); +struct gfConnection *conn = gfConnect(blatHost, port, FALSE); gfAlignStrand(conn, nibDir, seq, FALSE, 20, tFileCache, gvo, NULL, NULL); reverseComplement(seq->dna, seq->size); gfDisconnect(&conn); -conn = gfConnect(blatHost, port); +conn = gfConnect(blatHost, port, FALSE); gfAlignStrand(conn, nibDir, seq, TRUE, 20 , tFileCache, gvo, NULL, NULL); gfDisconnect(&conn); 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))