ad92b87bb799c20abccadf1f923972c1a9b6a6e7 max Fri Oct 17 06:39:16 2025 -0700 adding gnomad mpc track, ref #36531 diff --git src/hg/makeDb/doc/hg19.txt src/hg/makeDb/doc/hg19.txt index f1bfe9b38f0..05fd850a879 100644 --- src/hg/makeDb/doc/hg19.txt +++ src/hg/makeDb/doc/hg19.txt @@ -35926,15 +35926,22 @@ # - a few protein-coding genes not in the assembly # - pseudogenes and non-coding genes not in GENCODE or RefSeq for unknown reasons # - locus_group `other' types that are not mapped. # put the files into /gbdb mkdir -p /gbdb/hg19/hgnc ln -sf $(realpath hgnc.bb search.ix search.ixx) /gbdb/hg19/hgnc/ # double check ls -l /gbdb/hg19/hgnc/ file /gbdb/hg19/hgnc/* # if not already done for hg38: # edit trackDb/human/trackDb.ra to set filterValues.locus_type # from the locus_type.filter.txt file + +############################################################################# +# gnomad MPC 2025-10-17 max +cd /hive/data/genomes/hg19/bed/gnomad/mpc +wget https://storage.googleapis.com/gcp-public-data--gnomad/release/2.1.1/regional_missense_constraint/gnomAD_v2.1.1_transcripts_with_rmc.tsv +cat gnomAD_v2.1.1_transcripts_with_rmc.tsv | python ~/kent/src/hg/makeDb/gnomad/gnomadMpcToBed.py | sort -k1,2n > mpc.bed +bedToBigBed mpc.bed ../../../chrom.sizes mpc.bb -as=~/kent/src/hg/makeDb/gnomad/mpc.as -type=bed9+ -tab