def test_is_artifact(self): artifact = Cluster(read=PK, write=PK, rounding=True) artifact.read_line('chr1 1 111 5:1|100:4|6:2') almost_not_artifact = Cluster(read=PK, write=PK, rounding=True) almost_not_artifact.read_line('chr1 1 100 35:1|30:4|35:1') almost_artifact = Cluster(read=PK, write=PK, rounding=True) almost_artifact.read_line('chr1 1 99 45:1|25:4|30:1') almost_artifact._recalculate_end() self.assertTrue(artifact.is_artifact()) self.assertFalse(almost_not_artifact.is_artifact()) #with the changes for version 0.6.1, this line IS an artifact self.assertFalse(almost_artifact.is_artifact())
def test_bug_is_artifact(self): cluster = Cluster(rounding=True) cluster.read_line('chr1 10505760 10505928 8:4|3:5|5:6|8:8|2:10|1:11|5:12|3:15|2:16|2:16|3:18|7:20|3:19|57:20|21:19|8:15|3:14|5:13|8:11|2:9|1:8|5:7|3:4|2:3|2:3') self.assertTrue(cluster.is_artifact(0.3))