comparison mutation_analysis.py @ 76:becea91089ed draft

Uploaded
author davidvanzessen
date Mon, 09 May 2016 03:45:39 -0400
parents 13c3710604ef
children b523ce95d857
comparison
equal deleted inserted replaced
75:14749ced7ff2 76:becea91089ed
55 continue 55 continue
56 linecount += 1 56 linecount += 1
57 linesplt = line.split("\t") 57 linesplt = line.split("\t")
58 ID = linesplt[IDIndex] 58 ID = linesplt[IDIndex]
59 genedic[ID] = linesplt[best_matchIndex] 59 genedic[ID] = linesplt[best_matchIndex]
60 print line
61 try: 60 try:
62 mutationdic[ID + "_FR1"] = [mutationMatcher.match(x).groups() for x in linesplt[fr1Index].split("|") if x] if include_fr1 else [] 61 mutationdic[ID + "_FR1"] = [mutationMatcher.match(x).groups() for x in linesplt[fr1Index].split("|") if x] if include_fr1 else []
63 mutationdic[ID + "_CDR1"] = [mutationMatcher.match(x).groups() for x in linesplt[cdr1Index].split("|") if x] 62 mutationdic[ID + "_CDR1"] = [mutationMatcher.match(x).groups() for x in linesplt[cdr1Index].split("|") if x]
64 mutationdic[ID + "_FR2"] = [mutationMatcher.match(x).groups() for x in linesplt[fr2Index].split("|") if x] 63 mutationdic[ID + "_FR2"] = [mutationMatcher.match(x).groups() for x in linesplt[fr2Index].split("|") if x]
65 mutationdic[ID + "_CDR2"] = [mutationMatcher.match(x).groups() for x in linesplt[cdr2Index].split("|") if x] 64 mutationdic[ID + "_CDR2"] = [mutationMatcher.match(x).groups() for x in linesplt[cdr2Index].split("|") if x]