d20789efcf6e484753b9c8a21b1dd8db40961215
jcasper
Wed May 9 12:45:38 2018 -0700
Initial commit of gnomAD track on hg19, refs #18951
diff --git src/hg/hgc/vcfClick.c src/hg/hgc/vcfClick.c
index 16cb0f5..aaf91b3 100644
--- src/hg/hgc/vcfClick.c
+++ src/hg/hgc/vcfClick.c
@@ -474,30 +474,38 @@
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"). */
{
if (isNotEmpty(rec->name) && differentString(rec->name, "."))
printf("Name: %s
\n", rec->name);
if (sameString(tdb->track, "exacVariants"))
{
printf("ExAC: "
"%s:%d %s/%s
\n",
skipChr(rec->chrom), rec->chromStart+1, rec->alleles[0], rec->alleles[1],
skipChr(rec->chrom), rec->chromStart+1, rec->alleles[0], rec->alleles[1]);
}
+if (sameString(tdb->track, "gnomadVariants"))
+ {
+ printf("gnomAD: "
+ "%s:%d %s/%s
\n",
+ skipChr(rec->chrom), rec->chromStart+1, rec->alleles[0], rec->alleles[1],
+ skipChr(rec->chrom), rec->chromStart+1, rec->alleles[0], rec->alleles[1]);
+ }
// Since these are variants, if it looks like a dbSNP or dbVar ID, provide a link:
if (regexMatch(rec->name, "^rs[0-9]+$"))
{
printf("dbSNP: ");
printDbSnpRsUrl(rec->name, "%s", rec->name);
puts("
");
}
else if (regexMatch(rec->name, "^[en]ss?v[0-9]+$"))
{
printf("dbVar: ");
printf("%s
\n", rec->name, rec->name);
}
printCustomUrl(tdb, rec->name, TRUE);
boolean hapClustEnabled = cartOrTdbBoolean(cart, tdb, VCF_HAP_ENABLED_VAR, TRUE);