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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
BT=${BT-../../bin/bedtools}
check()
{
if diff $1 $2; then
echo ok
return 1
else
echo fail
return 0
fi
}
###########################################################
# test baseline complement
###########################################################
echo " complement.t1...\c"
echo "chr1 1 20" > exp
$BT complement -i <(echo -e "chr1\t0\t1") \
-g <(echo -e "chr1\t20") \
> obs
check obs exp
rm obs exp
###########################################################
# ends are covered
###########################################################
echo " complement.t2...\c"
echo "chr1 1 19" > exp
$BT complement -i <(echo -e "chr1\t0\t1\nchr1\t19\t20") \
-g <(echo -e "chr1\t20") \
> obs
check obs exp
rm obs exp
###########################################################
# middle is covered
###########################################################
echo " complement.t3...\c"
echo "chr1 0 10
chr1 15 20" > exp
$BT complement -i <(echo -e "chr1\t10\t15") \
-g <(echo -e "chr1\t20") \
> obs
check obs exp
rm obs exp
###########################################################
# entirety is covered
###########################################################
echo " complement.t4...\c"
touch exp
$BT complement -i <(echo -e "chr1\t0\t20") \
-g <(echo -e "chr1\t20") \
> obs
check obs exp
rm obs exp
###########################################################
# nothing is covered
###########################################################
echo " complement.t5...\c"
echo "chr1 0 20" > exp
$BT complement -i <(echo -e "chr2\t0\t20") \
-g <(echo -e "chr1\t20\nchr2\t20") \
> obs
check obs exp
rm obs exp
###########################################################
# Issue #356
###########################################################
echo " complement.t6...\c"
echo "chr1 10000 249250621" > exp
$BT complement -i <(echo -e "chr1\t0\t10000\ttelomere") \
-g <(echo -e "chr1\t249250621") \
> obs
check obs exp
rm obs exp
###########################################################
# Multiple chroms
###########################################################
echo " complement.t7...\c"
echo "chr1 0 10
chr2 0 10" > exp
$BT complement -i <(echo -e "chr1\t10\t20\nchr2\t10\t20") \
-g <(echo -e "chr1\t20\nchr2\t20") \
> obs
check obs exp
rm obs exp
###########################################################
# Multiple chroms, chr1 is covered
###########################################################
echo " complement.t8...\c"
echo "chr2 0 10" > exp
$BT complement -i <(echo -e "chr1\t0\t20\nchr2\t10\t20") \
-g <(echo -e "chr1\t20\nchr2\t20") \
> obs
check obs exp
rm obs exp
###########################################################
# record exceeds chrom length
###########################################################
echo " complement.t9...\c"
echo -e "***** WARNING: chr1:90-110 exceeds the length of chromosome (chr1)\nchr1\t0\t90" > exp
$BT complement -i <(echo -e "chr1\t90\t110") \
-g <(echo -e "chr1\t100") \
&> obs
check obs exp
rm obs exp
|
