src/hg/makeDb/doc/hg19.txt 1.108
1.108 2010/05/12 23:46:11 angie
Updated gwasCatalog.
Index: src/hg/makeDb/doc/hg19.txt
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/doc/hg19.txt,v
retrieving revision 1.107
retrieving revision 1.108
diff -b -B -U 4 -r1.107 -r1.108
--- src/hg/makeDb/doc/hg19.txt 12 May 2010 23:34:42 -0000 1.107
+++ src/hg/makeDb/doc/hg19.txt 12 May 2010 23:46:11 -0000 1.108
@@ -8357,27 +8357,29 @@
featureBits hg19 vegaPseudoGene
# 6885145 bases of 2897316137 (0.238%) in intersection
########################################################################
-# NHGRI GWAS CATALOG (DONE 4/1/10 angie)
+# NHGRI GWAS CATALOG (DONE 5/12/10 angie)
+# Updated 4/1/10
# Updated 3/1/10
# Originally done 1/19/10
mkdir /hive/data/genomes/hg19/bed/gwasCatalog
cd /hive/data/genomes/hg19/bed/gwasCatalog
+ # Done once, don't need to redo:
cut -f 1-4 ../snp130/snp130.bed \
| sort -k4,4 \
> snp130Coords.bed
- mkdir /hive/data/genomes/hg19/bed/gwasCatalog/100401
- cd /hive/data/genomes/hg19/bed/gwasCatalog/100401
+ mkdir /hive/data/genomes/hg19/bed/gwasCatalog/100512
+ cd /hive/data/genomes/hg19/bed/gwasCatalog/100512
# Mapping to hg19 by joining hg19 SNP coords with catalog flatfile (see hg18.txt)
- join -t " " -1 4 ../snp130Coords.bed /hive/data/genomes/hg18/bed/gwasCatalog/100401/noCoords.tab \
+ join -t " " -1 4 ../snp130Coords.bed /hive/data/genomes/hg18/bed/gwasCatalog/100512/noCoords.tab \
-o 1.1,1.2,1.3,1.4,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,2.10,2.11,2.12,2.13,2.14,2.15,2.16,2.17,2.18,2.19 \
| sort -k1,1 -k2n,2n \
> gwasCatalog.bed
hgLoadBed hg19 gwasCatalog gwasCatalog.bed \
-tab -sqlTable=$HOME/kent/src/hg/lib/gwasCatalog.sql -notItemRgb -allowStartEqualEnd
-#Loaded 2895 elements of size 22
+#Loaded 3277 elements of size 22
########################################################################
# ailMel1 Panda alignment (DONE - 2010-02-04 - Hiram)