def test_no_event_with_one_two_skipped_exon(self):
self.graph.add_path([str(self.ex1), str(self.ex2),
str(self.ex3), str(self.ex4), str(self.ex5)])
self.graph.add_path([str(self.ex1), str(self.ex2), str(self.ex5)])
self.events = list(find_MXE(self.graph, self.exonsDB))
self.assertEqual(len(self.events), 0)
def test_no_event(self):
pass
self.graph.add_path([str(self.ex1), str(self.ex2), str(self.ex3),
str(self.ex4), str(self.ex5)])
self.events = list(find_MXE(self.graph, self.exonsDB))
self.assertEqual(len(self.events), 0)
def test_one_event_one_mx_exon(self):
self.graph.add_path([str(self.ex1), str(self.ex2),
str(self.ex3), str(self.ex6),
str(self.ex7)])
self.graph.add_path([str(self.ex1), str(self.ex2),
str(self.ex4), str(self.ex6),
str(self.ex7)])
self.events = find_MXE(self.graph, self.exonsDB)
self.assertEqual(len(self.events), 1)
self.assertEqual(len(self.events[0]), 2)
self.assertEqual(self.exonsDB[self.events[0][0][0]].start, 300)
self.assertEqual(self.exonsDB[self.events[0][0][0]].end, 400)
self.assertEqual(self.exonsDB[self.events[0][0][-1]].start, 1100)
self.assertEqual(self.exonsDB[self.events[0][0][-1]].end, 1200)