def test_add_to_read_groups_doesNothingWhenNotInCommonKeys(self): readA = SplitRead(**initParams({'name':'readA'})) common_keys = set([readA.key()]) read_groups = {} readB = SplitRead(**initParams({'name':'readB'})) readB.add_to_read_groups(common_keys, read_groups) self.assertEqual(0, len(read_groups))
def test_add_to_read_groups_addsToCorrectSide(self): left10 = SplitRead(**initParams({'name':'readA', 'side':"L", 'position':10, 'split_len':40, 'original_read_len':100})) left15 = SplitRead(**initParams({'name':'readA', 'side':"L", 'position':15, 'split_len':40, 'original_read_len':100})) right30 = SplitRead(**initParams({'name':'readA', 'side':"R", 'position':30, 'split_len':60, 'original_read_len':100})) read_group_key = left10.key() common_keys = set([read_group_key]) read_groups = {} left10.add_to_read_groups(common_keys, read_groups) left15.add_to_read_groups(common_keys, read_groups) right30.add_to_read_groups(common_keys, read_groups) self.assertEqual(1, len(read_groups)) self.assertEqual(([left10, left15], [right30]), read_groups[read_group_key])