Mercurial > repos > davidvanzessen > mutation_analysis
diff gene_identification.py @ 64:0fdd90f7c654 draft
Uploaded
author | davidvanzessen |
---|---|
date | Fri, 01 Apr 2016 08:54:24 -0400 |
parents | 4262e880472d |
children | b869a126e2c4 |
line wrap: on
line diff
--- a/gene_identification.py Fri Mar 25 07:50:12 2016 -0400 +++ b/gene_identification.py Fri Apr 01 08:54:24 2016 -0400 @@ -25,13 +25,14 @@ with open(infile, 'r') as f: #read all sequences into a dictionary as key = ID, value = sequence for line in f: total += 1 + linesplt = line.split("\t") if first: - linesplt = line.split("\t") + print "linesplt", linesplt IDIndex = linesplt.index("Sequence ID") seqIndex = linesplt.index("Sequence") first = False continue - linesplt = line.split("\t") + ID = linesplt[IDIndex] if len(linesplt) < 28: #weird rows without a sequence dic[ID] = ""