def parse_iedb_rec(file, cutoffLOWER, cutoffUPPER):
fileHandle = open(file)
patternHash = {}
lineCount = 0
HLA_Set = set([])
for line in fileHandle:
if (lineCount < 4):
lineCount += 1
continue
record = line.rstrip().split('\t')
if (len(record) > 1):
if (float(record[7]) > cutoffLOWER):
continue
if (float(record[7]) < cutoffUPPER):
continue
if (record[5] in patternHash):
patternHash[record[5]].append(
CF.Epitope_Info(record[0], record[7], record[2],
record[3]))
else:
patternHash[record[5]] = [
CF.Epitope_Info(record[0], record[7], record[2], record[3])
]
if (record[0] not in HLA_Set):
HLA_Set.add(record[0])
return (patternHash, HLA_Set)
