src/utils/qa/checkMetaData.csh 1.17
1.17 2009/04/13 16:13:28 kuhn
fixed bug where genomeClade comparison breaks because original tables not sorted
Index: src/utils/qa/checkMetaData.csh
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/utils/qa/checkMetaData.csh,v
retrieving revision 1.16
retrieving revision 1.17
diff -b -B -U 4 -r1.16 -r1.17
--- src/utils/qa/checkMetaData.csh 11 Apr 2009 00:37:01 -0000 1.16
+++ src/utils/qa/checkMetaData.csh 13 Apr 2009 16:13:28 -0000 1.17
@@ -202,15 +202,16 @@
# get lookup for clade check
# filter out "/" when it appears in genome name - to avoid e.g, Dog/Human
hgsql $host1 -Ne 'SELECT * FROM genomeClade WHERE genome LIKE "%'$genome'"' \
- $centdb1 | grep -v "/" > $metatable.$db.$out1
+ $centdb1 | grep -v "/" | sort > $metatable.$db.$out1
hgsql $host2 -Ne 'SELECT * FROM genomeClade WHERE genome LIKE "%'$genome'"' \
- $centdb2 | grep -v "/" > $metatable.$db.$out2
+ $centdb2 | grep -v "/" | sort > $metatable.$db.$out2
set metatable=""
# compare and print results
+# should replace with commTiro.csh, but for now simply sorting genomeClade above
foreach table ( `echo $metatables` )
comm -23 $table.$db.$out1 $table.$db.$out2 > $table.$db.${out1}Only
comm -13 $table.$db.$out1 $table.$db.$out2 > $table.$db.${out2}Only
comm -12 $table.$db.$out1 $table.$db.$out2 > $table.$db.common