6433773a86994d104ca22309fefea208573c29d4 max Fri Apr 16 08:42:26 2021 -0700 fixing terrible trackDb typo that lead to cadd being dropped on the RR at some point, no ticket, automatic email from trackDb checker diff --git src/hg/makeDb/doc/hg19.txt src/hg/makeDb/doc/hg19.txt index 86a0199..78a48fd 100644 --- src/hg/makeDb/doc/hg19.txt +++ src/hg/makeDb/doc/hg19.txt @@ -34496,15 +34496,29 @@ | CDBox | 269 | | HAcaBox | 112 | | miRNA | 1918 | | scaRna | 21 | +---------+------+ # backup old table and compare to new one before rename. Should be only miRNA update: hgsql -Ne "select * from wgRna" hg19 | cut -f2- > wgRna.backup comm -23 <(sort -k1 -k2n wgRna2021-02-18.bed) <(sort -k1 -k2n wgRna.backup) | cut -f9 | sort -u miRNA # rename wgRnaNew table: hgsqlSwapTables -dropTable3 hg19 wgRnaNew wgRna wgRnaOld ############################################################################## +# Add Revel track, Max, Thu Apr 15 03:33:49 PDT 2021 +cd /hive/data/genomes/hg19/bed/revel +aria2c https://rothsj06.u.hpc.mssm.edu/revel_grch38_all_chromosomes.csv.zip +unzip revel_grch38_all_chromosomes.csv.zip +# make sure that no weird chroms and no weird multi-nucleotide changes +cat revel_grch38_all_chromosomes.csv | tr ',' '\t' | cut -f5 | grep -v alt | awk '{print length($1)}' | uniq > lens.tmp & +cut -d, revel_grch38_all_chromosomes.csv -f1 | uniq | sort | uniq -c > chroms.tmp & +pigz revel_grch38_all_chromosomes.csv +time python revelToWig.py hg19 +wigToBigWig t.wig ../chrom.sizes t.bw & +wigToBigWig c.wig ../chrom.sizes c.bw & +wigToBigWig a.wig ../chrom.sizes a.bw & +wigToBigWig g.wig ../chrom.sizes g.bw & +##############################################################################