zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Better support for reading truth file 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@2307 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
depristo 2009-12-10 12:16:05 +00:00
parent 0fae798b3a
commit b2dfe85648
1 changed files with 0 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
python/snpSelector.py
View File

@ -621,12 +621,6 @@ def main():
truthVCF = OPTIONS.truth
TRUTH_CALLS = readTruth(truthVCF)
if len(args) > 1:
truthVCF = args[1]
TRUTH_CALLS = readTruth(truthVCF)
#for v in TRUTH_CALLS.itervalues(): print v.getField("NA12878")
#sys.exit(1)
if OPTIONS.recalLog <> None:
RECAL_LOG = open(OPTIONS.recalLog, "w")
print >> RECAL_LOG, "# optimized vcf", args[0]