96180f7e8f52b8b33e3f6e960839f55bcaac4665
hiram
Tue Feb 9 13:51:27 2016 -0800
need to make hte protein click links for the ncbiRefSeq resf #13673
diff --git src/hg/hgc/hgc.c src/hg/hgc/hgc.c
index 81a481b..befa85c 100644
--- src/hg/hgc/hgc.c
+++ src/hg/hgc/hgc.c
@@ -11508,30 +11508,49 @@
if (pslList)
{
int start = cartInt(cart, "o");
printf("
mRNA/Genomic Alignments
");
printAlignments(pslList, start, "ncbiRefSeqPsl", "ncbiRefSeqPsl", itemName);
}
else
{
printf ("there is NO alignment for %s
\n", itemName);
}
htmlHorizontalLine();
showGenePos(itemName, tdb);
+printf("Links xyz to sequence:
\n");
+printf("\n");
+printf("- need to figure out how to make these links, these don't work.\n");
+
+// might need something like this: seq = hPepSeq(database, buffer);
+
+puts("
- \n");
+hgcAnchorSomewhere("htcTranslatedProtein", itemName, "seqNcbiRefSeq",
+ nrl->mrnaAcc);
+printf("Predicted Protein \n");
+puts("
\n");
+puts("- \n");
+hgcAnchorSomewhere("htcGeneMrna", itemName, tdb->track, nrl->mrnaAcc);
+printf("%s may be different from the genomic sequence.\n",
+ "Predicted mRNA");
+puts("
\n");
+
+printf("
\n");
+
printTrackHtml(tdb);
hFreeConn(&conn);
}
void doRefGene(struct trackDb *tdb, char *rnaName)
/* Process click on a known RefSeq gene. */
{
struct sqlConnection *conn = hAllocConn(database);
int start = cartInt(cart, "o");
int left = cartInt(cart, "l");
int right = cartInt(cart, "r");
char *chrom = cartString(cart, "c");
boolean isXeno = sameString(tdb->table, "xenoRefGene");
if (isXeno)