0e016ca51069bf70ea5c9e08d7ef2fa4dd81a838 max Wed Feb 10 06:11:16 2021 -0800 adding makedoc for cadd track, refs #18492 diff --git src/hg/makeDb/doc/hg19.txt src/hg/makeDb/doc/hg19.txt index fd37668..08a2824 100644 --- src/hg/makeDb/doc/hg19.txt +++ src/hg/makeDb/doc/hg19.txt @@ -34365,15 +34365,29 @@ hgBbiDbLink hg19 fetalGeneAtlasCellTypes /gbdb/hg19/bbi/fetalGeneAtlas/cell_type.bb hgBbiDbLink hg19 fetalGeneAtlasDonor /gbdb/hg19/bbi/fetalGeneAtlas/donor.bb hgBbiDbLink hg19 fetalGeneAtlasAssay /gbdb/hg19/bbi/fetalGeneAtlas/Assay.bb hgBbiDbLink hg19 fetalGeneAtlasExperiment /gbdb/hg19/bbi/fetalGeneAtlas/Experiment_batch.bb hgBbiDbLink hg19 fetalGeneAtlasOrgan /gbdb/hg19/bbi/fetalGeneAtlas/Organ.bb hgBbiDbLink hg19 fetalGeneAtlasOrganCellLineage /gbdb/hg19/bbi/fetalGeneAtlas/Organ_cell_lineage.bb hgBbiDbLink hg19 fetalGeneAtlasRtGroup /gbdb/hg19/bbi/fetalGeneAtlas/RG_group.bb hgBbiDbLink hg19 fetalGeneAtlasSex /gbdb/hg19/bbi/fetalGeneAtlas/sex.bb # Add the bits from clust/*.ra and clust/*.colors to hg19/trackDb.ra and you should be good. foreach s (cell_type Assay Experiment_batch donor Organ Organ_cell_lineage RT_group sex) echo >> clust/$s.ra end cat clust/*.ra > tracks.ra +############################################################################# +# CADD, max Wed Feb 10 06:08:01 PST 2021 +cd /hive/data/genomes/hg19/bed/cadd/ +wget https://krishna.gs.washington.edu/download/CADD/v1.6/GRCh37/InDels.tsv.gz +wget https://krishna.gs.washington.edu/download/CADD/v1.6/GRCh37/whole_genome_SNVs.tsv.gz +time python ~/kent/src/hg/makeDb/cadd/caddToBed.py InDels.tsv.gz +time python ~/kent/src/hg/makeDb/cadd/caddToWig.py +wigToBigWig a.wig ../../chrom.sizes a.bw & +wigToBigWig c.wig ../../chrom.sizes c.bw & +wigToBigWig a.wig ../../chrom.sizes a.bw & +wigToBigWig g.wig ../../chrom.sizes g.bw & +bedToBigBed ins.bed ../../chrom.sizes ins.bb -type=bed9+ -tab -as=${HOME}/kent/src/hg/makeDb/cadd/cadd.as +bedToBigBed del.bed ../../chrom.sizes del.bb -type=bed9+ -tab -as=${HOME}/kent/src/hg/makeDb/cadd/cadd.as +rm -f *.wig *.bed