88483f54e928d2341e942d5f617ea935b8e5dbf4 galt Wed May 1 13:06:59 2019 -0700 Small improvement in hgBlat All Genomes output explanation in the header. diff --git src/hg/hgBlat/hgBlat.c src/hg/hgBlat/hgBlat.c index 4f40209..8b15a2d 100644 --- src/hg/hgBlat/hgBlat.c +++ src/hg/hgBlat/hgBlat.c @@ -1921,35 +1921,36 @@ // Combine pairs with a query and its RC. if (!(debuggingGfResults) && (sameString(pfdDone->xType,"dna") || sameString(pfdDone->xType,"dnax"))) { slSort(&pfdDone, rcPairsCmp); hideWeakerOfQueryRcPairs(pfdDone); } // requires db for chromSize, do database after multi-threading done. changeMaxGenePositionToPositiveStrandCoords(pfdDone); // sort by maximum hits slSort(&pfdDone, genomeHitsCmp); // Print instructions - printf("The approximate results below are sorted by number of matching 'tiles', " + printf("The single best alignment found for each assembly is shown below.<br>\n" + "The approximate results below are sorted by number of matching 'tiles', " "perfectly matching sub-sequences of length 11 (DNA) " "or 4 (protein). <br>"); printf("Click the 'assembly' link to trigger a full BLAT alignment for that genome. \n"); - printf("If its alignment score is < 20 bp, no match will be found.<br>\n"); + printf("If its alignment score is < 20, including gaps and mismatches, no match will be found.<br>\n"); // Print report // TODO move to final report at the end of ALL Assemblies int lastSeqNumber = -1; int idCount = 0; char id[256]; struct genomeHits *gH = NULL; for (gH = pfdDone; gH; gH = gH->next) { if (lastSeqNumber != gH->seqNumber) { if (lastSeqNumber != -1) // end previous table { printf("</TABLE><br><br>\n"); } lastSeqNumber = gH->seqNumber;