7
7
class pyGenoSNPTests (unittest .TestCase ):
8
8
9
9
def setUp (self ):
10
+ try :
11
+ B .importGenome ("Human.GRCh37.75_Y-Only.tar.gz" )
12
+ except KeyError :
13
+ print "--> Seems to already exist in db"
14
+
15
+ try :
16
+ B .importSNPs ("Human_agnostic.dummySRY.tar.gz" )
17
+ except KeyError :
18
+ print "--> Seems to already exist in db"
10
19
self .ref = Genome (name = 'GRCh37.75_Y-Only' )
11
20
12
21
def tearDown (self ):
13
22
pass
14
23
15
- #@unittest.skip("skipping")
24
+ #@unittest.skip("skipping")
16
25
def test_vanilla (self ) :
17
26
dummy = Genome (name = 'GRCh37.75_Y-Only' , SNPs = 'dummySRY_AGN' )
18
27
persProt = dummy .get (Protein , id = 'ENSP00000438917' )[0 ]
19
28
refProt = self .ref .get (Protein , id = 'ENSP00000438917' )[0 ]
20
-
29
+
21
30
self .assertEqual ('ATGCAATCATATGCTTCTGC' , refProt .transcript .cDNA [:20 ])
22
31
self .assertEqual ('HTGCAATCATATGCTTCTGC' , persProt .transcript .cDNA [:20 ])
23
32
24
- #@unittest.skip("skipping")
33
+ #@unittest.skip("skipping")
25
34
def test_noModif (self ) :
26
35
from pyGeno .SNPFiltering import SNPFilter
27
36
@@ -37,7 +46,7 @@ def filter(self, chromosome, dummySRY_AGN) :
37
46
refProt = self .ref .get (Protein , id = 'ENSP00000438917' )[0 ]
38
47
39
48
self .assertEqual (persProt .transcript .cDNA [:20 ], refProt .transcript .cDNA [:20 ])
40
-
49
+
41
50
def test_insert (self ) :
42
51
from pyGeno .SNPFiltering import SNPFilter
43
52
@@ -47,7 +56,7 @@ def __init__(self) :
47
56
48
57
def filter (self , chromosome , dummySRY_AGN ) :
49
58
from pyGeno .SNPFiltering import SequenceInsert
50
-
59
+
51
60
refAllele = chromosome .refSequence [dummySRY_AGN .start ]
52
61
return SequenceInsert ('TCA' )
53
62
@@ -58,7 +67,7 @@ def filter(self, chromosome, dummySRY_AGN) :
58
67
self .assertEqual ('ATGCAATCATATGCTTCTGC' , refProt .transcript .cDNA [:20 ])
59
68
self .assertEqual ('TGAATGCAATCATATGCTTC' , persProt .transcript .cDNA [:20 ])
60
69
61
- #@unittest.skip("skipping")
70
+ #@unittest.skip("skipping")
62
71
def test_SNP (self ) :
63
72
from pyGeno .SNPFiltering import SNPFilter
64
73
@@ -68,17 +77,17 @@ def __init__(self) :
68
77
69
78
def filter (self , chromosome , dummySRY_AGN ) :
70
79
from pyGeno .SNPFiltering import SequenceSNP
71
-
80
+
72
81
return SequenceSNP (dummySRY_AGN .alt )
73
-
82
+
74
83
dummy = Genome (name = 'GRCh37.75_Y-Only' , SNPs = 'dummySRY_AGN' , SNPFilter = MyFilter ())
75
84
persProt = dummy .get (Protein , id = 'ENSP00000438917' )[0 ]
76
-
85
+
77
86
refProt = self .ref .get (Protein , id = 'ENSP00000438917' )[0 ]
78
87
self .assertEqual ('M' , refProt .sequence [0 ])
79
88
self .assertEqual ('L' , persProt .sequence [0 ])
80
89
81
- #@unittest.skip("skipping")
90
+ #@unittest.skip("skipping")
82
91
def test_deletion (self ) :
83
92
from pyGeno .SNPFiltering import SNPFilter
84
93
@@ -88,35 +97,35 @@ def __init__(self) :
88
97
89
98
def filter (self , chromosome , dummySRY_AGN ) :
90
99
from pyGeno .SNPFiltering import SequenceDel
91
-
100
+
92
101
refAllele = chromosome .refSequence [dummySRY_AGN .start ]
93
102
return SequenceDel (1 )
94
-
103
+
95
104
dummy = Genome (name = 'GRCh37.75_Y-Only' , SNPs = 'dummySRY_AGN' , SNPFilter = MyFilter ())
96
105
persProt = dummy .get (Protein , id = 'ENSP00000438917' )[0 ]
97
106
refProt = self .ref .get (Protein , id = 'ENSP00000438917' )[0 ]
98
107
99
108
self .assertEqual ('ATGCAATCATATGCTTCTGC' , refProt .transcript .cDNA [:20 ])
100
109
self .assertEqual ('TGCAATCATATGCTTCTGCT' , persProt .transcript .cDNA [:20 ])
101
110
102
- #@unittest.skip("skipping")
111
+ #@unittest.skip("skipping")
103
112
def test_bags (self ) :
104
113
dummy = Genome (name = 'GRCh37.75_Y-Only' )
105
114
self .assertEqual (dummy .wrapped_object , self .ref .wrapped_object )
106
-
107
- #@unittest.skip("skipping")
115
+
116
+ #@unittest.skip("skipping")
108
117
def test_prot_find (self ) :
109
118
prot = self .ref .get (Protein , id = 'ENSP00000438917' )[0 ]
110
119
needle = prot .sequence [:10 ]
111
120
self .assertEqual (0 , prot .find (needle ))
112
121
needle = prot .sequence [- 10 :]
113
122
self .assertEqual (len (prot )- 10 , prot .find (needle ))
114
123
115
- #@unittest.skip("skipping")
124
+ #@unittest.skip("skipping")
116
125
def test_trans_find (self ) :
117
126
trans = self .ref .get (Transcript , name = "SRY-001" )[0 ]
118
127
self .assertEqual (0 , trans .find (trans [:5 ]))
119
-
128
+
120
129
# @unittest.skip("remote server down")
121
130
# def test_import_remote_genome(self) :
122
131
# self.assertRaises(KeyError, B.importRemoteGenome, "Human.GRCh37.75_Y-Only.tar.gz")
@@ -130,13 +139,13 @@ def runTests() :
130
139
B .importGenome ("Human.GRCh37.75_Y-Only.tar.gz" )
131
140
except KeyError :
132
141
print "--> Seems to already exist in db"
133
-
142
+
134
143
try :
135
144
B .importSNPs ("Human_agnostic.dummySRY.tar.gz" )
136
145
except KeyError :
137
146
print "--> Seems to already exist in db"
138
147
# import time
139
- # time.sleep(10)
148
+ # time.sleep(10)
140
149
unittest .main ()
141
150
142
151
if __name__ == "__main__" :
0 commit comments