diff --git a/data/aluY.chr1.bed.gz b/data/aluY.chr1.bed.gz new file mode 100644 index 0000000000000000000000000000000000000000..8e03a623eb8d50da534a93fb70855a32a2c5998a Binary files /dev/null and b/data/aluY.chr1.bed.gz differ diff --git a/data/gerp.chr1.bed.gz b/data/gerp.chr1.bed.gz new file mode 100644 index 0000000000000000000000000000000000000000..d8a7e6b7144d8d56c179ad74b9f920e21ec3d1cb Binary files /dev/null and b/data/gerp.chr1.bed.gz differ diff --git a/data/refseq.chr1.exons.bed.gz b/data/refseq.chr1.exons.bed.gz new file mode 100644 index 0000000000000000000000000000000000000000..adbc418946dad20bd0cea89a7688370ebec3c76f Binary files /dev/null and b/data/refseq.chr1.exons.bed.gz differ diff --git a/data/simpleRepeats.chr1.bed.gz b/data/simpleRepeats.chr1.bed.gz new file mode 100644 index 0000000000000000000000000000000000000000..6beb2708f4cfc4e723bab8dcaf1482aca1fcaee8 Binary files /dev/null and b/data/simpleRepeats.chr1.bed.gz differ diff --git a/test/reldist/test-reldist.sh b/test/reldist/test-reldist.sh index 1d7d9669061cd1b22b46b8ecef4f1f80c9cc95d7..eebde8bda5c3608f9477bd0e91b2028d627e131b 100644 --- a/test/reldist/test-reldist.sh +++ b/test/reldist/test-reldist.sh @@ -16,9 +16,9 @@ check() echo " reldist.t01...\c" echo \ "reldist count total fraction -0.00 420889 420889 1.000" > exp -$BT reldist -a ../../data/refseq.hg19.exons.bed.gz \ - -b ../../data/refseq.hg19.exons.bed.gz > obs +0.00 43424 43424 1.000" > exp +$BT reldist -a ../../data/refseq.chr1.exons.bed.gz \ + -b ../../data/refseq.chr1.exons.bed.gz > obs check obs exp rm obs exp @@ -29,59 +29,58 @@ rm obs exp echo " reldist.t02...\c" echo \ "reldist count total fraction -0.00 1930 420584 0.005 -0.01 4937 420584 0.012 -0.02 6021 420584 0.014 -0.03 6774 420584 0.016 -0.04 7403 420584 0.018 -0.05 7152 420584 0.017 -0.06 7824 420584 0.019 -0.07 7818 420584 0.019 -0.08 7679 420584 0.018 -0.09 8275 420584 0.020 -0.10 7886 420584 0.019 -0.11 8308 420584 0.020 -0.12 8335 420584 0.020 -0.13 8532 420584 0.020 -0.14 8486 420584 0.020 -0.15 8335 420584 0.020 -0.16 8207 420584 0.020 -0.17 8674 420584 0.021 -0.18 8541 420584 0.020 -0.19 8795 420584 0.021 -0.20 8718 420584 0.021 -0.21 8731 420584 0.021 -0.22 8464 420584 0.020 -0.23 8627 420584 0.021 -0.24 8766 420584 0.021 -0.25 8927 420584 0.021 -0.26 8860 420584 0.021 -0.27 8876 420584 0.021 -0.28 9285 420584 0.022 -0.29 8882 420584 0.021 -0.30 9092 420584 0.022 -0.31 8882 420584 0.021 -0.32 8887 420584 0.021 -0.33 9296 420584 0.022 -0.34 8895 420584 0.021 -0.35 8964 420584 0.021 -0.36 9106 420584 0.022 -0.37 9377 420584 0.022 -0.38 9199 420584 0.022 -0.39 9240 420584 0.022 -0.40 9268 420584 0.022 -0.41 9443 420584 0.022 -0.42 9541 420584 0.023 -0.43 9287 420584 0.022 -0.44 8998 420584 0.021 -0.45 9238 420584 0.022 -0.46 8892 420584 0.021 -0.47 8859 420584 0.021 -0.48 9033 420584 0.021 -0.49 9028 420584 0.021 -0.50 11 420584 0.000" > exp -$BT reldist -a ../../data/refseq.hg19.exons.bed.gz \ - -b ../../data/aluY.hg19.bed.gz > obs +0.00 164 43408 0.004 +0.01 551 43408 0.013 +0.02 598 43408 0.014 +0.03 637 43408 0.015 +0.04 793 43408 0.018 +0.05 688 43408 0.016 +0.06 874 43408 0.020 +0.07 765 43408 0.018 +0.08 685 43408 0.016 +0.09 929 43408 0.021 +0.10 876 43408 0.020 +0.11 959 43408 0.022 +0.12 860 43408 0.020 +0.13 851 43408 0.020 +0.14 903 43408 0.021 +0.15 893 43408 0.021 +0.16 883 43408 0.020 +0.17 828 43408 0.019 +0.18 917 43408 0.021 +0.19 875 43408 0.020 +0.20 897 43408 0.021 +0.21 986 43408 0.023 +0.22 903 43408 0.021 +0.23 944 43408 0.022 +0.24 904 43408 0.021 +0.25 867 43408 0.020 +0.26 943 43408 0.022 +0.27 933 43408 0.021 +0.28 1132 43408 0.026 +0.29 881 43408 0.020 +0.30 851 43408 0.020 +0.31 963 43408 0.022 +0.32 950 43408 0.022 +0.33 965 43408 0.022 +0.34 907 43408 0.021 +0.35 884 43408 0.020 +0.36 965 43408 0.022 +0.37 944 43408 0.022 +0.38 911 43408 0.021 +0.39 939 43408 0.022 +0.40 921 43408 0.021 +0.41 950 43408 0.022 +0.42 935 43408 0.022 +0.43 919 43408 0.021 +0.44 915 43408 0.021 +0.45 934 43408 0.022 +0.46 843 43408 0.019 +0.47 850 43408 0.020 +0.48 1006 43408 0.023 +0.49 937 43408 0.022" > exp +$BT reldist -a ../../data/refseq.chr1.exons.bed.gz \ + -b ../../data/aluY.chr1.bed.gz > obs check obs exp rm obs exp @@ -93,58 +92,58 @@ rm obs exp echo " reldist.t03...\c" echo \ "reldist count total fraction -0.00 156776 396736 0.395 -0.01 33579 396736 0.085 -0.02 17808 396736 0.045 -0.03 11826 396736 0.030 -0.04 8961 396736 0.023 -0.05 7302 396736 0.018 -0.06 5994 396736 0.015 -0.07 5357 396736 0.014 -0.08 5032 396736 0.013 -0.09 4415 396736 0.011 -0.10 4411 396736 0.011 -0.11 4062 396736 0.010 -0.12 3741 396736 0.009 -0.13 3760 396736 0.009 -0.14 3501 396736 0.009 -0.15 3394 396736 0.009 -0.16 3539 396736 0.009 -0.17 3442 396736 0.009 -0.18 3500 396736 0.009 -0.19 3072 396736 0.008 -0.20 3713 396736 0.009 -0.21 3249 396736 0.008 -0.22 3328 396736 0.008 -0.23 3271 396736 0.008 -0.24 3017 396736 0.008 -0.25 3756 396736 0.009 -0.26 3242 396736 0.008 -0.27 3057 396736 0.008 -0.28 3519 396736 0.009 -0.29 3222 396736 0.008 -0.30 3529 396736 0.009 -0.31 3284 396736 0.008 -0.32 3387 396736 0.009 -0.33 3231 396736 0.008 -0.34 3263 396736 0.008 -0.35 3421 396736 0.009 -0.36 3529 396736 0.009 -0.37 3376 396736 0.009 -0.38 3414 396736 0.009 -0.39 3183 396736 0.008 -0.40 3855 396736 0.010 -0.41 3402 396736 0.009 -0.42 3803 396736 0.010 -0.43 3458 396736 0.009 -0.44 3594 396736 0.009 -0.45 3654 396736 0.009 -0.46 3489 396736 0.009 -0.47 3394 396736 0.009 -0.48 3862 396736 0.010 -0.49 2957 396736 0.007 -0.50 805 396736 0.002" > exp -$BT reldist -a ../../data/refseq.hg19.exons.bed.gz \ - -b ../../data/29way_pi.conservation.hg19.bed.gz > obs +0.00 20629 43422 0.475 +0.01 2629 43422 0.061 +0.02 1427 43422 0.033 +0.03 985 43422 0.023 +0.04 897 43422 0.021 +0.05 756 43422 0.017 +0.06 667 43422 0.015 +0.07 557 43422 0.013 +0.08 603 43422 0.014 +0.09 487 43422 0.011 +0.10 461 43422 0.011 +0.11 423 43422 0.010 +0.12 427 43422 0.010 +0.13 435 43422 0.010 +0.14 375 43422 0.009 +0.15 367 43422 0.008 +0.16 379 43422 0.009 +0.17 371 43422 0.009 +0.18 346 43422 0.008 +0.19 389 43422 0.009 +0.20 377 43422 0.009 +0.21 411 43422 0.009 +0.22 377 43422 0.009 +0.23 352 43422 0.008 +0.24 334 43422 0.008 +0.25 315 43422 0.007 +0.26 370 43422 0.009 +0.27 330 43422 0.008 +0.28 330 43422 0.008 +0.29 280 43422 0.006 +0.30 309 43422 0.007 +0.31 326 43422 0.008 +0.32 287 43422 0.007 +0.33 294 43422 0.007 +0.34 306 43422 0.007 +0.35 307 43422 0.007 +0.36 309 43422 0.007 +0.37 271 43422 0.006 +0.38 293 43422 0.007 +0.39 311 43422 0.007 +0.40 331 43422 0.008 +0.41 320 43422 0.007 +0.42 299 43422 0.007 +0.43 327 43422 0.008 +0.44 321 43422 0.007 +0.45 326 43422 0.008 +0.46 306 43422 0.007 +0.47 354 43422 0.008 +0.48 365 43422 0.008 +0.49 336 43422 0.008 +0.50 38 43422 0.001" > exp +$BT reldist -a ../../data/refseq.chr1.exons.bed.gz \ + -b ../../data/gerp.chr1.bed.gz > obs check obs exp rm obs exp \ No newline at end of file