c41282bd6cf17dd6d19072af594382ef7a33ee75
angie
Thu Jun 9 10:57:54 2011 -0700
Track #1656 (GWAS Catalog) update: not the usual 5-minute reload,because NHGRI interspersed a bunch of new columns into their download
file, and hg19 is now using snp132.
diff --git src/hg/hgc/hgc.c src/hg/hgc/hgc.c
index 456892f..b53ae33 100644
--- src/hg/hgc/hgc.c
+++ src/hg/hgc/hgc.c
@@ -17128,31 +17128,32 @@
boolean first = TRUE;
char **row;
while ((row = sqlNextRow(sr)) != NULL)
{
if (first)
first = FALSE;
else
printf("
\n");
struct gwasCatalog *gc = gwasCatalogLoad(row+rowOffset);
printCustomUrl(tdb, item, FALSE);
printPos(gc->chrom, gc->chromStart, gc->chromEnd, NULL, TRUE, gc->name);
printf("Reported region: %s
\n", gc->region);
printf("Publication: %s et al. "
"author);
printEntrezPubMedUidAbstractUrl(stdout, gc->pubMedID);
- printf("\" TARGET=_BLANK>%s. %s. %s
\n", gc->title, gc->journal, gc->pubDate);
+ printf("\" TARGET=_BLANK>%s%s %s. %s
\n",
+ gc->title, (endsWith(gc->title, ".") ? "" : "."), gc->journal, gc->pubDate);
printf("Disease or trait: %s
\n", subNrNs(gc->trait));
printf("Initial sample size: %s
\n", subNrNs(gc->initSample));
printf("Replication sample size: %s
\n", subNrNs(gc->replSample));
printf("Reported gene(s): %s
\n", subNrNs(gc->genes));
char *strongAllele = NULL, *strongRsID = splitSnpAndAllele(gc->riskAllele, &strongAllele);
if (strongRsID)
printf("Strongest SNP-Risk allele: "
"%s-%s
\n", strongRsID, strongRsID, strongAllele);
else
printf("Strongest SNP-Risk allele: %s
\n", subNrNs(gc->riskAllele));
gwasCatalogCheckSnpAlleles(tdb, gc);
printf("Risk Allele Frequency: %s
\n", subNrNs(gc->riskAlFreq));
if (isEmpty(gc->pValueDesc) || sameString(gc->pValueDesc, "NS"))
printf("p-Value: %s
\n", subNrNs(gc->pValue));