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);
sqlSafef(query, sizeof 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", getH1n1StructUrl(), 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, NULL, &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);
}
void doTransRegCode(struct trackDb *tdb, char *item, char *motifTable)
/* Display detailed info on a transcriptional regulatory code item. */
{
struct dnaMotif *motif = loadDnaMotif(item, motifTable);
int start = cartInt(cart, "o");
struct dnaSeq *seq = NULL;
char *table = tdb->table;
int rowOffset = hOffsetPastBin(database, seqName, table);
char query[256];
struct sqlResult *sr;
char **row;
struct sqlConnection *conn = hAllocConn(database);
struct transRegCode *trc = NULL;
cartWebStart(cart, database, "Regulatory Code Info");
sqlSafef(query, sizeof query,
"select * from %s where name = '%s' and chrom = '%s' and chromStart = %d",
table, item, seqName, start);
sr = sqlGetResult(conn, query);
row = sqlNextRow(sr);
if (row != NULL)
trc = transRegCodeLoad(row+rowOffset);
sqlFreeResult(&sr);
if (trc != NULL)
{
char strand[2];
seq = hDnaFromSeq(database, trc->chrom, trc->chromStart, trc->chromEnd, dnaLower);
if (seq->size != motif->columnCount)
{
printf("WARNING: seq->size = %d, motif->colCount=%d<BR>\n",
seq->size, motif->columnCount);
strand[0] = '?';
seq = NULL;
}
else
{
strand[0] = dnaMotifBestStrand(motif, seq->dna);
if (strand[0] == '-')
reverseComplement(seq->dna, seq->size);
}
strand[1] = 0;
printf("<B>Name:</B> ");
sacCerHgGeneLinkName(conn, trc->name);
printf("<BR>\n");
printf("<B>ChIP-chip Evidence:</B> %s<BR>\n", trc->chipEvidence);
printf("<B>Species conserved in:</B> %d of 2<BR>\n", trc->consSpecies);
if (seq != NULL)
printf("<B>Bit Score of Motif Hit:</B> %4.2f<BR>\n",
dnaMotifBitScore(motif, seq->dna));
printf("<B>Item score:</B> %d<BR>\n", trc->score);
printPosOnChrom(trc->chrom, trc->chromStart, trc->chromEnd, strand, TRUE, trc->name);
}
motifHitSection(seq, motif);
printTrackHtml(tdb);
}
static void singleBamDetails(const bam1_t *bam)
/* Print out the properties of this alignment. */
{
const bam1_core_t *core = &bam->core;
char *itemName = bam1_qname(bam);
int tLength = bamGetTargetLength(bam);
int tStart = core->pos, tEnd = tStart+tLength;
boolean isRc = useStrand && bamIsRc(bam);
printPosOnChrom(seqName, tStart, tEnd, NULL, FALSE, itemName);
if (!skipQualityScore)
printf("<B>Alignment Quality: </B>%d<BR>\n", core->qual);
printf("<B>CIGAR string: </B><tt>%s</tt> (", bamGetCigar(bam));
bamShowCigarEnglish(bam);
printf(")<BR>\n");
printf("<B>Tags:</B>");
bamShowTags(bam);
puts("<BR>");
printf("<B>Flags: </B><tt>0x%02x:</tt><BR>\n ", core->flag);
bamShowFlagsEnglish(bam);
puts("<BR>");
if (bamIsRc(bam))
printf("<em>Note: although the read was mapped to the reverse strand of the genome, "
"the sequence and CIGAR in BAM are relative to the forward strand.</em><BR>\n");
puts("<BR>");
struct dnaSeq *genoSeq = hChromSeq(database, seqName, tStart, tEnd);
char *qSeq = bamGetQuerySequence(bam, FALSE);
if (isNotEmpty(qSeq) && !sameString(qSeq, "*"))
{
char *qSeq = NULL;
struct ffAli *ffa = bamToFfAli(bam, genoSeq, tStart, useStrand, &qSeq);
printf("<B>Alignment of %s to %s:%d-%d%s:</B><BR>\n", itemName,
seqName, tStart+1, tEnd, (isRc ? " (reverse complemented)" : ""));
ffShowSideBySide(stdout, ffa, qSeq, 0, genoSeq->dna, tStart, tLength, 0, tLength, 8, isRc,
FALSE);
}
if (!skipQualityScore && core->l_qseq > 0)
{
printf("<B>Sequence quality scores:</B><BR>\n<TT><TABLE><TR>\n");
UBYTE *quals = bamGetQueryQuals(bam, useStrand);
int i;
for (i = 0; i < core->l_qseq; i++)
{
if (i > 0 && (i % 24) == 0)
printf("</TR>\n<TR>");
printf("<TD>%c<BR>%d</TD>", qSeq[i], quals[i]);
}
printf("</TR></TABLE></TT>\n");
}
}
static void vcfRecordDetails(struct trackDb *tdb, struct vcfRecord *rec)
/* Display the contents of a single line of VCF, assumed to be from seqName
* (using seqName instead of rec->chrom because rec->chrom might lack "chr"). */
{
printf("<B>Name:</B> %s<BR>\n", rec->name);
// Since these are variants, if it looks like a dbSNP or dbVar ID, provide a link:
if (regexMatch(rec->name, "^rs[0-9]+$"))
{
printf("<B>dbSNP:</B> ");
printDbSnpRsUrl(rec->name, "%s", rec->name);
puts("<BR>");
}
else if (regexMatch(rec->name, "^[en]ss?v[0-9]+$"))
{
printf("<B>dbVar:</B> ");
printf("<A HREF=\"http://www.ncbi.nlm.nih.gov/dbvar/variants/%s/\" "
"TARGET=_BLANK>%s</A><BR>\n", rec->name, rec->name);
}
printCustomUrl(tdb, rec->name, TRUE);
boolean hapClustEnabled = cartOrTdbBoolean(cart, tdb, VCF_HAP_ENABLED_VAR, TRUE);
if (hapClustEnabled)
{
static char *formName = "vcfCfgHapCenter";
printf("<FORM NAME=\"%s\" ACTION=\"%s\">\n", formName, hgTracksName());
cartSaveSession(cart);
vcfCfgHaplotypeCenter(cart, tdb, tdb->track, FALSE, rec->file, rec->name,
seqName, rec->chromStart, formName);
printf("</FORM>\n");
}
char leftBase = rec->alleles[0][0];
unsigned int vcfStart = vcfRecordTrimIndelLeftBase(rec);
boolean showLeftBase = (rec->chromStart == vcfStart+1);
(void)vcfRecordTrimAllelesRight(rec);
char *displayAls[rec->alleleCount];
makeDisplayAlleles(rec, showLeftBase, leftBase, 20, TRUE, FALSE, displayAls);
printPosOnChrom(seqName, rec->chromStart, rec->chromEnd, NULL, FALSE, rec->name);
printf("<B>Reference allele:</B> %s<BR>\n", displayAls[0]);
vcfAltAlleleDetails(rec, displayAls);
vcfQualDetails(rec);
vcfFilterDetails(rec);
vcfInfoDetails(rec);
pgSnpCodingDetail(rec);
makeDisplayAlleles(rec, showLeftBase, leftBase, 5, FALSE, TRUE, displayAls);
vcfGenotypesDetails(rec, tdb, displayAls);
}
static void bamPairDetails(const bam1_t *leftBam, const bam1_t *rightBam)
/* Print out details for paired-end reads. */
{
if (leftBam && rightBam)
{
const bam1_core_t *leftCore = &leftBam->core, *rightCore = &rightBam->core;
int leftLength = bamGetTargetLength(leftBam), rightLength = bamGetTargetLength(rightBam);
int start = min(leftCore->pos, rightCore->pos);
int end = max(leftCore->pos+leftLength, rightCore->pos+rightLength);
char *itemName = bam1_qname(leftBam);
printf("<B>Paired read name:</B> %s<BR>\n", itemName);
printPosOnChrom(seqName, start, end, NULL, FALSE, itemName);
puts("<P>");
}
showOverlap(leftBam, rightBam);
printf("<TABLE><TR><TD valign=top><H4>Left end read</H4>\n");
singleBamDetails(leftBam);
printf("</TD><TD valign=top><H4>Right end read</H4>\n");
singleBamDetails(rightBam);
printf("</TD></TR></TABLE>\n");
}
void doTransRegCodeProbe(struct trackDb *tdb, char *item,
char *codeTable, char *motifTable,
char *tfToConditionTable, char *conditionTable)
/* Display detailed info on a ChIP-chip probe from transRegCode experiments. */
{
char query[256];
struct sqlResult *sr;
char **row;
int rowOffset = hOffsetPastBin(database, seqName, tdb->table);
struct sqlConnection *conn = hAllocConn(database);
struct transRegCodeProbe *probe = NULL;
cartWebStart(cart, database, "ChIP-chip Probe Info");
sqlSafef(query, sizeof(query), "select * from %s where name = '%s'",
tdb->table, item);
sr = sqlGetResult(conn, query);
if ((row = sqlNextRow(sr)) != NULL)
probe = transRegCodeProbeLoad(row+rowOffset);
sqlFreeResult(&sr);
if (probe != NULL)
{
struct tfData *tfList = NULL, *tf;
struct hash *tfHash = newHash(0);
struct transRegCode *trc;
int i;
/* Print basic info. */
printf("<B>Name:</B> %s<BR>\n", probe->name);
printPosOnChrom(probe->chrom, probe->chromStart, probe->chromEnd,
NULL, TRUE, probe->name);
/* Make up list of all transcriptionFactors. */
for (i=0; i<probe->tfCount; ++i)
{
/* Parse out factor and condition. */
char *tfName = probe->tfList[i];
char *condition = strchr(tfName, '_');
struct tfCond *cond;
if (condition != NULL)
*condition++ = 0;
else
condition = "n/a";
tf = hashFindVal(tfHash, tfName);
if (tf == NULL)
{
AllocVar(tf);
hashAddSaveName(tfHash, tfName, tf, &tf->name);
slAddHead(&tfList, tf);
}
AllocVar(cond);
cond->name = cloneString(condition);
cond->binding = probe->bindVals[i];
slAddHead(&tf->conditionList, cond);
}
slSort(&tfList, tfDataCmpName);
/* Fold in motif hits in region. */
if (sqlTableExists(conn, codeTable))
{
sr = hRangeQuery(conn, codeTable,
probe->chrom, probe->chromStart, probe->chromEnd,
"chipEvidence != 'none'", &rowOffset);
while ((row = sqlNextRow(sr)) != NULL)
{
trc = transRegCodeLoad(row+rowOffset);
tf = hashFindVal(tfHash, trc->name);
if (tf != NULL)
slAddTail(&tf->trcList, trc);
}
sqlFreeResult(&sr);
}
if (tfList == NULL)
printf("No significant immunoprecipitation.");
else
{
tfBindLevelSection(tfList, conn, motifTable, tfToConditionTable);
}
transRegCodeProbeFree(&probe);
growthConditionSection(conn, conditionTable);
}
printf("\n<HR>\n");
printTrackHtml(tdb);
hFreeConn(&conn);
}
