a9947a9b812bb9ac68e6758a7976e00dae5f3aec galt Thu May 2 14:46:13 2019 -0700 Further improving hgBlat blat all genomes explanation. diff --git src/hg/hgBlat/hgBlat.c src/hg/hgBlat/hgBlat.c index a13fc42..2b520da 100644 --- src/hg/hgBlat/hgBlat.c +++ src/hg/hgBlat/hgBlat.c @@ -1924,33 +1924,33 @@ { 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 single best alignment found for each assembly is shown below.
\n" "The approximate results below are sorted by number of matching 'tiles', " "perfectly matching sub-sequences of length 11 (DNA) " - "or 4 (protein).
"); + "or 4 (protein). Using only tile hits, this speedy method can not see mismatches.
"); printf("Click the 'assembly' link to trigger a full BLAT alignment for that genome. \n"); - printf("If its alignment score is < 20, including gaps and mismatches, no match will be found.
\n"); + printf("If its alignment score which includes gaps and mismatches is less than 20, no match will be found.
\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("

\n"); } lastSeqNumber = gH->seqNumber;