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.
\n"
+ "The approximate results below are sorted by number of matching 'tiles', "
"perfectly matching sub-sequences of length 11 (DNA) "
"or 4 (protein).
");
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.
\n");
+ printf("If its alignment score is < 20, including gaps and mismatches, 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;