src/hg/makeDb/doc/ailMel1.txt 1.10

1.10 2010/04/16 17:25:05 hiram
add swap mention for human, dog, rat
Index: src/hg/makeDb/doc/ailMel1.txt
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/doc/ailMel1.txt,v
retrieving revision 1.9
retrieving revision 1.10
diff -b -B -U 4 -r1.9 -r1.10
--- src/hg/makeDb/doc/ailMel1.txt	12 Apr 2010 19:27:04 -0000	1.9
+++ src/hg/makeDb/doc/ailMel1.txt	16 Apr 2010 17:25:05 -0000	1.10
@@ -191,9 +191,9 @@
         etc/hgwdev.dbs
     cvs ci -m "Adding ailMel1 - Panda - Ailuropoda melanoleuca" \
 	etc/align.dbs etc/hgwdev.dbs
     make etc-update
-    #	done - 2010-01-19 - Hiram
+    #	done - 2010-04-12 - Hiram
 
 ############################################################################
 # reset position to RHO location as found from blat of hg19 RHO gene
     hgsql -e \
@@ -357,4 +357,61 @@
     #	real    84m58.382s
     #	no new gaps
 
 #########################################################################
+# lastz Human hg19 swap (DONE - 2010-02-05 - Hiram)
+    # original alignment
+    cd /hive/data/genomes/hg19/bed/lastzAilMel1.2010-02-04
+    cat fb.hg19.chainAilMel1Link.txt 
+    #	1453400264 bases of 2897316137 (50.164%) in intersection
+
+    #	and the swap
+    mkdir /hive/data/genomes/ailMel1/bed/blastz.hg19.swap
+    cd /hive/data/genomes/ailMel1/bed/blastz.hg19.swap
+    time nice -n +19 doBlastzChainNet.pl -verbose=2 \
+	/hive/data/genomes/hg19/bed/lastzAilMel1.2010-02-04/DEF \
+	-swap -noLoadChainSplit -syntenicNet \
+	-workhorse=hgwdev -smallClusterHub=memk -bigClusterHub=swarm \
+	-chainMinScore=3000 -chainLinearGap=medium > swap.log 2>&1 &
+    real    124m14.393s
+    cat fb.ailMel1.chainHg19Link.txt 
+    #	1411953704 bases of 2245312831 (62.884%) in intersection
+
+#########################################################################
+# lastz Dog canFam2 swap (DONE - 2010-02-05 - Hiram)
+    # original alignment
+    cd /hive/data/genomes/canFam2/bed/lastzAilMel1.2010-02-04
+    cat fb.canFam2.chainAilMel1Link.txt 
+    #	1791212709 bases of 2384996543 (75.103%) in intersection
+
+    #	and the swap
+    mkdir /hive/data/genomes/ailMel1/bed/blastz.canFam2.swap
+    cd /hive/data/genomes/ailMel1/bed/blastz.canFam2.swap
+    time doBlastzChainNet.pl -verbose=2 \
+	/hive/data/genomes/canFam2/bed/lastzAilMel1.2010-02-04/DEF \
+	-swap -noLoadChainSplit -workhorse=hgwdev -bigClusterHub=pk \
+	-smallClusterHub=memk -chainMinScore=3000 -chainLinearGap=medium \
+	> swap.log 2>&1 &
+    #	real    128m41.005s
+     cat fb.ailMel1.chainCanFam2Link.txt 
+    #	1788107935 bases of 2245312831 (79.637%) in intersection
+
+#########################################################################
+# lastz Rat rn4 swap (DONE - 2010-02-18 - Hiram)
+    # original alignment
+    cd /hive/data/genomes/rn4/bed/lastzAilMel1.2010-02-04
+    cat fb.rn4.chainAilMel1Link.txt 
+    #	708197812 bases of 2571531505 (27.540%) in intersection
+
+    #	and the swap
+    mkdir /hive/data/genomes/ailMel1/bed/blastz.rn4.swap
+    cd /hive/data/genomes/ailMel1/bed/blastz.rn4.swap
+    time nice -n +19 doBlastzChainNet.pl -verbose=2 \
+	/hive/data/genomes/rn4/bed/lastzAilMel1.2010-02-04/DEF \
+	-swap -syntenicNet \
+	-workhorse=hgwdev -smallClusterHub=memk -bigClusterHub=swarm \
+	-chainMinScore=3000 -chainLinearGap=medium > swap.log 2>&1 &
+    #	real    73m37.751s
+    cat fb.ailMel1.chainRn4Link.txt 
+    #	695366144 bases of 2245312831 (30.970%) in intersection
+
+#########################################################################