src/hg/hgc/virusClick.c 1.4

Index: src/hg/hgc/virusClick.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/hgc/virusClick.c,v
retrieving revision 1.3
retrieving revision 1.4
diff -b -B -U 4 -r1.3 -r1.4
--- src/hg/hgc/virusClick.c	14 Jun 2009 16:23:17 -0000	1.3
+++ src/hg/hgc/virusClick.c	15 Jun 2009 21:34:19 -0000	1.4
@@ -531,4 +531,63 @@
 sqlFreeResult(&sr);
 hFreeConn(&conn);
 }
 
+void doH1n1Gene(struct trackDb *tdb, char *item)
+/* Show details page for H1N1 Genes and Regions annotations track. */
+{
+struct sqlConnection *conn  = hAllocConn(database);
+struct sqlResult *sr;
+char query[256];
+char **row;
+char *chrom, *chromStart, *chromEnd;
+char *gene=NULL;
+
+genericHeader(tdb, item);
+
+gene = item;
+printf("<B>Gene: </B> %s\n<BR>", gene);
+sprintf(query, "select chrom, chromStart, chromEnd from h1n1Gene where name='%s';", gene);
+sr = sqlMustGetResult(conn, query);
+row = sqlNextRow(sr);
+if (row != NULL)
+   {
+   chrom      = row[0];
+   chromStart = row[1];
+   chromEnd   = row[2];
+   printPosOnChrom(chrom, atoi(chromStart), atoi(chromEnd), NULL, FALSE, item);
+   }
+sqlFreeResult(&sr);
+hFreeConn(&conn);
+htmlHorizontalLine();
+
+printf("<H3>Protein Structure Analysis and Prediction</H3>");
+printf("<B>3D Structure Prediction of consensus sequence (with variations of all selected sequences highlighted)");
+printf("<BR>PDB file:</B> ");
+
+char pdbUrl[PATH_LEN];
+safef(pdbUrl, sizeof(pdbUrl), "%s/%s/decoys/%s.try1-opt3.pdb.gz", h1n1StructUrl, item, item);
+
+// Modeller stuff
+char modelPdbUrl[PATH_LEN];
+if (getH1n1Model(gene, modelPdbUrl))
+    {
+    char *selectFile = cartOptionalString(cart, gisaidAaSeqList);
+    struct tempName imageFile, chimeraScript, chimerax;
+    mkH1n1StructData(gene, selectFile, &imageFile, &chimeraScript);
+    mkChimerax(gene, modelPdbUrl, chimeraScript.forCgi, &chimerax);
+    printf("<A HREF=\"%s\" TARGET=_blank>%s</A>, view with <A HREF=\"%s\">Chimera</A><BR>\n", 
+    	   modelPdbUrl, gene, chimerax.forHtml);
+    printf("<TABLE>\n");
+    printf("<TR>\n");
+    printf("<TD ALIGN=\"center\"><img src=\"%s\"></TD>", imageFile.forHtml);
+    printf("</TR>\n");
+    printf("</TABLE>\n");
+    }
+
+htmlHorizontalLine();
+printTrackHtml(tdb);
+
+sqlFreeResult(&sr);
+hFreeConn(&conn);
+}
+