3789b789440d37aa713c3160a8686fe30b59ee66 braney Fri Mar 16 17:50:03 2012 -0700 fix #5457, properly treat non-coding refSeq genes by printing out the "there is no protein here" message diff --git src/hg/hgc/hgc.c src/hg/hgc/hgc.c index 2a07d79..4a52b49 100644 --- src/hg/hgc/hgc.c +++ src/hg/hgc/hgc.c @@ -2616,35 +2616,36 @@ { puts("
  • \n"); hgcAnchorSomewhere(pepClick, pepName, pepTable, seqName); printf("Predicted Protein \n"); puts("
  • \n"); foundPep = TRUE; } if (!foundPep) { char *autoTranslate = trackDbSetting(tdb, "autoTranslate"); if (autoTranslate == NULL || differentString(autoTranslate, "0")) { puts("
  • \n"); /* put out correct message to describe translated mRNA */ if ((sameString(geneTable, "ensGene")) || (sameString(geneTable, "vegaGene")) || (sameString(geneTable, "vegaPseudoGene")) + || (genbankIsRefSeqNonCodingMRnaAcc(geneName)) || (sameString(geneTable, "lincRNAsTranscripts")) ) { printf("Non-protein coding gene or gene fragment, no protein prediction available."); } - else if (!genbankIsRefSeqNonCodingMRnaAcc(geneName)) + else { hgcAnchorSomewhere("htcTranslatedPredMRna", geneName, "translate", seqName); printf("Translated Protein from "); if (sameString(geneTable, "refGene") ) { printf("genomic DNA\n"); } else { printf("predicted mRNA \n"); } foundPep = TRUE; } puts("
  • \n");