5c6591502a68acc328ae00d0cbc31a61891fa9f9
jcasper
  Fri Aug 16 08:31:23 2024 -0700
Makedoc for knownGene VM36 on mm39, refs #34218

diff --git src/hg/makeDb/doc/ucscGenes/mm39.gencodeVM36.sh src/hg/makeDb/doc/ucscGenes/mm39.gencodeVM36.sh
new file mode 100644
index 0000000..bd49575
--- /dev/null
+++ src/hg/makeDb/doc/ucscGenes/mm39.gencodeVM36.sh
@@ -0,0 +1,52 @@
+# Jonathan - began 2024-08-14
+
+export db=mm39
+export GENCODE_VERSION=VM36
+export PREV_GENCODE_VERSION=VM35
+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
+