All.joiner problem

From genomewiki
Revision as of 03:03, 17 July 2010 by Rhead (talk | contribs) (all.joiner problem to be aware of)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigationJump to search

There is a problem with all.joiner rules and the -database parameter. This is taken from an email from Ann in 2007:

Bob and I just ran into something sneaky with running joinerCheck or runJoiner.csh. There are some rules that can be tricky and we are not exactly sure why. What we do know is the kind of output you can be on the lookout for to know that you're dealing with a "tricky" rule. The output sort of makes it look like the rule has run (with no problems), but it really hasn't run at all.

Example of rules that aren't really running from joinerCheck:

[hgwdev:~/joiner> joinerCheck -keys -database=fr2 -identifier=blastHg18KGPep04Id all.joiner
Checking keys on database fr2

from runJoiner.csh:

[hgwdev:~/joiner> runJoiner.csh fr2 blastHg18KG
found identifiers:

blastHg18KGRef04Id
blastHg18KGPep04Id

Checking keys on database fr2
Checking keys on database fr2


If the rule had really run there would be some output saying something like X hits out of Y ok.

The way to force these rules to run is to leave off the "database" parameter (you can only do this when using the joinerCheck (runJoiner.csh can't handle that)). This will run the rule for all databases that are listed -- just pick yours out of the output.

[hgwdev:~/joiner> joinerCheck -keys -identifier=blastHg18KGPep04Id all.joiner Checking keys on database (null)
 anoCar1.blastHg18KG.qName - hits 45332 of 45332 ok
 equCab1.blastHg18KG.qName - hits 61628 of 61628 ok
 ornAna1.blastHg18KG.qName - hits 35761 of 35761 ok
 oryLat1.blastHg18KG.qName - hits 38111 of 38111 ok
 bosTau2.blastHg18KG.qName - hits 56832 of 56832 ok
 strPur1.blastHg18KG.qName - hits 30907 of 30907 ok
 mm8.blastHg18KG.qName - hits 70754 of 70754 ok
 fr2.blastHg18KG.qName - hits 40510 of 40510 ok 
 dm2.blastHg18KG.qName - hits 22239 of 22239 ok
 felCat3.blastHg18KG.qName - hits 37715 of 37715 ok
 canFam2.blastHg18KG.qName - hits 62741 of 62741 ok
 danRer4.blastHg18KG.qName - hits 44692 of 44692 ok
 xenTro2.blastHg18KG.qName - hits 40291 of 40291 ok
 galGal3.blastHg18KG.qName - hits 35851 of 35851 ok
 monDom4.blastHg18KG.qName - hits 63264 of 63264 ok
 hg18.blastKGRef04.acc - hits 36727 of 36727 ok
 hg18.knownGene.name - hits 0 of 56722
Error: 56722 of 56722 elements of hg18.knownGene.name are not in key blastKGPep04.name line 1995 of all.joiner
Example miss: UC001AAA.1