832926811b278d5c656957e1565acdde48522def
kate
Mon Apr 27 17:59:12 2020 -0700
Update GTEx gene track click handler and track settings to handle V8 (add Kidney Medulla tissue, support TPM and RPKM). refs #25130
diff --git src/hg/hgc/gtexClick.c src/hg/hgc/gtexClick.c
index 4ac2c33..5886fd5 100644
--- src/hg/hgc/gtexClick.c
+++ src/hg/hgc/gtexClick.c
@@ -107,33 +107,34 @@
else
{
printf("%s
\n",
hgGeneName(), database, gtexGene->name, gtexGene->name);
printf("Description: %s
\n", desc);
}
printf("Ensembl gene ID: %s
\n", gtexGene->geneId);
// The actual transcript model is a union, so this identification is approximate
// (used just to find a transcript class)
char *geneClass = gtexGeneClass(gtexGene);
printf("GENCODE biotype: %s
\n", gtexGene->geneType);
printf("Gene class: %s
\n",
geneClassColorCode(geneClass), geneClass);
int tisId;
float highLevel = gtexGeneHighestMedianExpression(gtexGene, &tisId);
-printf("Highest median expression: %0.2f RPKM in %s
\n",
- highLevel, gtexGetTissueDescription(tisId, version));
-printf("Total median expression: %0.2f RPKM
\n", gtexGeneTotalMedianExpression(gtexGene));
+printf("Highest median expression: %0.2f %s in %s
\n",
+ highLevel, gtexExprUnit(version), gtexGetTissueDescription(tisId, version));
+printf("Total median expression: %0.2f %s
\n", gtexGeneTotalMedianExpression(gtexGene),
+ gtexExprUnit(version));
printf("Score: %d
\n", gtexGene->score);
printf("Genomic position: "
"%s %s:%d-%d
\n",
database, hgTracksPathAndSettings(), database,
gtexGene->chrom, gtexGene->chromStart+1, gtexGene->chromEnd,
gtexGene->chrom, gtexGene->chromStart+1, gtexGene->chromEnd);
puts("
"); // set gtexDetails (e.g. to 'log') to show log transformed details page // if hgTracks is log-transformed boolean doLogTransform = (trackDbSetting(tdb, "gtexDetails") && cartUsualBooleanClosestToHome(cart, tdb, FALSE, GTEX_LOG_TRANSFORM, GTEX_LOG_TRANSFORM_DEFAULT)); struct tempName pngTn;