f5915923089694356c860f93b6ce42a863989edd braney Fri Jun 10 13:54:48 2016 -0700 get longRange hgc page to put out information about the feet #15232 diff --git src/hg/hgc/hgc.c src/hg/hgc/hgc.c index 985255b..84fb619 100644 --- src/hg/hgc/hgc.c +++ src/hg/hgc/hgc.c @@ -3956,47 +3956,60 @@ struct aveStats *as = aveStatsCalc(doubleArray, count); printf("Item you clicked on:<BR>\n"); printf("<B>Score:</B> %g<BR>\n", ourLongRange->score); printf("<B>ID:</B> %u<BR>\n", ourLongRange->id); unsigned padding = (ourLongRange->e - ourLongRange->s) / 10; int s = ourLongRange->s - padding; int e = ourLongRange->e + padding; if (s < 0 ) s = 0; int chromSize = hChromSize(database, seqName); if (e > chromSize) e = chromSize; char num1Buf[1024],num2Buf[1024]; -sprintLongWithCommas(num1Buf, ourLongRange->s); -sprintLongWithCommas(num2Buf, ourLongRange->e); +char num3Buf[1024],num4Buf[1024]; +char num5Buf[1024],num6Buf[1024]; +char num7Buf[1024]; +sprintLongWithCommas(num1Buf, ourLongRange->s - ourLongRange->sw/2); +sprintLongWithCommas(num2Buf, ourLongRange->s + ourLongRange->sw/2); +sprintLongWithCommas(num3Buf, ourLongRange->e - ourLongRange->ew/2); +sprintLongWithCommas(num4Buf, ourLongRange->e + ourLongRange->ew/2); +sprintLongWithCommas(num5Buf, ourLongRange->sw); +sprintLongWithCommas(num6Buf, ourLongRange->ew); +sprintLongWithCommas(num7Buf, ourLongRange->ew + ourLongRange->e - ourLongRange->s); if (differentString(ourLongRange->sChrom, ourLongRange->eChrom)) { - printf("<B>Interchromosomal: </B> (Position on chrom you clicked on)</B>"); - printf("<A HREF=\"hgTracks?position=%s:%s-%s\" TARGET=_BLANK><B>%s:%s-%s</B></A><BR>\n", - ourLongRange->sChrom, num1Buf,num1Buf,ourLongRange->sChrom, num1Buf,num1Buf); - printf("<B>Position on other chromosome:</B>"); - printf("<A HREF=\"hgTracks?position=%s:%s-%s\" TARGET=_BLANK><B>%s:%s-%s<BR></A><BR>\n", - ourLongRange->eChrom, num2Buf, num2Buf, ourLongRange->eChrom, num2Buf, num2Buf); + printf("<B>Current region: </B>"); + printf("<A HREF=\"hgTracks?position=%s:%s-%s \" TARGET=_BLANK>%s:%s-%s (%sbp)</A><BR>\n", + ourLongRange->sChrom, num1Buf,num2Buf,ourLongRange->sChrom, num1Buf,num2Buf, num5Buf); + printf("<B>Paired region: </B>"); + printf("<A HREF=\"hgTracks?position=%s:%s-%s \" TARGET=_BLANK>%s:%s-%s (%sbp)<BR></A><BR>\n", + ourLongRange->eChrom, num3Buf, num4Buf, ourLongRange->eChrom, num3Buf, num4Buf, num6Buf); } else { - char num3Buf[1024]; - sprintLongWithCommas(num3Buf, ourLongRange->e - ourLongRange->s); - printf("<B>Intrachromosomal:</B><A HREF=\"hgTracks?position=%s:%s-%s\" TARGET=_BLANK><B> %s:%s-%s (%sbp) </B></A><BR>\n", - ourLongRange->sChrom, num1Buf, num2Buf, ourLongRange->sChrom, num1Buf, num2Buf, num3Buf); + printf("<B>Lower region: </B>"); + printf("<A HREF=\"hgTracks?position=%s:%s-%s \" TARGET=_BLANK>%s:%s-%s (%sbp)</A><BR>\n", + ourLongRange->sChrom, num1Buf,num2Buf, ourLongRange->sChrom, num1Buf,num2Buf, num5Buf); + printf("<B>Upper region: </B>"); + printf("<A HREF=\"hgTracks?position=%s:%s-%s \" TARGET=_BLANK>%s:%s-%s (%sbp)<BR></A><BR>\n", + ourLongRange->eChrom, num3Buf, num4Buf, ourLongRange->eChrom, num3Buf, num4Buf, num6Buf); + printf("<B>Intrachromasomal interaction region: </B>"); + printf("<A HREF=\"hgTracks?position=%s:%s-%s \" TARGET=_BLANK>%s:%s-%s (%sbp)<BR></A><BR>\n", + ourLongRange->eChrom, num1Buf, num4Buf, ourLongRange->eChrom, num1Buf, num4Buf, num7Buf); } printf("<BR>Statistics on the scores of all items in window (go to track controls to set minimum score to display):\n"); printf("<TABLE BORDER=1>\n"); printf("<TR><TD><B>Q1</B></TD><TD>%f</TD></TR>\n", as->q1); printf("<TR><TD><B>median</B></TD><TD>%f</TD></TR>\n", as->median); printf("<TR><TD><B>Q3</B></TD><TD>%f</TD></TR>\n", as->q3); printf("<TR><TD><B>average</B></TD><TD>%f</TD></TR>\n", as->average); printf("<TR><TD><B>min</B></TD><TD>%f</TD></TR>\n", as->minVal); printf("<TR><TD><B>max</B></TD><TD>%f</TD></TR>\n", as->maxVal); printf("<TR><TD><B>count</B></TD><TD>%d</TD></TR>\n", as->count); printf("<TR><TD><B>total</B></TD><TD>%f</TD></TR>\n", as->total); printf("<TR><TD><B>standard deviation</B></TD><TD>%f</TD></TR>\n", as->stdDev); printf("</TABLE>\n");