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

interim update 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@3173 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
depristo 2010-04-15 14:25:02 +00:00
parent b54031fc86
commit 7973806716
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
python/Verify1KGArchiveBAMs.py
View File

@ -200,6 +200,7 @@ def compareAlignmentIndices(remoteAlignmentIndex, alignmentIndex):
print ' md5s: local=%s remote=%s' % (raImd5, laImd5)
if raImd5 <> laImd5:
print ' [FAIL] -- alignment indices do not have the same hash!'
sys.exit(1)
else:
print ' [PASS] -- alignment indices are the same'
@ -250,6 +251,9 @@ if __name__ == "__main__":
FTPSERVER = FTP(ftpParsed[1])
FTPSERVER.login()
displayChangeLog(OPTIONS.remoteChangeLog)
compareAlignmentIndices(OPTIONS.remoteAlignmentIndex, OPTIONS.alignmentIndex)
results = dict()
for file in itertools.chain(readAlignmentIndex(OPTIONS.alignmentIndex), filesInLocalPath(root, OPTIONS.scanLocal )):
#print line
@ -258,9 +262,6 @@ if __name__ == "__main__":
compared = validateFile( file, root, ftpParsed[2] )
results[file] = compared
#localIndex
compareAlignmentIndices(OPTIONS.remoteAlignmentIndex, OPTIONS.alignmentIndex)
displayChangeLog(OPTIONS.remoteChangeLog)
printHeaderSep()
print 'SUMMARY: Total files examined', len(results)