def test_duplicateannotations_mixed_dups(self):
geno1 = self.createGeno(1)
annot1 = self.createAnnot(10, "term10")
annot2 = self.createAnnot(20, "term20")
annot2_2 = self.createAnnot(20, "term20")
annot3 = self.createAnnot(30, "term30")
ev1 = self.createEvidence(100)
ev2 = self.createEvidence(200)
ev3 = self.createEvidence(300)
ev4 = self.createEvidence(400)
ev5 = self.createEvidence(500)
geno1.mp_headers[0]['annots'].append(annot1)
geno1.mp_headers[0]['annots'].append(annot2)
geno1.mp_headers[0]['annots'].append(annot2_2)
geno1.mp_headers[0]['annots'].append(annot3)
annot1.evidences.append(ev1)
annot2.evidences.append(ev2)
annot2_2.evidences.append(ev3)
annot2_2.evidences.append(ev4)
annot3.evidences.append(ev5)
genotype_mp_hunter._collapseDuplicateAnnotations([geno1])
self.assertEquals([ev1], annot1.evidences)
self.assertEquals([ev2, ev3, ev4], annot2.evidences)
self.assertEquals([ev5], annot3.evidences)
self.assertEquals(3, len(geno1.mp_headers[0]['annots']))
def test_duplicateannotations_mixed_dups(self):
geno1 = self.createGeno(1)
annot1 = self.createAnnot(10, "term10")
annot2 = self.createAnnot(20, "term20")
annot2_2 = self.createAnnot(20, "term20")
annot3 = self.createAnnot(30, "term30")
ev1 = self.createEvidence(100)
ev2 = self.createEvidence(200)
ev3 = self.createEvidence(300)
ev4 = self.createEvidence(400)
ev5 = self.createEvidence(500)
geno1.mp_headers[0]['annots'].append(annot1)
geno1.mp_headers[0]['annots'].append(annot2)
geno1.mp_headers[0]['annots'].append(annot2_2)
geno1.mp_headers[0]['annots'].append(annot3)
annot1.evidences.append(ev1)
annot2.evidences.append(ev2)
annot2_2.evidences.append(ev3)
annot2_2.evidences.append(ev4)
annot3.evidences.append(ev5)
genotype_mp_hunter._collapseDuplicateAnnotations([geno1])
self.assertEquals([ev1], annot1.evidences)
self.assertEquals([ev2,ev3,ev4], annot2.evidences)
self.assertEquals([ev5], annot3.evidences)
self.assertEquals(3, len(geno1.mp_headers[0]['annots']))
def test_duplicateannotations_only_dups(self):
geno1 = self.createGeno(1)
annot1 = self.createAnnot(10, "term10")
annot1_1 = self.createAnnot(10, "term10")
ev1 = self.createEvidence(100)
ev2 = self.createEvidence(200)
geno1.mp_headers[0]['annots'].append(annot1)
geno1.mp_headers[0]['annots'].append(annot1_1)
annot1.evidences.append(ev1)
annot1_1.evidences.append(ev2)
genotype_mp_hunter._collapseDuplicateAnnotations([geno1])
self.assertEquals([ev1, ev2], annot1.evidences)
self.assertEquals(1, len(geno1.mp_headers[0]['annots']))
def test_duplicateannotations_only_dups(self):
geno1 = self.createGeno(1)
annot1 = self.createAnnot(10, "term10")
annot1_1 = self.createAnnot(10, "term10")
ev1 = self.createEvidence(100)
ev2 = self.createEvidence(200)
geno1.mp_headers[0]['annots'].append(annot1)
geno1.mp_headers[0]['annots'].append(annot1_1)
annot1.evidences.append(ev1)
annot1_1.evidences.append(ev2)
genotype_mp_hunter._collapseDuplicateAnnotations([geno1])
self.assertEquals([ev1,ev2], annot1.evidences)
self.assertEquals(1, len(geno1.mp_headers[0]['annots']))
