8787ed4e68be3d7ab55b6673a88e8b83884ec682
kate
Fri Nov 1 17:25:41 2013 -0700
Expand motif info in details top panel. refs #9092
diff --git src/hg/hgc/peakClusters.c src/hg/hgc/peakClusters.c
index 3d1a190..88ac3cc 100644
--- src/hg/hgc/peakClusters.c
+++ src/hg/hgc/peakClusters.c
@@ -473,34 +473,42 @@
{
struct dnaSeq *seq = hDnaFromSeq(database,
seqName, hit->chromStart, hit->chromEnd, dnaLower);
if(hit->strand[0] == '-')
reverseComplement(seq->dna, seq->size);
seqs[i] = seq;
// TODO: move to hgc.c (with other pos printers)
safef(posLink, sizeof(posLink),"%s:%d-%d",
hgTracksPathAndSettings(), database,
cluster->chrom, hit->chromStart+1, hit->chromEnd,
cluster->chrom, hit->chromStart+1, hit->chromEnd);
printf("Motif Score");
if (hitCount > 1)
+ {
printf("#%d", i + 1);
- printf(": %.2f (%s max: %.2f) at %s %c
",
+ printf(": %.2f (%s max: %.2f) at %s on %c strand
",
hit->score, cluster->name, maxScore, posLink, (int)hit->strand[0]);
}
+ else
+ {
+ printf(": %.2f (%s max: %.2f)
\n", hit->score, cluster->name, maxScore);
+ printf("Motif Position: %s
\n", posLink);
+ printf("Motif Strand: %c
\n", (int)hit->strand[0]);
+ }
+ }
}
if (seqs != NULL)
{
motifMultipleHitsSection(seqs, hitCount, motif);
}
/* Get list of tracks we'll look through for input. */
char *inputTrackTable = trackDbRequiredSetting(tdb, "inputTrackTable");
sqlSafef(query, sizeof(query),
"select tableName from %s where factor='%s' order by source", inputTrackTable,
cluster->name);
/* Next do the lists of hits and misses. We have the hits from the non-zero signals in
* cluster->expScores. We need to figure out the sources actually assayed though