simple monitor for watching pilot 1 call progress
git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@1769 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
parent
1bd0c3c145
commit
c1e1d910cb
|
|
@ -0,0 +1,37 @@
|
|||
import farm_commands
|
||||
import os.path
|
||||
import sys
|
||||
from optparse import OptionParser
|
||||
from datetime import date
|
||||
import glob
|
||||
import operator
|
||||
|
||||
if __name__ == "__main__":
|
||||
usage = "usage: %prog files"
|
||||
parser = OptionParser(usage=usage)
|
||||
parser.add_option("-s", "--size", dest="regionSize",
|
||||
type="int", default=None, help="")
|
||||
|
||||
(OPTIONS, args) = parser.parse_args()
|
||||
|
||||
def extract(line):
|
||||
s = line.split()
|
||||
t = s[0].split(":")
|
||||
if len(t) == 2:
|
||||
chr, pos = t
|
||||
return chr, int(pos)
|
||||
else:
|
||||
return None, None
|
||||
|
||||
for file in args:
|
||||
chr, lastPos, startPos, calledBases, progress = None, None, None, None, None
|
||||
lastLine = None
|
||||
for line in open(file):
|
||||
if startPos == None:
|
||||
chr, startPos = extract(line)
|
||||
lastLine = line
|
||||
if lastLine <> None and startPos <> None:
|
||||
lastPos = extract(lastLine)[1]
|
||||
calledBases = lastPos - startPos
|
||||
progress = "%.2f" % (float(calledBases) / OPTIONS.regionSize * 100)
|
||||
print file, chr, startPos, lastPos, calledBases, progress
|
||||
Loading…
Reference in New Issue