13bfd800b82e7439df8f5dce583676890139d3c4 hiram Mon Mar 28 09:56:25 2022 -0700 create bpt during addMask refs #28961 diff --git src/hg/utils/automation/doAssemblyHub.pl src/hg/utils/automation/doAssemblyHub.pl index 18e47a6..408562d 100755 --- src/hg/utils/automation/doAssemblyHub.pl +++ src/hg/utils/automation/doAssemblyHub.pl @@ -1452,31 +1452,34 @@ my $src2BitToMask = "../repeatMasker/$asmId.rmsk.2bit"; if ($noRmsk || ($wmMasked > $rmMasked)) { $src2BitToMask = "../windowMasker/$asmId.cleanWMSdust.2bit"; } $bossScript->add(<<_EOF_ export asmId=$asmId export accessionId=`echo \$asmId | cut -d'_' -f1-2` if [ ../simpleRepeat/trfMask.bed.gz -nt \$asmId.masked.faSize.txt ]; then twoBitMask $src2BitToMask -type=.bed \\ -add ../simpleRepeat/trfMask.bed.gz \$asmId.masked.2bit twoBitToFa \$asmId.masked.2bit stdout | faSize stdin > \$asmId.masked.faSize.txt touch -r \$asmId.masked.2bit \$asmId.masked.faSize.txt + bptForTwoBit \$asmId.masked.2bit \$asmId.masked.2bit.bpt + touch -r \$asmId.masked.2bit \$asmId.masked.2bit.bpt cp -p \$asmId.masked.faSize.txt ../../\$asmId.faSize.txt + cp -p \$asmId.masked.2bit.bpt ../../\$asmId.2bit.bpt size=`grep -w bases \$asmId.masked.faSize.txt | cut -d' ' -f1` if [ \$size -lt 4294967297 ]; then ln \$asmId.masked.2bit \$accessionId.2bit gfServer -trans index ../../\$accessionId.trans.gfidx \$accessionId.2bit & gfServer -stepSize=5 index ../../\$accessionId.untrans.gfidx \$accessionId.2bit wait rm \$accessionId.2bit touch -r \$asmId.masked.2bit ../../\$accessionId.trans.gfidx touch -r \$asmId.masked.2bit ../../\$accessionId.untrans.gfidx else printf "# genome \$asmId too large at \$size to make blat indexes\\n" 1>&2 fi else printf "# addMask step previously completed\\n" 1>&2 exit 0