Пример #1
0
 def test_regionkey(self):
     for chrom, startpos, endpos, strand, _, _, rk, _, _, _ in regionsTestData:
         h = variantkey.regionkey(chrom, startpos, endpos, strand)
         self.assertEqual(h, rk)
Пример #2
0
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)
Пример #3
0
 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)
Пример #4
0
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))