def test_regionkey(self): for chrom, startpos, endpos, strand, _, _, rk, _, _, _ in regionsTestData: h = variantkey.regionkey(chrom, startpos, endpos, strand) self.assertEqual(h, rk)
vk.extract_regionkey_startpos(0xc80001f400003e84) # 1000 vk.extract_regionkey_endpos(0xc80001f400003e84) # 2000 vk.extract_regionkey_strand(0xc80001f400003e84) # 2 vk.decode_regionkey(0xc80001f400003e84) # (25, 1000, 2000, 2) vk.reverse_regionkey(0xc80001f400003e84) # (b'MT', 1000, 2000, -1) vk.regionkey("MT", 1000, 2000, -1) # 14411520955069251204 vk.extend_regionkey(14411520955069251204, 100) # 14411520740320887204 vk.reverse_regionkey(14411520740320887204) # (b'MT', 900, 2100, -1) vk.regionkey_hex(0xc80001f400003e84) # b'c80001f400003e84' vk.parse_regionkey_hex("c80001f400003e84") # 14411520955069251204 vk.get_regionkey_chrom_startpos(0xc80001f400003e84)
def test_extend_regionkey(self): rk = variantkey.regionkey("X", 10000, 20000, -1) h = variantkey.extend_regionkey(rk, 1000) self.assertEqual(h, 13258616630331740228) h = variantkey.extend_regionkey(rk, 300000000) self.assertEqual(h, 13258597305126223868)
print(vk.extract_regionkey_startpos(0xc80001f400003e84)) # 1000 print(vk.extract_regionkey_endpos(0xc80001f400003e84)) # 2000 print(vk.extract_regionkey_strand(0xc80001f400003e84)) # 2 print(vk.decode_regionkey(0xc80001f400003e84)) # (25, 1000, 2000, 2) print(vk.reverse_regionkey(0xc80001f400003e84)) # (b'MT', 1000, 2000, -1) print(vk.regionkey("MT", 1000, 2000, -1)) # 14411520955069251204 print(vk.regionkey_hex(0xc80001f400003e84)) # b'c80001f400003e84' print(vk.parse_regionkey_hex("c80001f400003e84")) # 14411520955069251204 print(vk.get_regionkey_chrom_startpos(0xc80001f400003e84)) # 6710887400 print(vk.get_regionkey_chrom_endpos(0xc80001f400003e84)) # 6710888400 print(vk.are_overlapping_regions(5, 4, 6, 5, 3, 7))