[0473b3]: / tests / test_janggu_trim.py

Download this file

28 lines (19 with data), 732 Bytes

 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
import os
from pkg_resources import resource_filename
from janggu.data import GenomicIndexer
from janggu.janggutrim import trim_bed
def test_create_from_array(tmpdir):
inbed = resource_filename('janggu', 'resources/bed_test.bed')
outbed = os.path.join(tmpdir.strpath, 'out.bed')
trim_bed(inbed, outbed, 5)
# original file
gindexer = GenomicIndexer.create_from_file(inbed, None, None)
reg = gindexer[0]
assert (reg.start % 5) == 0
assert (reg.end % 5) > 0
# trimmed file
gindexer = GenomicIndexer.create_from_file(outbed, None, None)
gindexer = GenomicIndexer.create_from_file(outbed, None, None)
reg = gindexer[0]
assert (reg.start % 5) == 0
assert (reg.end % 5) == 0