b63c6cd3032d2395c889cb0a58b28464d9d7f6a6 braney Tue Mar 17 09:01:34 2020 -0700 fix some of the problems from #25160 diff --git src/hg/makeDb/doc/ucscGenes/hg38.ucscGenes20.sh src/hg/makeDb/doc/ucscGenes/hg38.ucscGenes20.sh index 8460609..f6bb475 100755 --- src/hg/makeDb/doc/ucscGenes/hg38.ucscGenes20.sh +++ src/hg/makeDb/doc/ucscGenes/hg38.ucscGenes20.sh @@ -1124,15 +1124,23 @@ ln -s `pwd`/$gp.$mz.exonAA.fa.gz $pd/$gp.exonAA.fa.gz ln -s `pwd`/$gp.$mz.exonNuc.fa.gz $pd/$gp.exonNuc.fa.gz cd $pd md5sum *.fa.gz > md5sum.txt #myGene2 mkdir $dir/myGene2 cd $dir/myGene2 # copy list of genes from https://mygene2.org/MyGene2/genes awk '{print $1}' | sort > genes.lst hgsql hg38 -Ne "select geneSymbol, kgId from kgXref" | sort > ids.txt join genes.lst ids.txt | awk '{print $2,$1}' | sort > knownToMyGene2.txt hgLoadSqlTab $db knownToMyGene2 ~/kent/src/hg/lib/knownTo.sql knownToMyGene2.txt + +# make gtexDistance table +mkdir $dir/hgNear +cd $dir/hgNear +hgsql hgFixed -e select 1.0, id, name from gtexTissue > gtex.weights +hgMapToGene hg38 -all -type=genePred gtexGeneModelV6 knownGene knownToGtex +hgExpDistance hg38 -verbose=2 -lookup=knownToGtex -weights=gtex.weights hgFixed.gtexTissueMedian hgFixed.gtexTissue gtexDistance +