76e480f51fbf36171248a26c9f24c3bd0e95af7f kate Thu Jun 4 11:38:18 2020 -0700 Save build script and schema for CCRE track. refs #24668 diff --git src/hg/makeDb/doc/hg38/reg.txt src/hg/makeDb/doc/hg38/reg.txt index 697ad22..75aa882 100644 --- src/hg/makeDb/doc/hg38/reg.txt +++ src/hg/makeDb/doc/hg38/reg.txt @@ -663,21 +663,23 @@ 10.000000 2 # noting that order of accessions in score file doesn't match bed file ;-( sort Human-maxZ-DNase.txt > Human-maxZ-DNase.sorted.txt sort -k 4 GRCh38-ccREs.bed > GRCh38-ccREs.sorted.bed paste GRCh38-ccREs.sorted.bed Human-maxZ-DNase.sorted.txt > ccres.prescored.bed # manually sanity check start and end of file to see that accessions match # score using zscore, min(zscore*100),1000), and reformat awk '{OFS="\t"; print $1, $2, $3, $4, ($13>10)? 1000 : int($13 * 100), $6, $7, $8, $9, $10, $13}' \ ccres.prescored.bed | bedSort stdin ccres.scored.bed # Reformat to add fields for filtering and mouseover, etc. set f = encodeCcreCombined -perl makeCcreCombined.pl < ccres.scored.bed > $f.bed -bedToBigBed -tab -type=bed9+6 -as=$f.as $f.bed /hive/data/genomes/hg38/chrom.sizes $f.bb +set bin = ~/kent/src/hg/makeDb/outside/encode3/ccre +perl $bin/makeCcreCombined.pl < ccres.scored.bed > $f.bed +set lib = ~/kent/src/hg/lib +bedToBigBed -tab -type=bed9+6 -as=$lib/$f.as $f.bed /hive/data/genomes/hg38/chrom.sizes $f.bb ln -s `pwd`/$f.bb /gbdb/hg38/encode3/ccre/ ###############################