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