src/hg/blastToPsl/tests/makefile 1.23
1.23 2010/02/08 03:04:24 markd
added blastXmlToPsl
Index: src/hg/blastToPsl/tests/makefile
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/blastToPsl/tests/makefile,v
retrieving revision 1.22
retrieving revision 1.23
diff -b -B -U 1000000 -r1.22 -r1.23
--- src/hg/blastToPsl/tests/makefile 6 Feb 2010 08:17:50 -0000 1.22
+++ src/hg/blastToPsl/tests/makefile 8 Feb 2010 03:04:24 -0000 1.23
@@ -1,149 +1,167 @@
-include ../../../inc/common.mk
+include ../../../../inc/common.mk
#blastToPsl = valgrind --tool=memcheck ../blastToPsl
-blastToPsl = ../blastToPsl
+blastToPsl = blastToPsl
+
+#blastXmlToPsl = valgrind --tool=memcheck ../blastToPsl
+blastXmlToPsl = blastXmlToPsl
all:
-test: protMrnaTest protDnaTest mrnaDnaTest mrnaDnaxTest \
+test: testBlastToPsl testBlastXmlToPsl
+
+###
+# blastToPsl
+###
+testBlastToPsl: protMrnaTest protDnaTest mrnaDnaTest mrnaDnaxTest \
blockBugTest identityBugTest protMrnaEval1Test protMrnaEval2Test \
longDbNameTest multiLineDescTest tblastxTest tblastnTest \
- blastpPsiTest psiNoHitsTest psiBadFormatTest blastpPsiWrapTest \
+ blastpPsiTest psiNoHitsTest psiBadFormatTest blastpPsiWrapTest
+
+###
+# blastXmlToPsl
+###
+testBlastXmlToPsl: mrna1Test prot1Test prot2Test
-# blast+ not implemented yet
-# blastn+Test
+mrna1Test: mkdirs
+ ${blastXmlToPsl} -scores=output/$@.bscores input/mrna1.blastn.xml output/$@.psl
+ pslCheck -verbose=0 output/$@.psl
+ diff -u expected/$@.psl output/$@.psl
+ diff -u expected/$@.bscores output/$@.bscores
+
+prot1Test: mkdirs
+ ${blastXmlToPsl} -scores=output/$@.bscores input/prot1.tblastn.xml output/$@.psl
+ pslCheck -verbose=0 output/$@.psl
+ diff -u expected/$@.psl output/$@.psl
+ diff -u expected/$@.bscores output/$@.bscores
-#tblastnDbFirstTest megaBlastTest
+prot2Test: mkdirs
+ ${blastXmlToPsl} -scores=output/$@.bscores input/prot2.psiblast.xml output/$@.psl
+ pslCheck -verbose=0 output/$@.psl
+ diff -u expected/$@.psl output/$@.psl
+ diff -u expected/$@.bscores output/$@.bscores
protMrnaTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/kgProtMRna1.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
protDnaTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/kgProtDna1.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
mrnaDnaTest: mkdirs
${blastToPsl} input/mrnaDna1.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
mrnaDnaxTest: mkdirs
${blastToPsl} -pslx input/mrnaDna1.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
# test case that cause huge target block start.
blockBugTest: mkdirs
${blastToPsl} input/blockBug.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
# handle id
identityBugTest: mkdirs
${blastToPsl} input/identityBug.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
# test E-value threshold filter
protMrnaEval1Test: mkdirs
${blastToPsl} -eVal=1 -scores=output/$@.bscores input/kgProtMRna1.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
protMrnaEval2Test: mkdirs
${blastToPsl} -eVal=1e-10 -scores=output/$@.bscores input/kgProtMRna1.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
longDbNameTest: mkdirs
${blastToPsl} input/longDbName.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
# multiline description broke parser
multiLineDescTest: mkdirs
${blastToPsl} input/multiLineDesc.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
# tblastx
tblastxTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/tblastx.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
# tblastn
tblastnTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/tblastn.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
# this test is currently disabled. This file has Database: before Query= and
# there are multiple queries for each Database:
tblastnDbFirstTest: mkdirs
${blastToPsl} input/tblastnDbFirst.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
# this test is currently disabled, mega-blast is painful because
# Database: comes before Query= and there are multiple queries for each
# Database:.
megaBlastTest: mkdirs
${blastToPsl} input/970.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
# PSI-BLASTP
blastpPsiTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/blastpPsi.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
# PSI-BLASTP with Identities line wrapping
blastpPsiWrapTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/blastpPsiWrap.blast output/$@.psl
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
# PSI-BLASTP with no hits
psiNoHitsTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/psiNoHits.blast output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
# PSI-BLASTP produces bogus blocks where query starts at 0 and there is not
# ending number on the line. Sbjct may be empty or have sequence. These
# are ignored
psiBadFormatTest: mkdirs
${blastToPsl} -scores=output/$@.bscores input/psiBadFormat.blast output/$@.psl >output/$@.err 2>&1
pslCheck -verbose=0 output/$@.psl
diff -u expected/$@.psl output/$@.psl
diff -u expected/$@.bscores output/$@.bscores
diff -u expected/$@.err output/$@.err
-# blast+ changes output format
-blastn+Test: mkdirs
- ${blastToPsl} -scores=output/$@.bscores input/blastn+.blast output/$@.psl >output/$@.err 2>&1
- pslCheck -verbose=0 output/$@.psl
- diff -u expected/$@.psl output/$@.psl
- diff -u expected/$@.bscores output/$@.bscores
-
-
clean::
rm -rf output
mkdirs:
@${MKDIR} output