|
a |
|
b/bin/DeepMod_scripts/myCom.py |
|
|
1 |
|
|
|
2 |
# |
|
|
3 |
# output level |
|
|
4 |
# |
|
|
5 |
OUTPUT_DEBUG=0; |
|
|
6 |
OUTPUT_INFO=1; |
|
|
7 |
OUTPUT_WARNING=2; |
|
|
8 |
OUTPUT_ERROR=3; |
|
|
9 |
|
|
|
10 |
# |
|
|
11 |
# base-pairing |
|
|
12 |
# 4 base types with case insensitive |
|
|
13 |
# |
|
|
14 |
na_bp = {"A":"T", \ |
|
|
15 |
"C":"G", \ |
|
|
16 |
"G":"C", \ |
|
|
17 |
"T":"A", \ |
|
|
18 |
"a":"t", \ |
|
|
19 |
"c":"g", \ |
|
|
20 |
"g":"c", \ |
|
|
21 |
"t":"a", \ |
|
|
22 |
"N":"N", \ |
|
|
23 |
"n":"n" \ |
|
|
24 |
} |
|
|
25 |
|
|
|
26 |
g_ACGT = ['A', 'C', 'G', 'T'] |
|
|
27 |
|
|
|
28 |
acgt = na_bp.keys(); |
|
|
29 |
|
|
|
30 |
# |
|
|
31 |
# get complementary base of a given base |
|
|
32 |
# |
|
|
33 |
def getComplementary(na): |
|
|
34 |
com_na = [] |
|
|
35 |
for i in range(len(na)): |
|
|
36 |
com_na.append(na_bp[na[i]]) |
|
|
37 |
''.join(com_na[::-1]) |
|
|
38 |
|
|
|
39 |
# |
|
|
40 |
# format the end of a folder |
|
|
41 |
# |
|
|
42 |
def format_last_letter_of_folder(cursub): |
|
|
43 |
if not cursub==None: |
|
|
44 |
if cursub[-1]=='/': return cursub; |
|
|
45 |
elif cursub[-1]=='\\': return cursub[:-1]+'/'; |
|
|
46 |
else: return cursub+'/'; |
|
|
47 |
|
|
|
48 |
# |
|
|
49 |
# default path of data in fast5 files |
|
|
50 |
# |
|
|
51 |
analyses_base = "Analyses" |
|
|
52 |
basecall_events_base = "Events" |
|
|
53 |
raw_base = 'Raw' |
|
|
54 |
reads_base = "Reads" |
|
|
55 |
signal_base = "Signal" |
|
|
56 |
basecall_fastq_base = "Fastq" |
|
|
57 |
|
|
|
58 |
|
|
|
59 |
|
|
|
60 |
|