cc84507a07203d9cf4aeec4b6335e4f0d9b93216 jcasper Mon Mar 3 10:45:32 2025 -0800 Makedoc for Gencode VM37 on mm39, refs #35819 diff --git src/hg/makeDb/doc/ucscGenes/mm39.gencodeVM37.sh src/hg/makeDb/doc/ucscGenes/mm39.gencodeVM37.sh new file mode 100644 index 00000000000..51b57fe0f24 --- /dev/null +++ src/hg/makeDb/doc/ucscGenes/mm39.gencodeVM37.sh @@ -0,0 +1,52 @@ +# Jonathan - began 2025-03-01 + +export db=mm39 +export GENCODE_VERSION=VM37 +export PREV_GENCODE_VERSION=VM36 +screen -S knownGene${GENCODE_VERSION} +mkdir /hive/data/genomes/$db/bed/gencode$GENCODE_VERSION/build +cd /hive/data/genomes/$db/bed/gencode$GENCODE_VERSION/build + +PATH=$HOME/kent/src/hg/utils/otto/knownGene":$PATH" +cp /hive/data/genomes/${db}/bed/gencode${PREV_GENCODE_VERSION}/build/buildEnv.sh buildEnv.sh + +# edit buildEnv.sh + . buildEnv.sh + +cp ${oldGeneDir}/${PREV_GENCODE_VERSION}.files.txt . + +cp ${oldGeneDir}/${PREV_GENCODE_VERSION}.tables.txt . + +hgsql ${oldKnownDb} -Ne "show tables" > ${oldKnownDb}.tables.txt +diff <(sort ${PREV_GENCODE_VERSION}.tables.txt) <(sort ${oldKnownDb}.tables.txt) +# no difference + +buildKnown.sh & +# wait for completion + +tail -n 1 *.log +# ==> doBioCyc.log <== +# BuildBioCyc successfully finished +# +# ==> doBlast.log <== +# BuildBlast successfully finished +# +# ==> doFoldUtr.log <== +# BuildFoldUtr successfully finished +# +# ==> doKnown.log <== +# BuildKnown successfully finished +# +# ==> doKnownTo.log <== +# BuildKnownTo successfully finished +# +# ==> doPfamScop.log <== +# BuildPfamScop successfully finished + +# Manual GTF step +cd /hive/data/genomes/$db/goldenPath/bigZips/genes +genePredToGtf -utr ${tempDb} knownGene ${db}.knownGene.gtf +rm -f ${db}.knownGene.gtf.gz +gzip ${db}.knownGene.gtf +cd $dir +