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

Minor updates to correctly handle emitting FN calls 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@2231 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
depristo 2009-12-02 22:53:17 +00:00
parent c776f9fb90
commit dbb8b86ed1
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
python/snpSelector.py
View File

@ -571,10 +571,13 @@ def evaluateTruth(header, callVCF, truth, truthVCF):
writeRecalibratedCalls(callVCF, header, calls) writeRecalibratedCalls(callVCF, header, calls)
def isFN(v):
return isVariantInSample(v, OPTIONS.useSample) and not v.hasField("FN")
if truth <> None and OPTIONS.FNoutputVCF: if truth <> None and OPTIONS.FNoutputVCF:
f = open(OPTIONS.FNoutputVCF, 'w') f = open(OPTIONS.FNoutputVCF, 'w')
#print 'HEADER', header #print 'HEADER', header
for line in formatVCF(header, filter( lambda x: not x.hasField("FN"), truth.itervalues())): for line in formatVCF(header, filter( isFN, truth.itervalues())):
print >> f, line print >> f, line
f.close() f.close()