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
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
set -e;
BT=${BT-../../bin/bedtools}
FAILURES=0;
check()
{
if diff $1 $2; then
echo ok
else
FAILURES=$(expr $FAILURES + 1);
echo fail
fi
}
# randData.bed was generated on OSX Sierra using the command below. Different
# systems may generate different results, so the file required for the tests is
# stored in the repo.
cat > /dev/null << EOCOMMENT
cat /dev/zero |tr "\0" "\n" | head -n 10000 |\
awk '{S=int(rand()*1000.0);E=S+int(rand()*100000); printf("chrX\t%d\t%d\n",S,E);}' > tmp.bed
EOCOMMENT
echo -e " split.01.size...\c"
rm -f _tmp.*.bed
echo "_tmp.00001.bed 9943540 200
_tmp.00002.bed 9943482 201
_tmp.00003.bed 9943541 200
_tmp.00004.bed 9943561 200
_tmp.00005.bed 9943471 200
_tmp.00006.bed 9943475 200
_tmp.00007.bed 9943468 200
_tmp.00008.bed 9943487 200
_tmp.00009.bed 9943539 200
_tmp.00010.bed 9943531 200" > exp
${BT} split -i randData.bed -p _tmp -n 50 -a size | head > _tmp.size.tsv
check exp _tmp.size.tsv
rm exp
echo -e " split.02.simple...\c"
rm -f _tmp.*.bed
echo "_tmp.00001.bed 9952674 200
_tmp.00002.bed 9751661 200
_tmp.00003.bed 9649058 200
_tmp.00004.bed 9929508 200
_tmp.00005.bed 9556713 200
_tmp.00006.bed 10298876 200
_tmp.00007.bed 10043102 200
_tmp.00008.bed 9781861 200
_tmp.00009.bed 9502188 200
_tmp.00010.bed 9991229 200" > exp
${BT} split -i randData.bed -p _tmp -n 50 -a simple | head > _tmp.simple.tsv
check exp _tmp.simple.tsv
rm exp
echo -e " spliit.03.simple...\c"
rm -f _tmp.*.bed
echo "_tmp.00001.bed 414150 10
_tmp.00002.bed 586843 10
_tmp.00003.bed 503604 10
_tmp.00004.bed 410044 10
_tmp.00005.bed 499400 10
_tmp.00006.bed 537341 10
_tmp.00007.bed 698581 10
_tmp.00008.bed 555258 10
_tmp.00009.bed 474511 10
_tmp.00010.bed 633012 10" > exp
${BT} split -i randData.bed -p _tmp -n 1000 -a simple | head > _tmp.simple.tsv
check exp _tmp.simple.tsv
rm exp
rm -f _tmp.*.bed jeter.bed _tmp.simple.tsv _tmp.size.tsv
[[ ! -f tmp.bed ]] || rm tmp.bed
[[ $FAILURES -eq 0 ]] || exit 1;
|
