Mercurial > repos > davidvanzessen > shm_csr
diff shm_csr.py @ 40:ca2512e1e3ab draft
Uploaded
author | davidvanzessen |
---|---|
date | Thu, 29 Dec 2016 07:05:45 -0500 |
parents | ce25fb581ca3 |
children | 77a7ac76c7b9 |
line wrap: on
line diff
--- a/shm_csr.py Thu Dec 22 09:39:27 2016 -0500 +++ b/shm_csr.py Thu Dec 29 07:05:45 2016 -0500 @@ -38,7 +38,7 @@ cdr1LengthDic = {} cdr2LengthDic = {} -with open(infile, 'r') as i: +with open(infile, 'ru') as i: for line in i: if first: linesplt = line.split("\t") @@ -80,6 +80,8 @@ IDlist += [ID] +print mutationList, linecount + AALength = (int(max(mutationList, key=lambda i: int(i[4]) if i[4] else 0)[4]) + 1) # [4] is the position of the AA mutation, None if silent if AALength < 60: AALength = 64 @@ -173,7 +175,7 @@ aggctatIndex = 0 atagcctIndex = 0 first = True -with open(infile, 'r') as i: +with open(infile, 'ru') as i: for line in i: if first: linesplt = line.split("\t")