src/hg/makeDb/doc/hg18.txt 1.379

1.379 2009/09/03 15:39:14 hartera
Added a vegaPep table for the Vega Genes update.
Index: src/hg/makeDb/doc/hg18.txt
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/doc/hg18.txt,v
retrieving revision 1.378
retrieving revision 1.379
diff -b -B -U 4 -r1.378 -r1.379
--- src/hg/makeDb/doc/hg18.txt	23 Aug 2009 04:14:05 -0000	1.378
+++ src/hg/makeDb/doc/hg18.txt	3 Sep 2009 15:39:14 -0000	1.379
@@ -28626,8 +28626,10 @@
 # can be shared with Ensembl to create the links to Vega transcript, gene
 # and protein reports on the details pages. 
 # 2009-08-22 - Finished code for adding Vega report URLs to the details pages.
 # Loaded the vegaGtp table.
+# 2009-09-01 - 2009-09-02 (hartera). Loaded a vegaPep table for the protein
+# sequence link on the details pages.
 
     mkdir /hive/data/genomes/hg18/bed/vega35
     cd /hive/data/genomes/hg18/bed/vega35
     wget --timestamping "ftp://ftp.sanger.ac.uk/pub/vega/human/*" \
@@ -28698,6 +28700,23 @@
     hgLoadSqlTab hg18 vegaGtp ensGtp.sql vegaGtp.tab
     # Loaded succesfully
     # Added code to hgc.c to use printEnsemblOrVegaCustomUrl() in 
     # doVegaGene() to add the links to Vega reports on the details pages.
-
+    # Code was added so that there is no protein sequence link on the details
+    # page if it there is none available e.g. noncoding.
+    # 2009-09-01 - 2009-09-02 (hartera)
+    # Coding genes are displaying the message that there is no protein
+    # prediction available. Need to add a vegaPep table.
+    cd /hive/data/genomes/hg18/bed/vega35
+    # from the Ensembl process:
+    zcat Homo_sapiens.VEGA.mar.pep.tot.fa.gz  \
+        | sed -e 's/^>.* Transcript:/>/;' | gzip > vegaPep.txt.gz
+    zcat vegaPep.txt.gz \
+         | ~/kent/src/utils/faToTab/faToTab.pl /dev/null /dev/stdin \
+         | sed -e '/^$/d; s/*$//' | sort > vegaPep.hg18.fa.tab
+    # Load table
+    hgPepPred hg18 tab vegaPep vegaPep.hg18.fa.tab
+    # Add vegaPep to the trackDb.ra entry for the vegaGeneComposite track 
+    # in the type line for src/hg/makeDb/trackDb/human/hg18/trackDb.ra.
+    # Check that the vegaPep table looks ok and then check protein-coding and 
+    # noncoding transcript details pages for protein links.
 ############################################################################