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))
