def test_get_profile(self):
r = Region(start=1, end=1999)
c = Cluster(read=BED)
c.read_line('chr4 1 40')
r.add_tags(c, True)
c = Cluster(read=BED, read_half_open=True)
c.read_line('chr4 400 500')
r.add_tags(c, True)
meta = r.get_metacluster()
self.assertEqual(meta._levels, [[40, 1.0], [360, 0.0], [100, 1.0]])
def test_subtract_with_gaps(self):
"""
238c238
< chr17 147132 147296 5:1.00|22:2.00|37:1.00|37:2.00|64:1.00 2.0 . 147180 224.0
---
> chr17 147132 147260 5:1.00|22:2.00|37:1.00|37:2.00|28:1.00 2.0 . 147180 188.0
experiment
chr17 146970 147070
chr17 147058 147158
chr17 147132 147232
chr17 147196 147296
chr17 147303 147403
chr17 147445 147545
chr17 147461 147561
control
chr17 146904 147004
chr17 147036 147136
chr17 147261 147361
chr17 147437 147537
chr17 147472 147572
chr17 147472 147572
"""
exp_cached = Cluster(cached=True)
control0 = Cluster(read=BED, cached=True)
control1 = Cluster(read=BED, cached=True)
control2 = Cluster(read=BED, cached=True)
control3 = Cluster(read=BED, cached=True)
control4 = Cluster(read=BED, cached=True)
#experiment.bed
exp_cached.read_line("chr17 146970 147296 88:1.00|13:2.00|61:1.00|27:2.00|37:1.00|37:2.00|64:1.00 2.0 . 147139 404.0")
#print exp_cached.write_line()
#control.bed
r = Region("chr17", 1, 2000)
control0.read_line("chr17 146904 147004")
control1.read_line("chr17 147036 147136")
control2.read_line("chr17 147261 147361")
control3.read_line("chr17 147437 147537")
control4.read_line("chr17 147472 147572")
r.add_tags([control1, control2])
meta = r.get_metacluster()
#print meta._profile
exp_cached -= meta
