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");