1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
# gene sequence example from @yoch, see
# https://github.com/ilanschnell/bitarray/pull/54
from random import choice
from timeit import timeit
from bitarray import bitarray
trans = {
"A": bitarray("00"),
"T": bitarray("01"),
"G": bitarray("10"),
"C": bitarray("11")
}
N = 10_000
seq = [choice("ATGC") for _ in range(N)]
arr = bitarray()
arr.encode(trans, seq)
assert list(arr.decode(trans)) == seq
# decodage
t = timeit(lambda: list(arr.decode(trans)), number=1000)
print(t)
|
