9a2f1da0c665b5b5afb248b89e2a18418f993aac hiram Thu May 13 13:15:15 2021 -0700 better naming for archived versions refs #24735 diff --git src/hg/utils/automation/doNcbiRefSeq.pl src/hg/utils/automation/doNcbiRefSeq.pl index 1cbc3d0..2f93677 100755 --- src/hg/utils/automation/doNcbiRefSeq.pl +++ src/hg/utils/automation/doNcbiRefSeq.pl @@ -559,35 +559,41 @@ &HgAutomate::mustMkdir($runDir); my $whatItDoes = "load NCBI processed files into track files."; my $bossScript = newBash HgRemoteScript("$runDir/doLoad.bash", $dbHost, $runDir, $whatItDoes); my $gbdbDir = "$HgAutomate::gbdb/\$db/ncbiRefSeq"; my $dbTwoBit = "$HgAutomate::clusterData/$db/$db.2bit"; $dbTwoBit = $target2bit if (-s "$target2bit"); my $genePredCheckDb = "genePredCheck -db=\$db"; if (! $dbExists) { $genePredCheckDb = "genePredCheck"; } + my $verString = `cat $buildDir/process/ncbiRefSeqVersion.txt`; + chomp $verString; + $verString =~ s/.*elease //; + $verString =~ s/ .*//; + $bossScript->add(<<_EOF_ # establish all variables to use here -export db=$db -export asmId=$asmId +export db="$db" +export asmId="$asmId" +export verString="$verString" _EOF_ ); if (! $dbExists) { $bossScript->add(<<_EOF_ export target2bit=$dbTwoBit twoBitInfo \$target2bit stdout | sort -k2,2nr > \$db.chrom.sizes wget -O bigGenePred.as 'http://genome-source.soe.ucsc.edu/gitlist/kent.git/raw/master/src/hg/lib/bigGenePred.as' wget -O bigPsl.as 'http://genome-source.soe.ucsc.edu/gitlist/kent.git/raw/master/src/hg/lib/bigPsl.as' ### overall gene track with both predicted and curated genePredToBigGenePred process/\$db.ncbiRefSeq.gp stdout | sort -k1,1 -k2,2n > \$db.ncbiRefSeq.bigGp genePredToBed -tab -fillSpace process/\$db.ncbiRefSeq.gp stdout \\ | bedToExons stdin stdout | bedSingleCover.pl stdin > \$asmId.exons.bed @@ -851,35 +857,34 @@ | pslMismatchGapToBed -cdsFile=stdin -db=\$db -ignoreQNamePrefix=X \\ process/\$asmId.\$db.psl.gz $dbTwoBit \\ \$db.rna.fa ncbiRefSeqGenomicDiff || true rm -f $gbdbDir/ncbiRefSeqGenomicDiff.bb if [ -s ncbiRefSeqGenomicDiff.bed ]; then bedToBigBed -type=bed9+ -tab -as=\${HOME}/kent/src/hg/lib/txAliDiff.as \\ ncbiRefSeqGenomicDiff.bed process/\$db.chrom.sizes ncbiRefSeqGenomicDiff.bb ln -s `pwd`/ncbiRefSeqGenomicDiff.bb $gbdbDir/ncbiRefSeqGenomicDiff.bb else rm -f ncbiRefSeqGenomicDiff.bed fi fi if [ -d "/usr/local/apache/htdocs-hgdownload/goldenPath/archive" ]; then - gtfFile=`ls \$db.*.ncbiRefSeq.gtf.gz` - mkdir -p /usr/local/apache/htdocs-hgdownload/goldenPath/archive/\$db/ncbiRefSeq - rm -f /usr/local/apache/htdocs-hgdownload/goldenPath/archive/\$db/ncbiRefSeq/\$gtfFile + mkdir -p /usr/local/apache/htdocs-hgdownload/goldenPath/archive/\$db/ncbiRefSeq/\$verString + rm -f /usr/local/apache/htdocs-hgdownload/goldenPath/archive/\$db/ncbiRefSeq/\$verString/\$db.\$verString.ncbiRefSeq.gtf.gz ln -s `pwd`/\$db.*.ncbiRefSeq.gtf.gz \\ - /usr/local/apache/htdocs-hgdownload/goldenPath/archive/\$db/ncbiRefSeq/ + /usr/local/apache/htdocs-hgdownload/goldenPath/archive/\$db/ncbiRefSeq/\$verString/\$db.\$verString.ncbiRefSeq.gtf.gz fi rm -f /usr/local/apache/htdocs-hgdownload/goldenPath/\$db/bigZips/genes/\$db.ncbiRefSeq.gtf.gz mkdir -p /usr/local/apache/htdocs-hgdownload/goldenPath/\$db/bigZips/genes ln -s `pwd`/\$db.*.ncbiRefSeq.gtf.gz \\ /usr/local/apache/htdocs-hgdownload/goldenPath/\$db/bigZips/genes/\$db.ncbiRefSeq.gtf.gz featureBits \$db ncbiRefSeq > fb.ncbiRefSeq.\$db.txt 2>&1 cat fb.ncbiRefSeq.\$db.txt 2>&1 _EOF_ ); } # if ($dbExists) $bossScript->execute(); } # doLoad