aa8e2e87c3b54e9662d87cd5b921d5607fc846ad max Mon Nov 10 05:51:08 2025 -0800 adding korea to varFreqs supertrack, refs #36642 diff --git src/hg/makeDb/doc/hg19.txt src/hg/makeDb/doc/hg19.txt index 05fd850a879..acba2869c66 100644 --- src/hg/makeDb/doc/hg19.txt +++ src/hg/makeDb/doc/hg19.txt @@ -35933,15 +35933,25 @@ # 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 +############################################################################# +# variant frequencies, 2025-11-07, Max +cd /hive/data/genomes/hg19/bed/varFreqs/mexbb/ +# got files as-is from Carmina +cd /hive/data/genomes/hg19/bed/varFreqs/sgdp/ +https://sharehost.hms.harvard.edu/genetics/reich_lab/sgdp/vcf_variants/vcfs.variants.public_samples.279samples.tar +tar xvfz *.tar +ls *.vcf.gz > vcf.list +bcftools merge -l vcf.list -O z -o SGDP.nh2.vcf.gz --threads 20 +tabix -p vcf SGDP.nh2.vcf.gz