def test_no_retained_introns(self):
self.path1 = [str(self.ex1), str(self.ex2), str(self.ex3)]
self.path2 = [str(self.ex1), str(self.ex3), str(self.ex4)]
self.graph.add_path(self.path1)
self.graph.add_path(self.path2)
self.events = list(find_RI(self.graph, self.tree, self.exonsDB))
self.assertEqual(len(self.events), 0)
def test_one_retained_introns(self):
self.ex5 = Exon('chrX', 3000, 6000, 'ex1.1', '+')
self.exonsDB[str(self.ex5)] = self.ex5
self.tree.add_interval(self.ex5)
self.path1 = [str(self.ex1), str(self.ex2),
str(self.ex3), str(self.ex4)]
self.path2 = [str(self.ex1), str(self.ex5), str(self.ex4)]
self.graph.add_path(self.path1)
self.graph.add_path(self.path2)
self.events = list(find_RI(self.graph, self.tree, self.exonsDB))
self.assertEqual(len(self.events), 1)