-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcount2.py
53 lines (45 loc) · 1.22 KB
/
count2.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
import re
chr1 = 'BABa'
chr2 = 'ABaBABaB'
chr3 = 'BABABaBaaBabABaBA'
chr4 = 'bABBABa'
chr5 = 'BAbABAbA'
chr6 = 'BBbaBabaaBabaBaBABA'
chr7 = 'aB'
chr8 = 'BBAbABaBABABABa'
chr9 = 'AB'
chr10 = 'BBABa'
chr11 = 'Aabbb'
chr12 = 'aBAB'
chr13 = 'BaaBaBABaBa'
chr14 = 'B'
chr15 = 'aBaBBABaBaBBAB'
chr16 = 'BABA'
chr17 = 'BBbAaBBbaBaaBBBa'
chr17b = 'aaaBbbbbaaB'
chr18 = 'bBaBABABAba'
chr20 = 'BbAbA'
chrX = 'BAaBbbAaBBba'
chrY = 'AaaaaAAaaAaaAA'
print("pJa conserved")
print("SF1")
for i in [chr1, chr3, chr5, chr6, chr7, chr10, chr12, chr16]:
print("chr <-"+str(re.findall('A+', i)))
print(i.count('A'))
for i in [chr1, chr3, chr5, chr6, chr7, chr10, chr12, chr16]:
print("chr <-"+str(re.findall('B+', i)))
print(i.count('B'))
print("SF2")
for i in [chr2, chr4, chr8, chr9, chr13, chr15, chr18, chr20]:
print("chr <-"+str(re.findall('A+', i)))
print(i.count('A'))
for i in [chr2, chr4, chr8, chr9, chr13, chr15, chr18, chr20]:
print("chr <-"+str(re.findall('B+', i)))
print(i.count('B'))
print("SF3")
for i in [chr11, chr17, chr17b, chrX]:
print("chr <-"+str(re.findall('A+', i)))
print(i.count('A'))
for i in [chr11, chr17, chr17b, chrX]:
print("chr <-"+str(re.findall('B+', i)))
print(i.count('B'))