b8d36ce63fd4a9803ad058357a8aaf16529837ed max Tue Jun 25 02:42:41 2024 -0700 Setting back the hgBlat internal minscore, for some reason it had no effect on hg38, worked only on hg18. refs #32918 diff --git src/hg/hgBlat/hgBlat.c src/hg/hgBlat/hgBlat.c index 5e4d2fa..7dafa61 100644 --- src/hg/hgBlat/hgBlat.c +++ src/hg/hgBlat/hgBlat.c @@ -1349,32 +1349,32 @@ int alphaBetSize; if (isProt) { alphaBetSize = 20; genomeSize = genomeSize / 3; } else { alphaBetSize = 4; } int k = 1; double expected = genomeSize; for (k=1; k<36; k++) { expected /= alphaBetSize; - // set this to .05 to allow 18bp searches on hg38. - if (expected < .19) + // Set this to .19 to allow 18bp searches on hg18 (no effect on hg38?). + if (expected < .004) break; } return k; } long findGenomeSize(char *database) // get genomeSize from database. { struct sqlConnection *conn = hAllocConn(database); char query[256]; sqlSafef(query, sizeof query, "select sum(size) from chromInfo"); long genomeSize = sqlQuickLongLong(conn, query); hFreeConn(&conn); if (genomeSize == 0) {