6091558bab452de0ba020f3e08a9d20ef9aa3d81 hiram Wed Mar 18 19:09:58 2020 -0700 lastz runs for 12-way alignment on danRer10 per research request no redmine diff --git src/hg/makeDb/doc/danRer10/lastzRuns.txt src/hg/makeDb/doc/danRer10/lastzRuns.txt index 19863a2..484d18d 100644 --- src/hg/makeDb/doc/danRer10/lastzRuns.txt +++ src/hg/makeDb/doc/danRer10/lastzRuns.txt @@ -1,362 +1,584 @@ ######################################################################### # SWAP lastz zebrafish danRer10 (DONE - 2015-09-18 - Hiram) # the primary alignment on mm10 cd /hive/data/genomes/mm10/bed/lastzDanRer10.2015-09-11 cat fb.mm10.chainDanRer10Link.txt # 73464192 bases of 2652783500 (2.769%) in intersection # and for the swap mkdir /hive/data/genomes/danRer10/bed/blastz.mm10.swap cd /hive/data/genomes/danRer10/bed/blastz.mm10.swap time (doBlastzChainNet.pl -verbose=2 \ /hive/data/genomes/mm10/bed/lastzDanRer10.2015-09-11/DEF \ -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ -swap -chainMinScore=5000 -chainLinearGap=loose) > swap.log 2>&1 & # real 16m8.387s cat fb.danRer10.chainMm10Link.txt # 71611488 bases of 1369683683 (5.228%) in intersection time (doRecipBest.pl -buildDir=`pwd` danRer10 mm10) > rbest.log 2>&1 # real 7m34.259s ######################################################################### # LASTZ Xenopus xenTro7 (DONE - 2017-02-09 - Hiram) mkdir /hive/data/genomes/danRer10/bed/lastzXenTro7.2017-02-09 cd /hive/data/genomes/danRer10/bed/lastzXenTro7.2017-02-09 printf '# zebrafish vs. X. tropicalis BLASTZ=/cluster/bin/penn/lastz-distrib-1.03.66/bin/lastz BLASTZ_M=254 BLASTZ_H=2000 BLASTZ_Y=3400 BLASTZ_L=6000 BLASTZ_K=2200 BLASTZ_Q=/scratch/data/blastz/HoxD55.q # TARGET: Zebrafish danRer10 SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes SEQ1_CHUNK=20000000 SEQ1_LAP=10000 SEQ1_LIMIT=5 # QUERY: X. tropicalis xenTro7 SEQ2_DIR=/hive/data/genomes/xenTro7/xenTro7.2bit SEQ2_LEN=/hive/data/genomes/xenTro7/chrom.sizes SEQ2_CHUNK=10000000 SEQ2_LAP=0 SEQ2_LIMIT=40 BASE=/hive/data/genomes/danRer10/bed/lastzXenTro7.2017-02-09 TMPDIR=/dev/shm ' > DEF # establish a screen to control this job screen time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ -chainMinScore=5000 -chainLinearGap=loose \ -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ -bigClusterHub=ku) > do.log 2>&1 # lost tty connection time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ -chainMinScore=5000 -chainLinearGap=loose \ -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ -continue=cat -bigClusterHub=ku) > cat.log 2>&1 # real 136m6.704s cat fb.danRer10.chainXenTro7Link.txt # 90862965 bases of 1369683683 (6.634%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` danRer10 xenTro7) \ > rbest.log 2>&1 & # real 387m16.949s # and for the swap: mkdir /hive/data/genomes/xenTro7/bed/blastz.danRer10.swap cd /hive/data/genomes/xenTro7/bed/blastz.danRer10.swap time (doBlastzChainNet.pl -verbose=2 \ /hive/data/genomes/danRer10/bed/lastzXenTro7.2017-02-09/DEF \ -chainMinScore=5000 -chainLinearGap=loose \ -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ -syntenicNet -swap) > swap.log 2>&1 # real 30m3.142s cat fb.xenTro7.chainDanRer10Link.txt # 96433154 bases of 1365936747 (7.060%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` xenTro7 danRer10) \ > rbest.log 2>&1 # real 408m19.596s ######################################################################## # LASTZ green spotted puffer tetNig2 (DONE - 2017-02-09 - Hiram) screen -S danRer10TetNig2 mkdir /hive/data/genomes/danRer10/bed/lastzTetNig2.2017-02-09 cd /hive/data/genomes/danRer10/bed/lastzTetNig2.2017-02-09 printf '# zebrafish vs. green spotted puffer BLASTZ=/cluster/bin/penn/lastz-distrib-1.03.66/bin/lastz BLASTZ_M=254 BLASTZ_H=2000 BLASTZ_Y=9400 BLASTZ_L=3000 BLASTZ_K=3000 BLASTZ_Q=/scratch/data/blastz/HoxD55.q # TARGET: Zebrafish danRer10 SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes SEQ1_CHUNK=20000000 SEQ1_LAP=10000 SEQ1_LIMIT=1 # QUERY: green spotted puffer tetNig2 SEQ2_DIR=/hive/data/genomes/tetNig2/tetNig2.2bit SEQ2_LEN=/hive/data/genomes/tetNig2/chrom.sizes SEQ2_CHUNK=10000000 SEQ2_LAP=0 SEQ2_LIMIT=1 BASE=/hive/data/genomes/danRer10/bed/lastzTetNig2.2017-02-09 TMPDIR=/dev/shm ' > DEF # establish a screen to control this job screen time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ -bigClusterHub=ku) > do.log 2>&1 # real 97m22.745s cat fb.danRer10.chainTetNig2Link.txt # 83044644 bases of 1369683683 (6.063%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` danRer10 tetNig2) \ > rbest.log 2>&1 & # real 254m45.837s mkdir /hive/data/genomes/tetNig2/bed/blastz.danRer10.swap cd /hive/data/genomes/tetNig2/bed/blastz.danRer10.swap time (doBlastzChainNet.pl -verbose=2 \ /hive/data/genomes/danRer10/bed/lastzTetNig2.2017-02-09/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ -swap) > swap.log 2>&1 # real 6m59.881s cat fb.tetNig2.chainDanRer10Link.txt # 70573446 bases of 302314788 (23.344%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` tetNig2 danRer10) \ > rbest.log 2>&1 # real 228m19.188s time (doBlastzChainNet.pl -verbose=2 \ /hive/data/genomes/danRer10/bed/lastzTetNig2.2017-02-09/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ -continue=syntenicNet -syntenicNet -swap) > synNet.log 2>&1 # real 0m43.708s ######################################################################## # LASTZ fugu fr3 (DONE - 2017-02-09 - Hiram) screen -S danRer10Fr3 mkdir /hive/data/genomes/danRer10/bed/lastzFr3.2017-02-09 cd /hive/data/genomes/danRer10/bed/lastzFr3.2017-02-09 printf '# zebrafish vs. fugu BLASTZ=/cluster/bin/penn/lastz-distrib-1.03.66/bin/lastz BLASTZ_M=254 BLASTZ_H=2000 BLASTZ_Y=9400 BLASTZ_L=3000 BLASTZ_K=3000 BLASTZ_Q=/scratch/data/blastz/HoxD55.q # TARGET: Zebrafish danRer10 SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes SEQ1_CHUNK=20000000 SEQ1_LAP=10000 SEQ1_LIMIT=6 # QUERY: fugu fr3 SEQ2_DIR=/hive/data/genomes/fr3/fr3.2bit SEQ2_LEN=/hive/data/genomes/fr3/chrom.sizes SEQ2_CHUNK=10000000 SEQ2_LAP=0 SEQ2_LIMIT=30 BASE=/hive/data/genomes/danRer10/bed/lastzFr3.2017-02-09 TMPDIR=/dev/shm ' > DEF # establish a screen to control this job screen time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ -bigClusterHub=ku) > do.log 2>&1 # real 133m18.224s cat fb.danRer10.chainFr3Link.txt # 96978660 bases of 1369683683 (7.080%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` danRer10 fr3) \ > rbest.log 2>&1 & # real 200m15.132s mkdir /hive/data/genomes/fr3/bed/blastz.danRer10.swap cd /hive/data/genomes/fr3/bed/blastz.danRer10.swap time (doBlastzChainNet.pl -verbose=2 \ /hive/data/genomes/danRer10/bed/lastzFr3.2017-02-09/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ -syntenicNet -swap) > swap.log 2>&1 # real 9m31.601s cat fb.fr3.chainDanRer10Link.txt # 79663999 bases of 350961831 (22.699%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` fr3 danRer10) \ > rbest.log 2>&1 # real 181m55.550s ######################################################################## # LASTZ medaka oryLat2 (DONE - 2017-02-09 - Hiram) screen -S oryLat2DanRer10 mkdir /hive/data/genomes/danRer10/bed/lastzOryLat2.2017-02-09 cd /hive/data/genomes/danRer10/bed/lastzOryLat2.2017-02-09 printf '# zebrafish vs. medaka BLASTZ=/cluster/bin/penn/lastz-distrib-1.03.66/bin/lastz BLASTZ_M=254 BLASTZ_H=2000 BLASTZ_Y=9400 BLASTZ_L=3000 BLASTZ_K=3000 BLASTZ_Q=/scratch/data/blastz/HoxD55.q # TARGET: Zebrafish danRer10 SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes SEQ1_CHUNK=20000000 SEQ1_LAP=10000 SEQ1_LIMIT=6 # QUERY: medaka oryLat2 SEQ2_DIR=/hive/data/genomes/oryLat2/oryLat2.2bit SEQ2_LEN=/hive/data/genomes/oryLat2/chrom.sizes SEQ2_CHUNK=10000000 SEQ2_LAP=0 SEQ2_LIMIT=30 BASE=/hive/data/genomes/danRer10/bed/lastzOryLat2.2017-02-09 TMPDIR=/dev/shm ' > DEF # establish a screen to control this job screen time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ -bigClusterHub=ku) > do.log 2>&1 # real 202m52.634s cat fb.danRer10.chainOryLat2Link.txt # 132042172 bases of 1369683683 (9.640%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` danRer10 oryLat2) \ > rbest.log 2>&1 & # real 388m46.359s mkdir /hive/data/genomes/oryLat2/bed/blastz.danRer10.swap cd /hive/data/genomes/oryLat2/bed/blastz.danRer10.swap time (doBlastzChainNet.pl -verbose=2 \ /hive/data/genomes/danRer10/bed/lastzOryLat2.2017-02-09/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ -syntenicNet -swap) > swap.log 2>&1 # real 18m21.238s cat fb.oryLat2.chainDanRer10Link.txt # 123853592 bases of 700386597 (17.684%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` oryLat2 danRer10) \ > rbest.log 2>&1 # real 370m26.030s ######################################################################## # LASTZ stickleback gasAcu1 (DONE - 2017-02-09 - Hiram) screen -S gasAcu1DanRer10 mkdir /hive/data/genomes/danRer10/bed/lastzGasAcu1.2017-02-09 cd /hive/data/genomes/danRer10/bed/lastzGasAcu1.2017-02-09 printf '# zebrafish vs. stickleback BLASTZ=/cluster/bin/penn/lastz-distrib-1.03.66/bin/lastz BLASTZ_M=254 BLASTZ_H=2000 BLASTZ_Y=9400 BLASTZ_L=3000 BLASTZ_K=3000 BLASTZ_Q=/scratch/data/blastz/HoxD55.q # TARGET: Zebrafish danRer10 SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes SEQ1_CHUNK=20000000 SEQ1_LAP=10000 SEQ1_LIMIT=6 # QUERY: stickleback gasAcu1 # 30,000,000 is enough to work with the biggest chrom as one item SEQ2_DIR=/hive/data/genomes/gasAcu1/gasAcu1.2bit SEQ2_LEN=/hive/data/genomes/gasAcu1/chrom.sizes SEQ2_CHUNK=30000000 SEQ2_LAP=0 SEQ2_LIMIT=30 BASE=/hive/data/genomes/danRer10/bed/lastzGasAcu1.2017-02-09 TMPDIR=/dev/shm ' > DEF # establish a screen to control this job screen time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ -bigClusterHub=ku) > do.log 2>&1 # real 766m31.436s cat fb.danRer10.chainGasAcu1Link.txt # 134422269 bases of 1369683683 (9.814%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` danRer10 gasAcu1) \ > rbest.log 2>&1 & # real 466m55.498s mkdir /hive/data/genomes/gasAcu1/bed/blastz.danRer10.swap cd /hive/data/genomes/gasAcu1/bed/blastz.danRer10.swap time (doBlastzChainNet.pl -verbose=2 \ /hive/data/genomes/danRer10/bed/lastzGasAcu1.2017-02-09/DEF \ -chainMinScore=2000 -chainLinearGap=medium \ -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ -syntenicNet -swap) > swap.log 2>&1 # real 26m37.898s cat fb.gasAcu1.chainDanRer10Link.txt # 123076224 bases of 446627861 (27.557%) in intersection time (doRecipBest.pl -workhorse=hgwdev -buildDir=`pwd` gasAcu1 danRer10) \ > rbest.log 2>&1 # real 441m12.898s ######################################################################## +# LASTZ Coelacanth latCha1 (DONE - 2020-02-25 - Hiram) + mkdir /hive/data/genomes/danRer10/bed/lastzLatCha1.2020-02-25 + cd /hive/data/genomes/danRer10/bed/lastzLatCha1.2020-02-25 + + printf '# zebrafish vs. Coelacanth +BLASTZ=/cluster/bin/penn/lastz-distrib-1.04.00/bin/lastz +BLASTZ_M=254 + +BLASTZ_H=2000 +BLASTZ_Y=9400 +BLASTZ_L=3000 +BLASTZ_K=3000 +BLASTZ_Q=/cluster/data/blastz/HoxD55.q + +# TARGET: Zebrafish danRer10 +SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit +SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes +SEQ1_CHUNK=20000000 +SEQ1_LAP=10000 +SEQ1_LIMIT=6 + +# QUERY: Coelacanth latCha1 +# 14,000,000 is enough to work with the biggest chrom as one item +SEQ2_DIR=/hive/data/genomes/latCha1/latCha1.2bit +SEQ2_LEN=/hive/data/genomes/latCha1/chrom.sizes +SEQ2_CHUNK=14000000 +SEQ2_LAP=0 +SEQ2_LIMIT=30 + +BASE=/hive/data/genomes/danRer10/bed/lastzLatCha1.2020-02-25 +TMPDIR=/dev/shm +' > DEF + + # establish a screen to control this job + screen + time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ + -chainMinScore=5000 -chainLinearGap=loose \ + -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ + -bigClusterHub=ku) > do.log 2>&1 + # real 177m11.851s + + cat fb.danRer10.chainLatCha1Link.txt + # 86552427 bases of 1369683683 (6.319%) in intersection + cat fb.danRer10.chainSynLatCha1Link.txt + # 15892801 bases of 1369683683 (1.160%) in intersection + + time (doRecipBest.pl -load -workhorse=hgwdev -buildDir=`pwd` \ + danRer10 latCha1) > rbest.log 2>&1 & + # real 71m59.526s + + cat fb.danRer10.chainRBest.LatCha1.txt + # 56607816 bases of 1369683683 (4.133%) in intersection + + mkdir /hive/data/genomes/latCha1/bed/blastz.danRer10.swap + cd /hive/data/genomes/latCha1/bed/blastz.danRer10.swap + time (doBlastzChainNet.pl -verbose=2 \ + /hive/data/genomes/danRer10/bed/lastzLatCha1.2020-02-25/DEF \ + -chainMinScore=5000 -chainLinearGap=loose \ + -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ + -syntenicNet -swap) > swap.log 2>&1 + # real 12m58.361s + + cat fb.latCha1.chainDanRer10Link.txt + # 98338936 bases of 2183592768 (4.504%) in intersection + cat fb.latCha1.chainSynDanRer10Link.txt + # 13550035 bases of 2183592768 (0.621%) in intersection + + time (doRecipBest.pl -load -workhorse=hgwdev -buildDir=`pwd` \ + latCha1 danRer10) > rbest.log 2>&1 + # real 67m29.683s + cat fb.latCha1.chainRBest.DanRer10.txt + # 56742837 bases of 2183592768 (2.599%) in intersection + +######################################################################## +# LASTZ spotted gar lepOcu1 (DONE - 2020-02-25 - Hiram) + mkdir /hive/data/genomes/danRer10/bed/lastzLepOcu1.2020-02-25 + cd /hive/data/genomes/danRer10/bed/lastzLepOcu1.2020-02-25 + + printf '# zebrafish vs. spotted gar +BLASTZ=/cluster/bin/penn/lastz-distrib-1.04.00/bin/lastz +BLASTZ_M=254 + +BLASTZ_H=2000 +BLASTZ_Y=9400 +BLASTZ_L=3000 +BLASTZ_K=3000 +BLASTZ_Q=/cluster/data/blastz/HoxD55.q + +# TARGET: Zebrafish danRer10 +SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit +SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes +SEQ1_CHUNK=20000000 +SEQ1_LAP=10000 +SEQ1_LIMIT=6 + +# QUERY: spotted gar lepOcu1 +SEQ2_DIR=/hive/data/genomes/lepOcu1/lepOcu1.2bit +SEQ2_LEN=/hive/data/genomes/lepOcu1/chrom.sizes +SEQ2_CHUNK=20000000 +SEQ2_LAP=0 +SEQ2_LIMIT=30 + +BASE=/hive/data/genomes/danRer10/bed/lastzLepOcu1.2020-02-25 +TMPDIR=/dev/shm +' > DEF + + # establish a screen to control this job + screen + time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ + -chainMinScore=5000 -chainLinearGap=loose \ + -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ + -bigClusterHub=ku) > do.log 2>&1 + # real 118m59.600s + + cat fb.danRer10.chainLepOcu1Link.txt + # 115180491 bases of 1369683683 (8.409%) in intersection + cat fb.danRer10.chainSynLepOcu1Link.txt + # 58800217 bases of 1369683683 (4.293%) in intersection + + time (doRecipBest.pl -load -workhorse=hgwdev -buildDir=`pwd` danRer10 \ + lepOcu1) > rbest.log 2>&1 & + # real 88m25.172s + + cat fb.danRer10.chainRBest.LepOcu1.txt + # 78985057 bases of 1369683683 (5.767%) in intersection + + mkdir /hive/data/genomes/lepOcu1/bed/blastz.danRer10.swap + cd /hive/data/genomes/lepOcu1/bed/blastz.danRer10.swap + time (doBlastzChainNet.pl -verbose=2 \ + /hive/data/genomes/danRer10/bed/lastzLepOcu1.2020-02-25/DEF \ + -chainMinScore=5000 -chainLinearGap=loose \ + -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ + -syntenicNet -swap) > swap.log 2>&1 + # real 7m34.876s + + cat fb.lepOcu1.chainDanRer10Link.txt + # 98189747 bases of 869414361 (11.294%) in intersection + cat fb.lepOcu1.chainSynDanRer10Link.txt + # 47688590 bases of 869414361 (5.485%) in intersection + + time (doRecipBest.pl -load -workhorse=hgwdev -buildDir=`pwd` \ + lepOcu1 danRer10) > rbest.log 2>&1 + # real 83m11.208s + + cat fb.lepOcu1.chainRBest.DanRer10.txt + # 78000762 bases of 869414361 (8.972%) in intersection + +######################################################################## +# LASTZ Elephant shark calMil1 (DONE - 2020-02-25 - Hiram) + mkdir /hive/data/genomes/danRer10/bed/lastzCalMil1.2020-02-25 + cd /hive/data/genomes/danRer10/bed/lastzCalMil1.2020-02-25 + + printf '# zebrafish vs. Elephant shark +BLASTZ=/cluster/bin/penn/lastz-distrib-1.04.00/bin/lastz +BLASTZ_M=254 + +BLASTZ_H=2000 +BLASTZ_Y=9400 +BLASTZ_L=3000 +BLASTZ_K=3000 +BLASTZ_Q=/cluster/data/blastz/HoxD55.q + +# TARGET: Zebrafish danRer10 +SEQ1_DIR=/hive/data/genomes/danRer10/danRer10.2bit +SEQ1_LEN=/hive/data/genomes/danRer10/chrom.sizes +SEQ1_CHUNK=20000000 +SEQ1_LAP=10000 +SEQ1_LIMIT=6 + +# QUERY: Elephant shark calMil1 +SEQ2_DIR=/hive/data/genomes/calMil1/calMil1.2bit +SEQ2_LEN=/hive/data/genomes/calMil1/chrom.sizes +SEQ2_CHUNK=20000000 +SEQ2_LAP=0 +SEQ2_LIMIT=100 + +BASE=/hive/data/genomes/danRer10/bed/lastzCalMil1.2020-02-25 +TMPDIR=/dev/shm +' > DEF + + # establish a screen to control this job + screen + time (doBlastzChainNet.pl -verbose=2 `pwd`/DEF \ + -chainMinScore=5000 -chainLinearGap=loose \ + -syntenicNet -workhorse=hgwdev -smallClusterHub=ku \ + -bigClusterHub=ku) > do.log 2>&1 + # real 48m39.821s + + cat fb.danRer10.chainCalMil1Link.txt + # 76690714 bases of 1369683683 (5.599%) in intersection + cat fb.danRer10.chainSynCalMil1Link.txt + # 18240388 bases of 1369683683 (1.332%) in intersection + + time (doRecipBest.pl -load -workhorse=hgwdev -buildDir=`pwd` danRer10 \ + calMil1) > rbest.log 2>&1 & + # real 466m55.498s + + cat fb.danRer10.chainRBest.CalMil1.txt + # 48435763 bases of 1369683683 (3.536%) in intersection + + mkdir /hive/data/genomes/calMil1/bed/blastz.danRer10.swap + cd /hive/data/genomes/calMil1/bed/blastz.danRer10.swap + time (doBlastzChainNet.pl -verbose=2 \ + /hive/data/genomes/danRer10/bed/lastzCalMil1.2020-02-25/DEF \ + -chainMinScore=5000 -chainLinearGap=loose \ + -workhorse=hgwdev -smallClusterHub=ku -bigClusterHub=ku \ + -syntenicNet -swap) > swap.log 2>&1 + # real 7m34.876s + + cat fb.calMil1.chainDanRer10Link.txt + # 66411165 bases of 936953458 (7.088%) in intersection + cat fb.calMil1.chainSynDanRer10Link.txt + # 15359860 bases of 936953458 (1.639%) in intersection + + time (doRecipBest.pl -load -workhorse=hgwdev -buildDir=`pwd` \ + calMil1 danRer10) > rbest.log 2>&1 + # real 57m38.923s + + cat fb.calMil1.chainRBest.DanRer10.txt + # 47381650 bases of 936953458 (5.057%) in intersection + +########################################################################