Mercurial > repos > davidvanzessen > mutation_analysis

comparison gene_identification.py @ 60:1558fbbc19ef draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Fri, 18 Mar 2016 08:02:50 -0400
parents c9f9623f1f76
children 64e6a7803e07
comparison
equal deleted inserted replaced
59:79696e273853 60:1558fbbc19ef
38 else: 38 else:
39 dic[ID] = linesplt[seqIndex] 39 dic[ID] = linesplt[seqIndex]
40 40
41 print "Number of input sequences:", len(dic) 41 print "Number of input sequences:", len(dic)
42 42
43 #old cm sequence: gggagtgcatccgccccaacccttttccccctcgtctcctgtgagaattccc
44
43 #lambda/kappa reference sequence 45 #lambda/kappa reference sequence
44 searchstrings = {"ca": "catccccgaccagccccaaggtcttcccgctgagcctctgcagcacccagccagatgggaacgtggtcatcgcctgcctgg", 46 searchstrings = {"ca": "catccccgaccagccccaaggtcttcccgctgagcctctgcagcacccagccagatgggaacgtggtcatcgcctgcctgg",
45 "cg": "ctccaccaagggcccatcggtcttccccctggcaccctcctccaagagcacctctgggggcacagcggccctgggctgcctggtcaaggactacttccccgaaccggtgacggtgtcgtggaactcaggcgccctgaccag", 47 "cg": "ctccaccaagggcccatcggtcttccccctggcaccctcctccaagagcacctctgggggcacagcggccctgggctgcctggtcaaggactacttccccgaaccggtgacggtgtcgtggaactcaggcgccctgaccag",
46 "cm": "gggagtgcatccgccccaacccttttccccctcgtctcctgtgagaattccc"} 48 "cm": "gggagtgcatccgccccaacc"} #new (shorter) cm sequence
47 49
48 compiledregex = {"ca": [], 50 compiledregex = {"ca": [],
49 "cg": [], 51 "cg": [],
50 "cm": []} 52 "cm": []}
51 53