src/hg/makeDb/doc/hg18.txt 1.414
1.414 2010/05/05 00:49:41 kent
Starting to put in instructions for the encode regulatory track. Will need polishing, but at least this will help me remember what I did.....
Index: src/hg/makeDb/doc/hg18.txt
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/doc/hg18.txt,v
retrieving revision 1.413
retrieving revision 1.414
diff -b -B -U 4 -r1.413 -r1.414
--- src/hg/makeDb/doc/hg18.txt 28 Apr 2010 22:35:21 -0000 1.413
+++ src/hg/makeDb/doc/hg18.txt 5 May 2010 00:49:41 -0000 1.414
@@ -16383,8 +16383,38 @@
see file: ucscGenes10.txt
#######################################################################
+# ENCODE Regulation track (In progress May 2010, very rough instructions
+# as am still defining how to build track...)
+
+ # Create the DNAse peak clusters subtrack.
+ # Get all of the narrowPeak format files for the wgEncodeUwDnaseSeq
+ # unzipped in directory /hive/users/kent/regulate/dnase/peaks
+ cd /hive/users/kent/regulate/dnase
+ ls -1 peaks/*.narrowPeak > peak.lst
+ regClusterMakeTableOfTables peak.lst peak.table
+ regCluster peak.table /dev/null peak.bed
+ awk '$4 > 1 || $5 >= 100' peak.bed > wgEncodeRegDnaseClustered.bed
+ hgLoadBed hg18 wgEncodeRegDnaseClustered wgEncodeRegDnaseClustered.bed
+
+ # Create the Transcription Factor Binding Site subtrack
+ # Get all of the narrowPeak file for wgEncodeYaleChIPseq into tfbs/yale
+ # directory and create all.lst by cutting and pasting
+ # the file list from
+ # http://hgdownload.cse.ucsc.edu/goldenPath/hg18/encodeDCC/wgEncodeHudsonalphaChipSeq/
+ cd /hive/users/kent/regulate/tfbs/yale
+ grep narrowPeak all. > peak.lst
+ cut -f 2 peak.lst > file.lst
+ foreach i (`cat file.lst`)
+ wget http://hgdownload.cse.ucsc.edu/goldenPath/hg18/encodeDCC/wgEncodeHudsonalphaChipSeq/$i
+ end
+ # Then edit peak.cfg to make sure all cell lines have uniq letters
+ regClusterBedExpCfg peak.lst -encodeList peak.cfg peak.bed peak.exps
+ hgLoadBed hg18 wgEncodeRegTfbsClustered peak.bed
+ hgLoadSqlTab hg18 wgEncodeRegTfbsCells ~/kent/src/hg/lib/expRecord.sql peak.exps
+
+#######################################################################
# CGAP SAGE (DONE 2007-04-17 Andy)
ssh hgwdev
bash