128526c4704a9a40aa4dd26e4bdb3662c4daebb0
braney
Mon Nov 28 17:57:25 2016 -0800
for BLAT custom track support, call bedToBigBed in a way that disables its checking for overlapping
exons, which BLAT is generating for some protein queries
diff --git src/hg/hgBlat/hgBlat.c src/hg/hgBlat/hgBlat.c
index ea8c389..f3d9b91 100644
--- src/hg/hgBlat/hgBlat.c
+++ src/hg/hgBlat/hgBlat.c
@@ -260,31 +260,31 @@
printf("", url);
printf("", url);
htmlEnd();
}
static void makeBigPsl(char *pslName, char *faName, char *db, char *outputBigBed)
/* Make a bigPsl with the blat results. */
{
struct tempName bigPslTn;
trashDirFile(&bigPslTn, "hgBlat", "bp", ".bigPsl");
char cmdBuffer[4096];
safef(cmdBuffer, sizeof(cmdBuffer), "loader/pslToBigPsl %s -fa=%s stdout | sort -k1,1 -k2,2n > %s", pslName, faName, bigPslTn.forCgi);
system(cmdBuffer);
-safef(cmdBuffer, sizeof(cmdBuffer), "loader/bedToBigBed -extraIndex=name -tab -as=loader/bigPsl.as -type=bed12+13 %s http://hgdownload.cse.ucsc.edu/goldenPath/%s/bigZips/%s.chrom.sizes %s",
+safef(cmdBuffer, sizeof(cmdBuffer), "loader/bedToBigBed -extraIndex=name -tab -as=loader/bigPsl.as -type=bed9+16 %s http://hgdownload.cse.ucsc.edu/goldenPath/%s/bigZips/%s.chrom.sizes %s",
bigPslTn.forCgi, db, db, outputBigBed);
system(cmdBuffer);
unlink(bigPslTn.forCgi);
}
void showAliPlaces(char *pslName, char *faName, char *customText, char *database,
enum gfType qType, enum gfType tType,
char *organism, boolean feelingLucky)
/* Show all the places that align. */
{
struct lineFile *lf = pslFileOpen(pslName);
struct psl *pslList = NULL, *psl;
char *browserUrl = hgTracksName();
char *hgcUrl = hgcName();
char uiState[64];