1 new commit in galaxy-central: https://bitbucket.org/galaxy/galaxy-central/changeset/0622b489a0f1/ changeset: 0622b489a0f1 user: guru date: 2011-11-03 15:49:45 summary: Added 3 tools: 1) to compute weighted average values based on interval intersections (under Regional variation), 2) to perform logistic regression (under Multiple regression), 3) to compute partial R-squared (under Multiple regression). affected #: 13 files diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 test-data/interpolate_result.bed --- /dev/null +++ b/test-data/interpolate_result.bed @@ -0,0 +1,219 @@ +chr1 27800000 29750669 6600000 0 1.7386 +chr1 29800669 34400000 6600000 0 1.7920 +chr1 46500000 51300000 4800000 0 1.4004 +chr1 102000000 103636494 5000000 0 1.2007 +chr1 103686494 103843838 5000000 0 NA +chr1 103893838 107000000 5000000 0 1.1883 +chr1 117600000 120498679 3100000 0 1.4007 +chr1 120548679 120700000 3100000 0 NA +chr1 153300000 154800000 1500000 0 1.9511 +chr1 163800000 171200000 7400000 0 1.4719 +chr1 174300000 184000000 9700000 0 1.5022 +chr10 6700000 12300000 5600000 0 1.5533 +chr10 42100000 45746970 11200000 0 1.4905 +chr10 45896970 46849175 11200000 0 NA +chr10 46999175 47262482 11200000 0 1.4600 +chr10 47412482 47575713 11200000 0 NA +chr10 47725713 48715542 11200000 0 1.3900 +chr10 48865542 50807416 11200000 0 1.7448 +chr10 50857416 51068851 11200000 0 NA +chr10 51118851 53300000 11200000 0 1.3877 +chr10 74600000 81231464 15000000 0 1.4494 +chr10 81241464 89600000 15000000 0 1.4475 +chr10 98000000 99400000 1400000 0 1.7700 +chr10 102000000 111800000 9800000 0 1.4549 +chr10 114900000 119100000 4200000 0 1.3184 +chr10 127400000 128606059 7974737 0 1.2600 +chr10 128656059 133271394 7974737 0 1.3583 +chr10 133281394 133527517 7974737 0 NA +chr10 133577517 135374737 7974737 0 1.4720 +chr11 50000 1152759 2800000 0 1.5200 +chr11 1169335 2800000 2800000 0 1.6101 +chr11 12600000 16100000 3500000 0 1.2212 +chr11 21600000 26000000 4400000 0 1.3320 +chr11 27200000 31000000 3800000 0 1.2483 +chr11 48800000 50740429 2600000 0 1.3100 +chr11 50947429 51400000 2600000 0 NA +chr11 56400000 68846377 12800000 0 1.7959 +chr11 68848982 69200000 12800000 0 1.6400 +chr11 70700000 85300000 14600000 0 1.5591 +chr11 87900000 92300000 4400000 0 1.2134 +chr11 110000000 115400000 5400000 0 1.4715 +chr11 120700000 134452384 13752384 0 1.4668 +chr12 36500000 44600000 8100000 0 1.3403 +chr12 53100000 56300000 3200000 0 1.6808 +chr12 69800000 74791940 8900000 0 1.2548 +chr12 75041940 78700000 8900000 0 1.2656 +chr12 91200000 94800000 3600000 0 1.6655 +chr12 100000000 107500000 7500000 0 1.6243 +chr13 18400000 32900000 14500000 0 1.5346 +chr13 34700000 39500000 4800000 0 1.1758 +chr13 52200000 57600000 5400000 0 1.2506 +chr13 60500000 72100000 11600000 0 1.2112 +chr13 100500000 109100000 8600000 0 1.3144 +chr14 23600000 31800000 8200000 0 1.3057 +chr14 36900000 41000000 4100000 0 1.2309 +chr14 43200000 48300000 5100000 0 1.2680 +chr14 69300000 72900000 3600000 0 1.6758 +chr14 78400000 88900000 10500000 0 1.2495 +chr15 23300000 24687593 2400000 0 1.2379 +chr15 24731593 25169606 2400000 0 1.2200 +chr15 25270606 25700000 2400000 0 NA +chr15 31400000 37900000 6500000 0 1.4543 +chr15 65300000 70400000 5100000 0 1.6950 +chr16 0 8576922 14700000 0 1.5411 +chr16 8594422 14700000 14700000 0 1.5935 +chr16 16700000 21700000 5000000 0 1.6366 +chr16 69400000 78200000 8800000 0 1.4012 +chr16 80500000 82700000 2200000 0 1.5200 +chr17 0 296854 11200000 0 NA +chr17 343376 11200000 11200000 0 1.7399 +chr17 28800000 31699961 6600000 0 1.4123 +chr17 31799961 33429230 6600000 0 1.7100 +chr17 33529230 35400000 6600000 0 1.8415 +chr17 47600000 54900000 7300000 0 1.4561 +chr17 55600000 59900000 4300000 0 1.5958 +chr18 17300000 31000000 13700000 0 1.3855 +chr18 35500000 50313134 16500000 0 1.3679 +chr18 50360134 52000000 16500000 0 1.3600 +chr19 11000 7297004 12600000 0 1.8232 +chr19 7302004 8593198 12600000 0 2.1300 +chr19 8598198 12600000 12600000 0 1.7010 +chr2 12800000 16171403 4200000 0 1.2641 +chr2 16221403 17000000 4200000 0 1.2505 +chr2 19100000 21012571 4800000 0 1.2500 +chr2 21037571 23900000 4800000 0 1.4405 +chr2 47600000 52700000 5100000 0 1.2384 +chr2 54800000 61100000 6300000 0 1.4715 +chr2 75400000 83700000 8300000 0 1.1639 +chr2 102100000 108600000 6500000 0 1.3381 +chr2 113800000 118600000 4800000 0 1.1949 +chr2 129600000 134800000 5200000 0 1.3332 +chr2 136600000 144700000 8100000 0 1.3447 +chr2 189100000 197100000 8000000 0 1.4590 +chr2 209100000 215100000 6000000 0 1.3963 +chr2 221300000 233711985 15700000 0 1.3701 +chr2 233731985 237000000 15700000 0 1.4560 +chr20 5000000 9000000 4000000 0 1.3515 +chr20 11900000 17800000 5900000 0 1.3648 +chr20 21200000 25700000 4500000 0 1.4588 +chr20 37100000 41100000 4000000 0 1.3980 +chr20 49200000 54400000 5200000 0 1.5272 +chr21 38600000 41877429 8344323 0 1.6542 +chr21 41878628 43505733 8344323 0 1.9800 +chr21 43507092 46944323 8344323 0 1.6567 +chr22 16300000 18889431 11600000 0 NA +chr22 18939431 27900000 11600000 0 1.6362 +chr22 30500000 39300000 8800000 0 1.7639 +chr22 42600000 43057958 7091432 0 1.5600 +chr22 43107958 47356150 7091432 0 1.5400 +chr22 47366250 48750436 7091432 0 1.4700 +chr22 48767136 49087576 7091432 0 NA +chr22 49089176 49107103 7091432 0 NA +chr22 49126803 49591432 7091432 0 NA +chr3 14700000 23800000 9100000 0 1.3438 +chr3 26400000 32100000 5700000 0 1.2751 +chr3 54400000 58500000 4100000 0 1.4670 +chr3 63700000 66115833 8100000 0 1.1571 +chr3 66375833 71800000 8100000 0 1.4022 +chr3 74200000 87200000 13000000 0 1.2224 +chr3 104400000 123400000 19000000 0 1.4382 +chr3 144400000 150400000 6000000 0 1.2844 +chr4 0 1413146 5200000 0 NA +chr4 1464146 3883456 5200000 0 1.8696 +chr4 3963456 5200000 5200000 0 1.4193 +chr4 40900000 45600000 4700000 0 1.3101 +chr4 59200000 59429326 17300000 0 NA +chr4 59479326 69275441 17300000 0 1.1989 +chr4 69375441 71536854 17300000 0 1.4800 +chr4 71711854 75641303 17300000 0 1.3813 +chr4 75671303 76500000 17300000 0 NA +chr4 102500000 120600000 18100000 0 1.2796 +chr4 124000000 151000000 27000000 0 1.3180 +chr5 42400000 45800000 3400000 0 1.6600 +chr5 135400000 138812257 11800000 0 1.6903 +chr5 138817257 147200000 11800000 0 1.4559 +chr6 5000 4100000 4100000 0 1.4960 +chr6 7000000 9199728 6500000 0 1.3742 +chr6 9249728 13500000 6500000 0 1.5531 +chr6 15500000 23500000 8000000 0 1.3678 +chr6 26100000 40600000 14500000 0 1.5998 +chr6 45200000 57200000 12000000 0 1.3561 +chr6 63500000 70000000 6500000 0 1.2305 +chr6 75900000 87500000 11600000 0 1.3037 +chr6 92100000 95737264 7800000 0 1.1620 +chr6 95937264 99900000 7800000 0 1.2897 +chr6 139100000 149100000 10000000 0 1.4172 +chr6 164400000 167766922 6499992 0 1.2400 +chr6 167784922 170021897 6499992 0 1.5383 +chr6 170171897 170896992 6499992 0 1.4900 +chr7 19500000 35600000 16100000 0 1.5139 +chr7 37500000 43300000 5800000 0 1.4283 +chr7 46600000 48167949 7300000 0 1.3900 +chr7 48207949 50338125 7300000 0 1.4468 +chr7 50378125 53900000 7300000 0 1.3027 +chr7 61100000 61314455 36800000 0 NA +chr7 61364455 61554592 36800000 0 NA +chr7 61604592 74353660 36800000 0 1.4805 +chr7 74603660 97900000 36800000 0 1.3240 +chr7 117200000 130100000 12900000 0 1.3323 +chr7 132400000 137300000 4900000 0 1.4019 +chr7 142800000 147500000 4700000 0 1.1759 +chr8 12700000 19100000 6400000 0 1.1690 +chr8 29700000 29732668 8800000 0 NA +chr8 29798768 38500000 8800000 0 1.4250 +chr8 55600000 66100000 10500000 0 1.5225 +chr8 74000000 86069241 13200000 0 1.3023 +chr8 86193341 86763703 13200000 0 NA +chr8 86851003 87200000 13200000 0 NA +chr8 99100000 101600000 2500000 0 1.6452 +chr8 106100000 117700000 11600000 0 1.2240 +chr8 127300000 140000000 12700000 0 1.4038 +chr9 9000000 14100000 5100000 0 1.2167 +chr9 40200000 40223029 6500000 0 NA +chr9 40273029 40415834 6500000 0 NA +chr9 40465834 40930341 6500000 0 NA +chr9 40980341 41133214 6500000 0 NA +chr9 41183214 41355793 6500000 0 NA +chr9 41405793 42603951 6500000 0 NA +chr9 42653951 43203694 6500000 0 NA +chr9 43253694 43886565 6500000 0 NA +chr9 43936565 44616642 6500000 0 NA +chr9 44666642 44848289 6500000 0 NA +chr9 44898289 45190199 6500000 0 NA +chr9 45240199 45705517 6500000 0 NA +chr9 45755517 46106426 6500000 0 NA +chr9 46156426 46351035 6500000 0 NA +chr9 46401035 46700000 6500000 0 NA +chr9 91000000 91533236 10600000 0 NA +chr9 91583236 91668616 10600000 0 NA +chr9 91718616 101600000 10600000 0 1.5580 +chrX 0 34821 6000000 0 NA +chrX 84821 171384 6000000 0 NA +chrX 201384 967557 6000000 0 NA +chrX 1017557 1054113 6000000 0 NA +chrX 1104113 1184234 6000000 0 NA +chrX 1274234 2028238 6000000 0 NA +chrX 2128238 6000000 6000000 0 1.4520 +chrX 9500000 37008177 37800000 0 1.4047 +chrX 37033177 47300000 37800000 0 1.4817 +chrX 65000000 67700000 2700000 0 1.2658 +chrX 76000000 76570348 22200000 0 NA +chrX 76590348 98200000 22200000 0 1.2642 +chrX 102500000 113403924 18200000 0 1.3809 +chrX 113473924 115596318 18200000 0 1.3900 +chrX 115616318 120700000 18200000 0 1.4886 +chrX 137800000 140100000 2300000 0 1.2054 +chrX 141900000 143335010 5000000 0 1.2700 +chrX 143365010 146900000 5000000 0 1.2517 +chrY 0 34821 11200000 0 NA +chrY 84821 171384 11200000 0 NA +chrY 201384 967557 11200000 0 NA +chrY 1017557 1054113 11200000 0 NA +chrY 1104113 1184234 11200000 0 NA +chrY 1274234 2028238 11200000 0 NA +chrY 2128238 8974955 11200000 0 1.3987 +chrY 9024955 9301322 11200000 0 NA +chrY 9901322 10714553 11200000 0 1.2200 +chrY 12500000 22310816 14700000 0 1.3632 +chrY 22360816 27200000 14700000 0 1.2357 \ No newline at end of file diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 test-data/interval_interpolate.bed --- /dev/null +++ b/test-data/interval_interpolate.bed @@ -0,0 +1,219 @@ +chr1 27800000 29750669 6600000 0 +chr1 29800669 34400000 6600000 0 +chr1 46500000 51300000 4800000 0 +chr1 102000000 103636494 5000000 0 +chr1 103686494 103843838 5000000 0 +chr1 103893838 107000000 5000000 0 +chr1 117600000 120498679 3100000 0 +chr1 120548679 120700000 3100000 0 +chr1 153300000 154800000 1500000 0 +chr1 163800000 171200000 7400000 0 +chr1 174300000 184000000 9700000 0 +chr10 6700000 12300000 5600000 0 +chr10 42100000 45746970 11200000 0 +chr10 45896970 46849175 11200000 0 +chr10 46999175 47262482 11200000 0 +chr10 47412482 47575713 11200000 0 +chr10 47725713 48715542 11200000 0 +chr10 48865542 50807416 11200000 0 +chr10 50857416 51068851 11200000 0 +chr10 51118851 53300000 11200000 0 +chr10 74600000 81231464 15000000 0 +chr10 81241464 89600000 15000000 0 +chr10 98000000 99400000 1400000 0 +chr10 102000000 111800000 9800000 0 +chr10 114900000 119100000 4200000 0 +chr10 127400000 128606059 7974737 0 +chr10 128656059 133271394 7974737 0 +chr10 133281394 133527517 7974737 0 +chr10 133577517 135374737 7974737 0 +chr11 50000 1152759 2800000 0 +chr11 1169335 2800000 2800000 0 +chr11 12600000 16100000 3500000 0 +chr11 21600000 26000000 4400000 0 +chr11 27200000 31000000 3800000 0 +chr11 48800000 50740429 2600000 0 +chr11 50947429 51400000 2600000 0 +chr11 56400000 68846377 12800000 0 +chr11 68848982 69200000 12800000 0 +chr11 70700000 85300000 14600000 0 +chr11 87900000 92300000 4400000 0 +chr11 110000000 115400000 5400000 0 +chr11 120700000 134452384 13752384 0 +chr12 36500000 44600000 8100000 0 +chr12 53100000 56300000 3200000 0 +chr12 69800000 74791940 8900000 0 +chr12 75041940 78700000 8900000 0 +chr12 91200000 94800000 3600000 0 +chr12 100000000 107500000 7500000 0 +chr13 18400000 32900000 14500000 0 +chr13 34700000 39500000 4800000 0 +chr13 52200000 57600000 5400000 0 +chr13 60500000 72100000 11600000 0 +chr13 100500000 109100000 8600000 0 +chr14 23600000 31800000 8200000 0 +chr14 36900000 41000000 4100000 0 +chr14 43200000 48300000 5100000 0 +chr14 69300000 72900000 3600000 0 +chr14 78400000 88900000 10500000 0 +chr15 23300000 24687593 2400000 0 +chr15 24731593 25169606 2400000 0 +chr15 25270606 25700000 2400000 0 +chr15 31400000 37900000 6500000 0 +chr15 65300000 70400000 5100000 0 +chr16 0 8576922 14700000 0 +chr16 8594422 14700000 14700000 0 +chr16 16700000 21700000 5000000 0 +chr16 69400000 78200000 8800000 0 +chr16 80500000 82700000 2200000 0 +chr17 0 296854 11200000 0 +chr17 343376 11200000 11200000 0 +chr17 28800000 31699961 6600000 0 +chr17 31799961 33429230 6600000 0 +chr17 33529230 35400000 6600000 0 +chr17 47600000 54900000 7300000 0 +chr17 55600000 59900000 4300000 0 +chr18 17300000 31000000 13700000 0 +chr18 35500000 50313134 16500000 0 +chr18 50360134 52000000 16500000 0 +chr19 11000 7297004 12600000 0 +chr19 7302004 8593198 12600000 0 +chr19 8598198 12600000 12600000 0 +chr2 12800000 16171403 4200000 0 +chr2 16221403 17000000 4200000 0 +chr2 19100000 21012571 4800000 0 +chr2 21037571 23900000 4800000 0 +chr2 47600000 52700000 5100000 0 +chr2 54800000 61100000 6300000 0 +chr2 75400000 83700000 8300000 0 +chr2 102100000 108600000 6500000 0 +chr2 113800000 118600000 4800000 0 +chr2 129600000 134800000 5200000 0 +chr2 136600000 144700000 8100000 0 +chr2 189100000 197100000 8000000 0 +chr2 209100000 215100000 6000000 0 +chr2 221300000 233711985 15700000 0 +chr2 233731985 237000000 15700000 0 +chr20 5000000 9000000 4000000 0 +chr20 11900000 17800000 5900000 0 +chr20 21200000 25700000 4500000 0 +chr20 37100000 41100000 4000000 0 +chr20 49200000 54400000 5200000 0 +chr21 38600000 41877429 8344323 0 +chr21 41878628 43505733 8344323 0 +chr21 43507092 46944323 8344323 0 +chr22 16300000 18889431 11600000 0 +chr22 18939431 27900000 11600000 0 +chr22 30500000 39300000 8800000 0 +chr22 42600000 43057958 7091432 0 +chr22 43107958 47356150 7091432 0 +chr22 47366250 48750436 7091432 0 +chr22 48767136 49087576 7091432 0 +chr22 49089176 49107103 7091432 0 +chr22 49126803 49591432 7091432 0 +chr3 14700000 23800000 9100000 0 +chr3 26400000 32100000 5700000 0 +chr3 54400000 58500000 4100000 0 +chr3 63700000 66115833 8100000 0 +chr3 66375833 71800000 8100000 0 +chr3 74200000 87200000 13000000 0 +chr3 104400000 123400000 19000000 0 +chr3 144400000 150400000 6000000 0 +chr4 0 1413146 5200000 0 +chr4 1464146 3883456 5200000 0 +chr4 3963456 5200000 5200000 0 +chr4 40900000 45600000 4700000 0 +chr4 59200000 59429326 17300000 0 +chr4 59479326 69275441 17300000 0 +chr4 69375441 71536854 17300000 0 +chr4 71711854 75641303 17300000 0 +chr4 75671303 76500000 17300000 0 +chr4 102500000 120600000 18100000 0 +chr4 124000000 151000000 27000000 0 +chr5 42400000 45800000 3400000 0 +chr5 135400000 138812257 11800000 0 +chr5 138817257 147200000 11800000 0 +chr6 5000 4100000 4100000 0 +chr6 7000000 9199728 6500000 0 +chr6 9249728 13500000 6500000 0 +chr6 15500000 23500000 8000000 0 +chr6 26100000 40600000 14500000 0 +chr6 45200000 57200000 12000000 0 +chr6 63500000 70000000 6500000 0 +chr6 75900000 87500000 11600000 0 +chr6 92100000 95737264 7800000 0 +chr6 95937264 99900000 7800000 0 +chr6 139100000 149100000 10000000 0 +chr6 164400000 167766922 6499992 0 +chr6 167784922 170021897 6499992 0 +chr6 170171897 170896992 6499992 0 +chr7 19500000 35600000 16100000 0 +chr7 37500000 43300000 5800000 0 +chr7 46600000 48167949 7300000 0 +chr7 48207949 50338125 7300000 0 +chr7 50378125 53900000 7300000 0 +chr7 61100000 61314455 36800000 0 +chr7 61364455 61554592 36800000 0 +chr7 61604592 74353660 36800000 0 +chr7 74603660 97900000 36800000 0 +chr7 117200000 130100000 12900000 0 +chr7 132400000 137300000 4900000 0 +chr7 142800000 147500000 4700000 0 +chr8 12700000 19100000 6400000 0 +chr8 29700000 29732668 8800000 0 +chr8 29798768 38500000 8800000 0 +chr8 55600000 66100000 10500000 0 +chr8 74000000 86069241 13200000 0 +chr8 86193341 86763703 13200000 0 +chr8 86851003 87200000 13200000 0 +chr8 99100000 101600000 2500000 0 +chr8 106100000 117700000 11600000 0 +chr8 127300000 140000000 12700000 0 +chr9 9000000 14100000 5100000 0 +chr9 40200000 40223029 6500000 0 +chr9 40273029 40415834 6500000 0 +chr9 40465834 40930341 6500000 0 +chr9 40980341 41133214 6500000 0 +chr9 41183214 41355793 6500000 0 +chr9 41405793 42603951 6500000 0 +chr9 42653951 43203694 6500000 0 +chr9 43253694 43886565 6500000 0 +chr9 43936565 44616642 6500000 0 +chr9 44666642 44848289 6500000 0 +chr9 44898289 45190199 6500000 0 +chr9 45240199 45705517 6500000 0 +chr9 45755517 46106426 6500000 0 +chr9 46156426 46351035 6500000 0 +chr9 46401035 46700000 6500000 0 +chr9 91000000 91533236 10600000 0 +chr9 91583236 91668616 10600000 0 +chr9 91718616 101600000 10600000 0 +chrX 0 34821 6000000 0 +chrX 84821 171384 6000000 0 +chrX 201384 967557 6000000 0 +chrX 1017557 1054113 6000000 0 +chrX 1104113 1184234 6000000 0 +chrX 1274234 2028238 6000000 0 +chrX 2128238 6000000 6000000 0 +chrX 9500000 37008177 37800000 0 +chrX 37033177 47300000 37800000 0 +chrX 65000000 67700000 2700000 0 +chrX 76000000 76570348 22200000 0 +chrX 76590348 98200000 22200000 0 +chrX 102500000 113403924 18200000 0 +chrX 113473924 115596318 18200000 0 +chrX 115616318 120700000 18200000 0 +chrX 137800000 140100000 2300000 0 +chrX 141900000 143335010 5000000 0 +chrX 143365010 146900000 5000000 0 +chrY 0 34821 11200000 0 +chrY 84821 171384 11200000 0 +chrY 201384 967557 11200000 0 +chrY 1017557 1054113 11200000 0 +chrY 1104113 1184234 11200000 0 +chrY 1274234 2028238 11200000 0 +chrY 2128238 8974955 11200000 0 +chrY 9024955 9301322 11200000 0 +chrY 9901322 10714553 11200000 0 +chrY 12500000 22310816 14700000 0 +chrY 22360816 27200000 14700000 0 \ No newline at end of file diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 test-data/logreg_inp.tabular --- /dev/null +++ b/test-data/logreg_inp.tabular @@ -0,0 +1,100 @@ +2.04 2.01 1070 1 +2.56 3.4 1254 1 +3.75 3.68 1466 1 +1.1 1.54 706 1 +3 3.32 1160 1 +0.05 0.33 756 1 +1.38 0.36 1058 1 +1.5 1.97 1008 0 +1.38 2.03 1104 1 +4.01 2.05 1200 1 +1.5 2.13 896 1 +1.29 1.34 848 1 +1.9 1.51 958 1 +3.11 3.12 1246 0 +1.92 2.14 1106 1 +0.81 2.6 790 1 +1.01 1.9 954 1 +3.66 3.06 1500 0 +2 1.6 1046 0 +2.05 1.96 1054 1 +2.6 1.96 1198 0 +2.55 1.56 940 1 +0.38 1.6 456 0 +2.48 1.92 1150 1 +2.74 3.09 636 0 +1.77 0.78 744 1 +1.61 2.12 644 0 +0.99 1.85 842 1 +1.62 1.78 852 1 +2.03 1.03 1170 0 +3.5 3.44 1034 1 +3.18 2.42 1202 1 +2.39 1.74 1018 1 +1.48 1.89 1180 1 +1.54 1.43 952 0 +1.57 1.64 1038 1 +2.46 2.69 1090 0 +2.42 1.79 694 0 +2.11 2.72 1096 0 +2.04 2.15 1114 0 +1.68 2.22 1256 1 +1.64 1.55 1208 0 +2.41 2.34 820 0 +2.1 2.92 1222 0 +1.4 2.1 1120 0 +2.03 1.64 886 0 +1.99 2.83 1126 0 +2.24 1.76 1158 0 +0.45 1.81 676 0 +2.31 2.68 1214 0 +2.41 2.55 1136 1 +2.56 2.7 1264 0 +2.5 1.66 1116 1 +2.92 2.23 1292 1 +2.35 2.01 604 1 +2.82 1.24 854 1 +1.8 1.95 814 0 +1.29 1.73 778 1 +1.68 1.08 800 0 +3.44 3.46 1424 0 +1.9 3.01 950 0 +2.06 0.54 1056 1 +3.3 3.2 956 1 +1.8 1.5 1352 1 +2 1.71 852 1 +1.68 1.99 1168 0 +1.94 2.76 970 1 +0.97 1.56 776 1 +1.12 1.78 854 1 +1.31 1.32 1232 0 +1.68 0.87 1140 0 +3.09 1.75 1084 0 +1.87 1.41 954 0 +2 2.77 1000 0 +2.39 1.78 1084 0 +1.5 1.34 1058 0 +1.82 1.52 816 0 +1.8 2.97 1146 0 +2.01 1.75 1000 1 +1.88 1.64 856 1 +1.64 1.8 798 1 +2.42 3.37 1324 1 +0.22 1.15 704 1 +2.31 1.72 1222 1 +0.95 2.27 948 0 +1.99 2.85 1182 0 +1.86 2.21 1000 1 +1.79 1.94 910 0 +3.02 4.25 1374 1 +1.85 1.83 1014 0 +1.98 2.75 1420 0 +2.15 1.71 400 0 +1.46 2.2 998 1 +2.29 2.13 776 1 +2.39 2.38 1134 0 +1.8 1.64 772 0 +2.64 1.87 1304 0 +2.08 2.53 1212 0 +0.7 1.78 818 1 +0.89 1.2 864 1 \ No newline at end of file diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 test-data/logreg_out2.tabular --- /dev/null +++ b/test-data/logreg_out2.tabular @@ -0,0 +1,19 @@ +response column c4 +predictor column(s) c1,c2,c3 +Y-intercept 0.9111624714 +p-value (Y-intercept) 0.3571052008 +Slope (c1) 0.057995684 +p-value (c1) 0.8677866885 +Slope (c2) -0.2319990287 +p-value (c2) 0.4986584837 +Slope (c3) -0.0004633556 +p-value (c3) 0.6785709433 +Null deviance 138.46939 +Residual deviance 137.44023 +pseudo R-squared 0.00743 + + +vif +c1 1.65649272465 +c2 1.47696547452 +c3 1.4307725027 diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 test-data/partialR_result.tabular --- /dev/null +++ b/test-data/partialR_result.tabular @@ -0,0 +1,4 @@ +#Model R-sq partial_R_Terms partial_R_Value +1 2 0.9388 - - +2 0.7280 1 0.7750 +1 0.9104 2 0.3167 \ No newline at end of file diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 test-data/value_interpolate.bed --- /dev/null +++ b/test-data/value_interpolate.bed @@ -0,0 +1,2794 @@ +chr9 138121811 138276705 0 +chr10 80472830 80642405 0 +chr10 112350437 112499683 0 +chr18 44397436 44552893 0 +chr19 40876527 41059171 0 +chr19 47423994 47633625 0 +chr4 189654626 189851415 0.98 +chr14 41341385 41500259 0.99 +chr10 68138624 68290770 1.02 +chr1 194486128 194637476 1.03 +chr2 79730231 79730948 1.04 +chr8 77233600 77394640 1.04 +chr3 164366661 164516629 1.05 +chr4 182025739 182177673 1.05 +chr7 78310705 78485567 1.05 +chr8 34287764 34455012 1.05 +chr9 28840514 28840768 1.05 +chr16 62921688 63071731 1.05 +chr18 47940868 48119474 1.05 +chr1 215306203 215473376 1.06 +chr3 77584830 77738883 1.06 +chr5 165677723 165841050 1.06 +chr7 79539152 79645155 1.06 +chr8 4596118 4755769 1.06 +chr9 105583839 105754361 1.06 +chr10 56709834 56882157 1.06 +chr14 25538290 25757198 1.06 +chr1 213686097 213890738 1.07 +chr4 93630873 93775737 1.07 +chr11 37673721 37838885 1.07 +chr13 34909940 35080870 1.07 +chrX 66474376 66475056 1.07 +chr2 116313042 116516750 1.08 +chr5 28750780 28917994 1.08 +chr7 83471581 83635351 1.08 +chr9 24090721 24243438 1.08 +chr13 66420966 66421707 1.08 +chr14 83853577 84001147 1.08 +chr18 67890358 68041609 1.08 +chr4 45608190 45780635 1.09 +chr6 93072999 93225231 1.09 +chr8 16444047 16618550 1.09 +chr8 16444054 16618572 1.09 +chr8 89218072 89400932 1.09 +chr13 86846880 86847238 1.09 +chr13 88476238 88642707 1.09 +chr21 21123964 21277726 1.09 +chr1 162445439 162629394 1.1 +chr3 19336535 19494336 1.1 +chr3 20636022 20828571 1.1 +chr4 28748334 28909642 1.1 +chr4 65766258 65968108 1.1 +chr5 25187412 25333169 1.1 +chr13 82592804 82779046 1.1 +chr1 72790142 72983871 1.11 +chr1 193197934 193377693 1.11 +chr2 165238485 165239216 1.11 +chr3 65182404 65361116 1.11 +chr4 31129762 31286231 1.11 +chr4 160913220 161055170 1.11 +chr6 68213170 68355088 1.11 +chr6 121373179 121565240 1.11 +chr7 14108578 14279180 1.11 +chr7 83230804 83414012 1.11 +chr7 145959694 146102483 1.11 +chr8 85136758 85319037 1.11 +chr9 121031753 121032407 1.11 +chr12 80260958 80413762 1.11 +chr13 84781269 84965380 1.11 +chr13 93229523 93403152 1.11 +chr18 56874822 57034619 1.11 +chr2 78739654 78929106 1.12 +chr4 28258977 28448528 1.12 +chr6 81474419 81671778 1.12 +chr8 5106183 5256062 1.12 +chr8 15557561 15699688 1.12 +chr9 120913556 121095236 1.12 +chr11 41835501 42018028 1.12 +chr11 54839338 55007563 1.12 +chr13 82035437 82201311 1.12 +chr14 29256331 29403102 1.12 +chr14 79340662 79520755 1.12 +chr15 52652616 52815257 1.12 +chrX 139430017 139573923 1.12 +chr2 22411649 22627999 1.13 +chr2 146591634 146775911 1.13 +chr2 165687734 165861526 1.13 +chr3 21574294 21742228 1.13 +chr3 84623846 84624672 1.13 +chr4 118641782 118817785 1.13 +chr4 121715867 121905993 1.13 +chr4 171952201 171952789 1.13 +chr5 113199522 113366010 1.13 +chr5 121061430 121226477 1.13 +chr6 19131250 19131847 1.13 +chr6 62586405 62725287 1.13 +chr6 96322030 96503179 1.13 +chr7 10359917 10517208 1.13 +chr8 88145759 88302638 1.13 +chr9 119255125 119430632 1.13 +chr12 72112997 72312988 1.13 +chr14 28380477 28530538 1.13 +chr14 43288620 43453116 1.13 +chr14 84168823 84169632 1.13 +chr15 45134914 45309877 1.13 +chr16 60359424 60359773 1.13 +chr1 73625862 73764288 1.14 +chr1 73625895 73764288 1.14 +chr1 103893882 104012508 1.14 +chr1 215209972 215373859 1.14 +chr2 50386937 50387456 1.14 +chr2 81618970 81763529 1.14 +chr3 74234096 74391234 1.14 +chr3 80133749 80297117 1.14 +chr3 95519627 95706411 1.14 +chr3 98524261 98714848 1.14 +chr4 117579710 117743718 1.14 +chr4 171958495 172124798 1.14 +chr6 153105001 153217269 1.14 +chr7 52666632 52867234 1.14 +chr7 151558264 151732070 1.14 +chr8 112515727 112683952 1.14 +chr9 120076427 120271081 1.14 +chr11 90725899 90902355 1.14 +chr12 33170553 33333489 1.14 +chr13 106527433 106708235 1.14 +chr14 82169838 82169979 1.14 +chr18 26696101 26878287 1.14 +chr21 20389919 20390568 1.14 +chrX 92132730 92310405 1.14 +chrX 127763208 127921922 1.14 +chrX 137822409 137909028 1.14 +chr1 185443401 185614023 1.15 +chr2 49302567 49513831 1.15 +chr2 66394115 66576561 1.15 +chr2 83447798 83650748 1.15 +chr2 83447850 83650731 1.15 +chr2 115793564 115936236 1.15 +chr3 1347065 1517440 1.15 +chr3 119242817 119243578 1.15 +chr4 19599685 19772537 1.15 +chr4 112315831 112504400 1.15 +chr5 4351733 4533614 1.15 +chr5 52909243 53065314 1.15 +chr5 119668835 119831950 1.15 +chr5 120133844 120134605 1.15 +chr6 58610235 58788598 1.15 +chr6 65092079 65280053 1.15 +chr6 127363622 127546913 1.15 +chr7 84729683 84904466 1.15 +chr8 112100512 112266163 1.15 +chr8 132628905 132799581 1.15 +chr11 29071526 29253144 1.15 +chr11 90903367 91076093 1.15 +chr12 76524989 76660180 1.15 +chr12 83887475 84059958 1.15 +chr13 37217921 37218551 1.15 +chr13 70296348 70477659 1.15 +chr13 71078490 71256801 1.15 +chr13 89314358 89462845 1.15 +chr16 63462907 63616441 1.15 +chr18 37305512 37482566 1.15 +chr1 61044529 61191204 1.16 +chr1 74413167 74572629 1.16 +chr1 175137086 175272030 1.16 +chr1 188252333 188442427 1.16 +chr1 195614240 195777468 1.16 +chr2 52903987 53101240 1.16 +chr2 122950238 123095058 1.16 +chr2 193770382 193940783 1.16 +chr3 145885187 146013616 1.16 +chr4 60813762 60999837 1.16 +chr4 67050357 67201898 1.16 +chr4 91753784 91923803 1.16 +chr4 93320220 93531288 1.16 +chr4 107736688 107943788 1.16 +chr4 179406134 179587711 1.16 +chr5 16193697 16358412 1.16 +chr6 8464242 8664833 1.16 +chr6 22580422 22736310 1.16 +chr6 100781592 100953408 1.16 +chr7 51595765 51731468 1.16 +chr7 121314792 121472652 1.16 +chr7 123711877 123721155 1.16 +chr8 40501073 40649244 1.16 +chr8 94214334 94424255 1.16 +chr9 1340595 1488472 1.16 +chr9 82142727 82324329 1.16 +chr9 84622895 84783002 1.16 +chr11 91893716 92071858 1.16 +chr13 52804333 52973849 1.16 +chr13 90342056 90524017 1.16 +chr14 32378912 32540754 1.16 +chr14 37010397 37213971 1.16 +chr14 86077200 86077996 1.16 +chr15 95508740 95699545 1.16 +chr17 13211602 13416771 1.16 +chr17 48565062 48664511 1.16 +chr17 49013829 49212688 1.16 +chr21 13462479 13626951 1.16 +chrX 82482413 82601150 1.16 +chrX 84848634 85055360 1.16 +chrX 143420958 143549753 1.16 +chrY 4917122 4917768 1.16 +chrY 22545464 22688210 1.16 +chr1 71927607 72094826 1.17 +chr1 188974250 189126140 1.17 +chr2 50398943 50560714 1.17 +chr2 153550169 153751686 1.17 +chr3 167343076 167517669 1.17 +chr3 189731929 189918719 1.17 +chr4 12410352 12586679 1.17 +chr4 18028014 18213766 1.17 +chr4 22911964 23061775 1.17 +chr4 44529019 44675676 1.17 +chr4 59533297 59706381 1.17 +chr4 63770893 63973350 1.17 +chr4 101813229 101963326 1.17 +chr4 126354725 126519739 1.17 +chr4 135464067 135640784 1.17 +chr4 188522259 188676888 1.17 +chr5 135915628 135916057 1.17 +chr6 78821469 78822201 1.17 +chr6 113213629 113374404 1.17 +chr7 10241793 10242601 1.17 +chr7 108884723 109038822 1.17 +chr8 113491855 113659908 1.17 +chr10 67500726 67674349 1.17 +chr11 97102811 97328868 1.17 +chr11 100145397 100299538 1.17 +chr11 126420648 126568009 1.17 +chr11 132215897 132391418 1.17 +chr13 35921756 36083776 1.17 +chr13 37797297 37951253 1.17 +chr13 60719187 60882412 1.17 +chr13 65294071 65461528 1.17 +chr16 61098268 61286985 1.17 +chr17 66194639 66343264 1.17 +chr18 35489814 35657078 1.17 +chrX 66577693 66716608 1.17 +chrX 140211864 140377242 1.17 +chr1 75234702 75388843 1.18 +chr1 102366355 102506882 1.18 +chr1 106102412 106261344 1.18 +chr1 118630891 118728136 1.18 +chr1 192296623 192424789 1.18 +chr2 2210217 2390008 1.18 +chr2 58959370 59141335 1.18 +chr2 72247356 72428550 1.18 +chr2 82554305 82716433 1.18 +chr2 164085637 164247934 1.18 +chr2 199939521 200121814 1.18 +chr3 29867269 30018470 1.18 +chr3 73854905 74025974 1.18 +chr3 80611850 80795918 1.18 +chr4 67263328 67431970 1.18 +chr5 20414933 20415717 1.18 +chr5 26353486 26539834 1.18 +chr5 109349873 109506373 1.18 +chr5 144610138 144777056 1.18 +chr5 165737022 165925395 1.18 +chr6 50775747 50907356 1.18 +chr6 77624011 77788328 1.18 +chr7 81437179 81635132 1.18 +chr7 114921735 115101378 1.18 +chr7 144927075 145051342 1.18 +chr8 3899657 4066797 1.18 +chr8 16323134 16472679 1.18 +chr8 18447618 18644291 1.18 +chr8 36363082 36516620 1.18 +chr8 78915842 79062005 1.18 +chr8 92382762 92553372 1.18 +chr8 117283875 117484166 1.18 +chr9 12225005 12396287 1.18 +chr9 13729630 13912935 1.18 +chr9 23376562 23557443 1.18 +chr11 15778995 15940672 1.18 +chr11 26201318 26402942 1.18 +chr12 58234674 58235464 1.18 +chr12 60211160 60370325 1.18 +chr13 58618179 58784565 1.18 +chr14 41996958 42171562 1.18 +chr16 61943900 62080675 1.18 +chr18 25599572 25763880 1.18 +chr18 33450045 33623547 1.18 +chr20 12493232 12650585 1.18 +chrX 98334220 98501957 1.18 +chr1 98681325 98839806 1.19 +chr1 106759920 106884446 1.19 +chr1 188892773 189058601 1.19 +chr2 14523373 14693180 1.19 +chr2 40346711 40526003 1.19 +chr2 52193265 52373212 1.19 +chr2 167445852 167644200 1.19 +chr2 192914919 193076744 1.19 +chr2 211871510 212039574 1.19 +chr2 221335235 221512023 1.19 +chr3 2181314 2367266 1.19 +chr3 26718031 26882772 1.19 +chr3 149511877 149686925 1.19 +chr3 166908960 167080203 1.19 +chr3 169544516 169732146 1.19 +chr4 73321382 73500901 1.19 +chr4 103999791 104161865 1.19 +chr4 121155218 121304760 1.19 +chr4 132769467 132770056 1.19 +chr5 12130086 12282083 1.19 +chr5 163117739 163118288 1.19 +chr6 93629580 93823810 1.19 +chr7 12919070 13079613 1.19 +chr8 15279388 15445713 1.19 +chr8 118297084 118467725 1.19 +chr8 136470371 136653380 1.19 +chr8 137773461 137919469 1.19 +chr9 8398615 8557610 1.19 +chr9 12819715 13004078 1.19 +chr10 106748194 106912778 1.19 +chr11 80418056 80571094 1.19 +chr12 43511484 43650184 1.19 +chr12 76975277 77145211 1.19 +chr12 81988506 82172342 1.19 +chr14 26561378 26744993 1.19 +chr14 32835734 32999496 1.19 +chr14 40363695 40563204 1.19 +chr14 61886249 62046831 1.19 +chr16 6165969 6339847 1.19 +chr18 1551143 1734122 1.19 +chr18 47274877 47443274 1.19 +chrX 28257745 28427558 1.19 +chrX 50147722 50314773 1.19 +chrX 83845334 83846061 1.19 +chrX 141360475 141476781 1.19 +chr1 76640828 76740449 1.2 +chr1 95507969 95685941 1.2 +chr2 5764246 5956172 1.2 +chr2 63195188 63358079 1.2 +chr2 76357535 76358153 1.2 +chr2 77030009 77195797 1.2 +chr2 103816337 103983810 1.2 +chr2 150592885 150779567 1.2 +chr2 185140095 185334504 1.2 +chr2 224018074 224018802 1.2 +chr3 35590216 35766097 1.2 +chr3 105180963 105347887 1.2 +chr3 116940533 117110558 1.2 +chr4 60219574 60371862 1.2 +chr4 64867965 65032881 1.2 +chr4 179834815 180006495 1.2 +chr5 98954115 99111442 1.2 +chr5 100991121 101181889 1.2 +chr5 102992405 103190088 1.2 +chr6 69438735 69596850 1.2 +chr6 72157311 72342266 1.2 +chr6 102119094 102294630 1.2 +chr7 136263965 136416914 1.2 +chr8 4187855 4386741 1.2 +chr8 24278979 24469357 1.2 +chr8 75850345 76026015 1.2 +chr9 28027075 28204449 1.2 +chr9 75584378 75757095 1.2 +chr9 80495074 80676552 1.2 +chr9 103481717 103640785 1.2 +chr9 118368426 118523858 1.2 +chr10 10034833 10202423 1.2 +chr10 16038936 16196794 1.2 +chr10 57554233 57708605 1.2 +chr11 79636827 79817896 1.2 +chr12 18499860 18558864 1.2 +chr12 71742532 71963509 1.2 +chr12 86171722 86361997 1.2 +chr13 61999015 62167070 1.2 +chr13 86241575 86431660 1.2 +chr13 87858551 88035241 1.2 +chr15 24370637 24550028 1.2 +chr18 14843504 14998636 1.2 +chr18 40044271 40244003 1.2 +chr18 73605047 73776512 1.2 +chr21 24509084 24509572 1.2 +chrX 27748283 27942635 1.2 +chrX 140718912 140719477 1.2 +chrY 24157833 24324267 1.2 +chrY 26443321 26640420 1.2 +chrY 26443368 26443796 1.2 +chr1 81828369 81949209 1.21 +chr2 4779724 4937997 1.21 +chr2 141388838 141564023 1.21 +chr2 229316646 229475755 1.21 +chr2 235485097 235485527 1.21 +chr3 22372670 22514731 1.21 +chr3 34981231 35131137 1.21 +chr3 64198721 64357750 1.21 +chr3 86088438 86268550 1.21 +chr4 29381036 29558458 1.21 +chr4 96994009 97153613 1.21 +chr4 96994043 97153591 1.21 +chr4 98832394 98983885 1.21 +chr4 149764955 149967710 1.21 +chr4 158576067 158755508 1.21 +chr4 189957368 190151886 1.21 +chr5 166550028 166647640 1.21 +chr6 23499391 23684090 1.21 +chr6 66876814 67021555 1.21 +chr6 85028192 85201298 1.21 +chr6 94783870 94898446 1.21 +chr7 85371156 85485924 1.21 +chr7 118629320 118806601 1.21 +chr8 122743719 122893632 1.21 +chr9 11170428 11341967 1.21 +chr9 17839221 18013839 1.21 +chr9 82951747 83098620 1.21 +chr9 105214273 105379246 1.21 +chr9 117911982 118081322 1.21 +chr10 24159166 24341427 1.21 +chr10 36972913 37144238 1.21 +chr10 53912343 54101552 1.21 +chr11 104950107 105093266 1.21 +chr12 28228135 28228812 1.21 +chr12 58917391 59098004 1.21 +chr12 65509085 65693926 1.21 +chr12 87561517 87715880 1.21 +chr12 113937914 114127188 1.21 +chr13 56523630 56717407 1.21 +chr13 64727675 64900032 1.21 +chr13 77789426 77997374 1.21 +chr13 91284445 91452121 1.21 +chr13 92499797 92681312 1.21 +chr13 95333180 95498792 1.21 +chr13 100313432 100480730 1.21 +chr13 105286223 105459675 1.21 +chr14 33897190 34061283 1.21 +chr14 37337239 37487767 1.21 +chr15 41409166 41591912 1.21 +chr17 47844680 48028972 1.21 +chr18 61594910 61752944 1.21 +chr18 68809466 69000766 1.21 +chr20 9197717 9340955 1.21 +chr20 26006126 26220344 1.21 +chr21 40630421 40630982 1.21 +chrX 28029871 28187101 1.21 +chrX 81388806 81389596 1.21 +chrX 143983780 144162387 1.21 +chrX 145778720 145884794 1.21 +chrY 23612040 23794544 1.21 +chr1 102758472 102908844 1.22 +chr1 175974940 176155256 1.22 +chr2 77934103 78088722 1.22 +chr2 155527554 155706737 1.22 +chr2 188855759 189025686 1.22 +chr3 82655862 82799401 1.22 +chr3 87569879 87730255 1.22 +chr3 111576135 111740734 1.22 +chr3 159590747 159779789 1.22 +chr3 168566872 168739931 1.22 +chr3 182297993 182470707 1.22 +chr3 183174949 183353241 1.22 +chr4 11429778 11591553 1.22 +chr4 34462310 34632915 1.22 +chr4 58185990 58349753 1.22 +chr4 96348969 96517904 1.22 +chr4 116748182 116937051 1.22 +chr4 128468219 128623706 1.22 +chr4 142401195 142587354 1.22 +chr4 157665396 157840215 1.22 +chr4 182510849 182656131 1.22 +chr5 13218128 13374982 1.22 +chr5 56921507 57065064 1.22 +chr5 59340905 59341702 1.22 +chr6 48824625 48993352 1.22 +chr6 122952226 123105805 1.22 +chr7 11772911 11940697 1.22 +chr7 110006153 110166536 1.22 +chr7 110006194 110166524 1.22 +chr7 114005319 114175715 1.22 +chr8 31908045 32086561 1.22 +chr8 50882181 51087101 1.22 +chr9 77399680 77400521 1.22 +chr9 117248721 117432842 1.22 +chr10 52898354 53069428 1.22 +chr10 66811696 66961629 1.22 +chr10 107683471 107864697 1.22 +chr10 109290751 109488706 1.22 +chr10 128798554 128968853 1.22 +chr11 20206660 20372501 1.22 +chr11 40928715 41089514 1.22 +chr13 55372320 55544178 1.22 +chr13 57629383 57794495 1.22 +chr13 69138664 69317439 1.22 +chr13 83766576 83924559 1.22 +chr14 27534095 27695886 1.22 +chr15 22736034 22894833 1.22 +chr15 24960345 25130043 1.22 +chr16 6846409 7019846 1.22 +chr16 72300075 72453658 1.22 +chr18 28808364 28955525 1.22 +chr18 63881662 63997076 1.22 +chr21 17336806 17489108 1.22 +chr21 22151924 22333707 1.22 +chr21 30673984 30855176 1.22 +chrX 76640611 76788971 1.22 +chrX 80307294 80471254 1.22 +chrX 90182126 90332983 1.22 +chrX 137286220 137458601 1.22 +chrY 10088705 10250875 1.22 +chr1 56595685 56724303 1.23 +chr1 57456434 57628847 1.23 +chr1 69215584 69364440 1.23 +chr1 105468369 105631855 1.23 +chr1 187292331 187459950 1.23 +chr1 213048945 213231437 1.23 +chr2 56202097 56396479 1.23 +chr2 108502515 108667821 1.23 +chr2 150865067 151039803 1.23 +chr3 27349643 27531253 1.23 +chr3 118153905 118326520 1.23 +chr4 15981564 16185166 1.23 +chr4 18737375 18858211 1.23 +chr4 32138117 32292373 1.23 +chr4 101510584 101672095 1.23 +chr4 109877070 110039180 1.23 +chr4 172208306 172379167 1.23 +chr5 8413255 8600646 1.23 +chr5 90806344 90964362 1.23 +chr5 152539684 152728888 1.23 +chr6 19833231 20018547 1.23 +chr6 57292519 57451664 1.23 +chr6 104541641 104626173 1.23 +chr6 142222399 142361416 1.23 +chr7 67546085 67703913 1.23 +chr7 124677379 124865352 1.23 +chr8 64465972 64651147 1.23 +chr8 69451201 69616969 1.23 +chr8 93160371 93308372 1.23 +chr8 118573281 118737738 1.23 +chr8 139307409 139453227 1.23 +chr9 25690187 25853227 1.23 +chr9 29493271 29639053 1.23 +chr9 108208167 108395870 1.23 +chr10 65699664 65700044 1.23 +chr10 84424513 84577760 1.23 +chr10 130997310 130997876 1.23 +chr11 42864060 43020941 1.23 +chr12 39661185 39833116 1.23 +chr12 126779004 126932042 1.23 +chr13 22108664 22109399 1.23 +chr13 108167626 108341331 1.23 +chr16 58177854 58178675 1.23 +chr16 58776587 58993431 1.23 +chr16 76343711 76539056 1.23 +chr21 23491444 23651727 1.23 +chrX 29011746 29143750 1.23 +chrX 36335465 36531290 1.23 +chrX 45140586 45320235 1.23 +chrX 95213828 95348733 1.23 +chrX 97508539 97716150 1.23 +chrX 128151204 128272268 1.23 +chr1 237591432 237780877 1.24 +chr2 52045737 52233820 1.24 +chr2 129638238 129801244 1.24 +chr2 142692260 142876374 1.24 +chr2 175978460 176145301 1.24 +chr2 186280439 186426350 1.24 +chr3 7560443 7760313 1.24 +chr3 112056227 112057022 1.24 +chr3 191651334 191835494 1.24 +chr4 33712127 33712843 1.24 +chr4 61650910 61798492 1.24 +chr4 150581367 150736943 1.24 +chr5 15191420 15315309 1.24 +chr5 146975536 147089688 1.24 +chr5 174474436 174643725 1.24 +chr6 165061379 165244688 1.24 +chr7 62325611 62464970 1.24 +chr7 62908796 63102874 1.24 +chr7 122549071 122716095 1.24 +chr7 122549101 122549448 1.24 +chr8 14278141 14461142 1.24 +chr8 18249326 18450010 1.24 +chr8 47816618 47982164 1.24 +chr8 51483937 51664779 1.24 +chr8 70223686 70378901 1.24 +chr8 106002673 106200542 1.24 +chr9 9657611 9823754 1.24 +chr9 9932073 10130653 1.24 +chr9 83969058 84168612 1.24 +chr9 110239919 110428437 1.24 +chr10 20561285 20741009 1.24 +chr10 54769832 54938912 1.24 +chr11 21566376 21747063 1.24 +chr11 40185000 40379506 1.24 +chr11 106083757 106252007 1.24 +chr11 132451658 132626345 1.24 +chr12 16742828 16908591 1.24 +chr12 82133389 82134182 1.24 +chr13 96143661 96297462 1.24 +chr15 22736020 22894833 1.24 +chr16 50431379 50432028 1.24 +chr18 7827406 7985296 1.24 +chrX 62268206 62462742 1.24 +chrX 92989831 93172474 1.24 +chrX 104392277 104584586 1.24 +chr1 48917318 49102938 1.25 +chr1 85868659 85959828 1.25 +chr1 184949732 185093551 1.25 +chr2 16245984 16435838 1.25 +chr2 19328431 19511727 1.25 +chr2 19328462 19511695 1.25 +chr2 115568405 115730309 1.25 +chr2 117162771 117350524 1.25 +chr2 166510124 166700135 1.25 +chr2 185796056 186001245 1.25 +chr3 30704315 30705070 1.25 +chr3 68472813 68645976 1.25 +chr3 149208193 149208677 1.25 +chr4 13490350 13655289 1.25 +chr4 20347912 20531936 1.25 +chr4 105288371 105449318 1.25 +chr4 127250479 127434554 1.25 +chr4 127434114 127434494 1.25 +chr5 155562586 155563351 1.25 +chr6 9606972 9733133 1.25 +chr6 82503161 82690511 1.25 +chr6 83405526 83590455 1.25 +chr7 94315045 94480285 1.25 +chr7 95121078 95310527 1.25 +chr8 84378331 84550071 1.25 +chr9 26681234 26681720 1.25 +chr9 81315514 81316298 1.25 +chr10 118689986 118853807 1.25 +chr11 24362912 24551940 1.25 +chr11 86920498 87069201 1.25 +chr11 99613168 99757459 1.25 +chr12 23609734 23784529 1.25 +chr12 30395911 30548839 1.25 +chr13 43148302 43305457 1.25 +chr13 55715887 55883767 1.25 +chr17 12920204 13085706 1.25 +chr18 33631807 33798667 1.25 +chr18 60284571 60448778 1.25 +chr19 23555073 23732302 1.25 +chrX 63351109 63519991 1.25 +chrX 124726353 124825239 1.25 +chrY 8634089 8808181 1.25 +chr1 68609865 68712960 1.26 +chr1 95756281 95929212 1.26 +chr1 115660513 115837500 1.26 +chr1 118390909 118391668 1.26 +chr1 206612616 206789340 1.26 +chr2 18490095 18666766 1.26 +chr2 106847377 107037217 1.26 +chr2 115568405 115730304 1.26 +chr2 130418007 130583404 1.26 +chr2 138069597 138255840 1.26 +chr2 150165935 150366790 1.26 +chr2 198816122 199013992 1.26 +chr3 30275518 30447552 1.26 +chr3 62016738 62182587 1.26 +chr3 107635524 107671459 1.26 +chr3 119676303 119822832 1.26 +chr3 147011398 147184818 1.26 +chr4 116321048 116472386 1.26 +chr4 138014454 138176309 1.26 +chr4 179303625 179403709 1.26 +chr5 38200579 38306406 1.26 +chr5 100711511 100906817 1.26 +chr7 34143563 34292666 1.26 +chr7 77710865 77899863 1.26 +chr7 110835679 110975058 1.26 +chr8 57918204 58101868 1.26 +chr8 140285430 140443678 1.26 +chr9 24877888 25069382 1.26 +chr9 29089125 29250390 1.26 +chr9 117088639 117257661 1.26 +chr10 37693680 37822602 1.26 +chr10 59005947 59241663 1.26 +chr10 128180787 128326333 1.26 +chr11 14830578 14997910 1.26 +chr11 27223566 27388824 1.26 +chr11 88739639 88927964 1.26 +chr12 28623872 28774457 1.26 +chr12 78967748 79140491 1.26 +chr12 82690981 82845109 1.26 +chr12 96065359 96245190 1.26 +chr13 68180810 68334608 1.26 +chr16 7933143 8120601 1.26 +chr16 77599573 77781796 1.26 +chr17 49753469 49922957 1.26 +chr17 67169020 67327134 1.26 +chr18 6752920 6933631 1.26 +chr18 22740630 22930994 1.26 +chr18 23968045 23968797 1.26 +chr18 64985278 65144962 1.26 +chr18 66753077 66929602 1.26 +chr21 26643246 26798261 1.26 +chrX 25247627 25248382 1.26 +chrX 94578207 94717584 1.26 +chrX 116471060 116589009 1.26 +chrY 6163949 6386452 1.26 +chr1 174297595 174457323 1.27 +chr1 236568362 236755350 1.27 +chr2 14744438 14891566 1.27 +chr2 89772211 89772694 1.27 +chr2 156410442 156580098 1.27 +chr4 108638904 108797073 1.27 +chr4 184232930 184404421 1.27 +chr5 49907782 50093906 1.27 +chr5 63325516 63499151 1.27 +chr5 98073104 98073634 1.27 +chr5 110059851 110181740 1.27 +chr5 125160995 125355179 1.27 +chr6 51267587 51472545 1.27 +chr6 80800588 80953422 1.27 +chr6 152324233 152411820 1.27 +chr7 80388111 80528144 1.27 +chr7 112279599 112435186 1.27 +chr7 112279641 112435146 1.27 +chr7 136606310 136780129 1.27 +chr8 20700584 20890068 1.27 +chr8 49937832 50097992 1.27 +chr8 77869835 78033428 1.27 +chr8 97773109 97944343 1.27 +chr8 108980084 109137354 1.27 +chr8 110436818 110612510 1.27 +chr8 111740331 111910034 1.27 +chr8 132474524 132638132 1.27 +chr9 11389658 11588107 1.27 +chr9 102748915 102917912 1.27 +chr11 38512156 38673888 1.27 +chr11 56238253 56397616 1.27 +chr11 78641711 78791818 1.27 +chr11 103803623 103964960 1.27 +chr12 74178839 74333757 1.27 +chr12 80827000 80973942 1.27 +chr13 25122500 25300811 1.27 +chr13 77954930 78112209 1.27 +chr17 22465099 22465647 1.27 +chr18 32884985 32885636 1.27 +chr18 62491384 62645196 1.27 +chr20 40680814 40792351 1.27 +chrX 67342298 67474154 1.27 +chrX 68455961 68651311 1.27 +chrX 87708952 87799321 1.27 +chrX 110343521 110498830 1.27 +chrX 120476674 120618187 1.27 +chrX 130385383 130512118 1.27 +chrX 142121643 142239873 1.27 +chrX 144849234 145030848 1.27 +chrX 150493632 150654782 1.27 +chr1 96868133 97055430 1.28 +chr1 97892548 98068736 1.28 +chr1 216532635 216718346 1.28 +chr1 238957519 239131143 1.28 +chr2 13501029 13665382 1.28 +chr2 94782159 94983765 1.28 +chr2 184329536 184509986 1.28 +chr2 209951701 210150189 1.28 +chr2 210167648 210331732 1.28 +chr3 3134172 3304965 1.28 +chr3 54612454 54778896 1.28 +chr4 62331997 62507133 1.28 +chr4 83182488 83348944 1.28 +chr4 110057781 110225193 1.28 +chr4 130001146 130195432 1.28 +chr4 134372150 134551414 1.28 +chr5 87524422 87715075 1.28 +chr5 152817276 152985131 1.28 +chr6 54720745 54910224 1.28 +chr6 91251335 91352253 1.28 +chr6 115381841 115506680 1.28 +chr6 128448619 128633208 1.28 +chr7 48821083 48911960 1.28 +chr7 67228859 67229276 1.28 +chr8 33027897 33028624 1.28 +chr8 73498872 73677480 1.28 +chr9 71447572 71636838 1.28 +chr9 79930264 79930787 1.28 +chr10 92100790 92284115 1.28 +chr11 81301979 81443218 1.28 +chr11 84737833 84908431 1.28 +chr12 37327276 37489074 1.28 +chr12 53730294 53914923 1.28 +chr12 69745683 69884764 1.28 +chr12 83409425 83566599 1.28 +chr13 80871787 81079018 1.28 +chr13 101318694 101507371 1.28 +chr14 46624364 46819048 1.28 +chr15 19791367 19970461 1.28 +chr15 23661899 23823333 1.28 +chr20 22476800 22597585 1.28 +chr20 53411918 53537915 1.28 +chr21 24478239 24645814 1.28 +chrX 72656424 72823653 1.28 +chrX 75925804 75926596 1.28 +chrX 123854281 123854885 1.28 +chr1 82426194 82603316 1.29 +chr1 217478892 217657878 1.29 +chr1 235542207 235725718 1.29 +chr2 17006354 17177187 1.29 +chr2 50950354 50950648 1.29 +chr2 139529650 139697746 1.29 +chr2 227039805 227205229 1.29 +chr3 28645912 28824317 1.29 +chr3 41172087 41354700 1.29 +chr4 5094075 5290909 1.29 +chr4 22163792 22309218 1.29 +chr4 68644475 68645263 1.29 +chr4 94662902 94838368 1.29 +chr4 94995921 94996745 1.29 +chr4 163467846 163611458 1.29 +chr5 7426758 7530112 1.29 +chr5 50061472 50168170 1.29 +chr5 117322499 117464139 1.29 +chr6 27422168 27563699 1.29 +chr6 103441042 103538793 1.29 +chr6 125389199 125541230 1.29 +chr6 141124549 141227062 1.29 +chr6 152125036 152333573 1.29 +chr7 54607124 54784717 1.29 +chr7 153721042 153901560 1.29 +chr8 73722232 73910787 1.29 +chr8 75573453 75722946 1.29 +chr8 122721235 122875986 1.29 +chr9 88883454 89058511 1.29 +chr10 20025423 20189475 1.29 +chr10 77812080 78010375 1.29 +chr10 117796864 117962944 1.29 +chr11 36670298 36852531 1.29 +chr12 85522910 85686077 1.29 +chr12 124939576 125073994 1.29 +chr13 53669880 53839191 1.29 +chr13 103214616 103420259 1.29 +chr14 38348741 38528993 1.29 +chr14 78800500 78957629 1.29 +chr15 44046046 44246628 1.29 +chr16 26107505 26268718 1.29 +chr16 32190271 32345186 1.29 +chr16 74775686 74929189 1.29 +chr17 11938351 12108967 1.29 +chr18 73908671 74017409 1.29 +chr20 58880457 58881237 1.29 +chrX 96349949 96509334 1.29 +chrY 16111429 16274219 1.29 +chr1 14511126 14690862 1.3 +chr1 55647414 55828567 1.3 +chr1 64206725 64303907 1.3 +chr1 161212962 161388982 1.3 +chr1 169324707 169484782 1.3 +chr1 230006922 230007507 1.3 +chr2 104867082 105049417 1.3 +chr2 132886154 133065879 1.3 +chr2 159108864 159300999 1.3 +chr2 179945358 180117962 1.3 +chr2 180322044 180501048 1.3 +chr2 188197236 188359167 1.3 +chr3 25551992 25552787 1.3 +chr3 28813297 28966841 1.3 +chr3 156642078 156810052 1.3 +chr3 191495503 191657007 1.3 +chr4 54822133 54996008 1.3 +chr4 125513760 125670332 1.3 +chr4 141603576 141774273 1.3 +chr4 177405907 177623780 1.3 +chr5 89564662 89728632 1.3 +chr5 154429161 154604443 1.3 +chr6 29321047 29495843 1.3 +chr6 49686804 49805309 1.3 +chr7 142857956 142929771 1.3 +chr8 2898662 3068843 1.3 +chr8 17782600 17783052 1.3 +chr9 2521111 2521805 1.3 +chr10 61902228 62069013 1.3 +chr10 110235647 110401173 1.3 +chr11 133984028 134156487 1.3 +chr12 66077114 66250763 1.3 +chr13 38450249 38608557 1.3 +chr13 46981034 46981776 1.3 +chr13 62944467 63099883 1.3 +chr13 67325775 67493542 1.3 +chr13 80711182 80918589 1.3 +chr14 86833840 86986583 1.3 +chr15 25800263 25980037 1.3 +chr15 35738732 35906861 1.3 +chr17 65314673 65421406 1.3 +chr18 5260008 5444854 1.3 +chr18 29304152 29472016 1.3 +chr20 16588770 16671907 1.3 +chrX 32142361 32277808 1.3 +chrX 86196923 86368724 1.3 +chrX 101574713 101575302 1.3 +chrX 138587336 138777284 1.3 +chr1 98681360 98839775 1.31 +chr2 15932536 16131069 1.31 +chr2 99955231 100127792 1.31 +chr2 206054394 206242414 1.31 +chr3 3156042 3326580 1.31 +chr3 110118090 110272403 1.31 +chr3 141263358 141264123 1.31 +chr3 149312730 149313543 1.31 +chr3 193387515 193550113 1.31 +chr4 80383883 80547099 1.31 +chr4 90923945 90924643 1.31 +chr4 103067218 103220633 1.31 +chr4 138445756 138622735 1.31 +chr4 166740143 166917702 1.31 +chr4 190520037 190668407 1.31 +chr5 3172060 3340431 1.31 +chr5 5500139 5686966 1.31 +chr5 40113144 40281744 1.31 +chr5 60735468 60897700 1.31 +chr6 96958620 97139258 1.31 +chr6 139955534 140131757 1.31 +chr7 46123519 46124125 1.31 +chr8 2303640 2411205 1.31 +chr8 2303675 2304133 1.31 +chr8 55504844 55667413 1.31 +chr9 16987010 17148443 1.31 +chr9 72488794 72637606 1.31 +chr9 78821714 78991548 1.31 +chr9 90604332 90759707 1.31 +chr10 33721598 33912463 1.31 +chr10 60813496 60976631 1.31 +chr11 50013045 50013630 1.31 +chr11 50013069 50182216 1.31 +chr12 22465615 22628113 1.31 +chr12 29550028 29733366 1.31 +chr12 29550028 29733343 1.31 +chr12 41495108 41495904 1.31 +chr12 99369399 99527535 1.31 +chr13 33995112 34159357 1.31 +chr13 111638158 111785816 1.31 +chr19 33293142 33399234 1.31 +chr20 7012254 7110690 1.31 +chr20 11812380 11979405 1.31 +chr20 15512841 15602650 1.31 +chr20 37522311 37666045 1.31 +chr20 50717376 50828721 1.31 +chrX 32347186 32479426 1.31 +chrX 111819716 111820515 1.31 +chrY 17166593 17345966 1.31 +chrY 20066103 20221627 1.31 +chr1 69875752 69876541 1.32 +chr1 75818828 76000027 1.32 +chr1 77147556 77248155 1.32 +chr1 197774987 197913697 1.32 +chr2 1612842 1783287 1.32 +chr2 1612887 1785220 1.32 +chr2 88939391 89105046 1.32 +chr2 131896303 132072566 1.32 +chr2 200012080 200209072 1.32 +chr3 104052002 104208021 1.32 +chr3 122378972 122541183 1.32 +chr3 138893130 139048105 1.32 +chr3 163463027 163648089 1.32 +chr4 36674866 36838139 1.32 +chr4 173240668 173410133 1.32 +chr4 181215989 181422694 1.32 +chr5 2570762 2724892 1.32 +chr5 39226026 39422569 1.32 +chr5 136365521 136366234 1.32 +chr6 101399800 101541206 1.32 +chr7 19402555 19587422 1.32 +chr7 26748780 26749335 1.32 +chr9 74808583 74977612 1.32 +chr9 89208558 89375340 1.32 +chr9 116920527 117122848 1.32 +chr9 116941815 117122920 1.32 +chr10 25969731 26132531 1.32 +chr10 87190070 87373214 1.32 +chr10 132635061 132635351 1.32 +chr12 73066799 73273410 1.32 +chr14 39369030 39537972 1.32 +chr17 50112691 50268063 1.32 +chr21 18224008 18383964 1.32 +chrX 5403668 5566016 1.32 +chrX 35271454 35446803 1.32 +chrX 43352640 43449221 1.32 +chrX 51220237 51394261 1.32 +chrX 56977002 57134049 1.32 +chrY 21140424 21306286 1.32 +chr1 91067460 91175162 1.33 +chr1 99670834 99835759 1.33 +chr1 170113616 170285659 1.33 +chr2 7661942 7829976 1.33 +chr2 22998605 23182948 1.33 +chr2 39502082 39691390 1.33 +chr2 48938583 49093125 1.33 +chr2 131163481 131364402 1.33 +chr2 137366984 137538813 1.33 +chr3 81809988 81810796 1.33 +chr3 96093190 96250469 1.33 +chr3 105815245 105971443 1.33 +chr3 198339454 198340252 1.33 +chr4 26746608 26747433 1.33 +chr4 75136723 75305982 1.33 +chr4 167845359 168019133 1.33 +chr5 66695993 66888015 1.33 +chr5 92986024 93173416 1.33 +chr5 107339380 107503471 1.33 +chr5 122463059 122626530 1.33 +chr5 146313344 146482537 1.33 +chr6 114513327 114699492 1.33 +chr6 122617580 122660549 1.33 +chr6 152094531 152188434 1.33 +chr6 154042939 154223658 1.33 +chr7 45903670 46077198 1.33 +chr7 89321351 89474116 1.33 +chr8 32530048 32715279 1.33 +chr8 104875758 105015139 1.33 +chr9 18579957 18743091 1.33 +chr9 112915772 113092453 1.33 +chr10 122467471 122467966 1.33 +chr11 112570898 112764095 1.33 +chr12 13505608 13642782 1.33 +chr12 114694711 114865871 1.33 +chr13 72606913 72777348 1.33 +chr13 74355387 74530080 1.33 +chr16 53823781 53975887 1.33 +chr18 6041920 6213636 1.33 +chr20 59426926 59538706 1.33 +chrX 74701380 74866100 1.33 +chrX 131183827 131343176 1.33 +chrX 136272151 136409212 1.33 +chrX 146504015 146682434 1.33 +chr1 50343064 50459755 1.34 +chr1 107386318 107478559 1.34 +chr1 119360541 119480610 1.34 +chr2 99112041 99262127 1.34 +chr2 148009604 148187956 1.34 +chr4 68459613 68647327 1.34 +chr4 187292822 187293270 1.34 +chr5 14496184 14655227 1.34 +chr6 44987777 45164728 1.34 +chr6 54092161 54254813 1.34 +chr7 8950609 9102382 1.34 +chr7 42647643 42810331 1.34 +chr7 70687274 70842877 1.34 +chr7 90960627 91079688 1.34 +chr7 115640211 115832430 1.34 +chr8 119798302 119961178 1.34 +chr9 16141186 16325481 1.34 +chr9 88787420 88940656 1.34 +chr9 116878175 117053374 1.34 +chr10 60095385 60270663 1.34 +chr10 81634066 81634692 1.34 +chr10 115833961 115994342 1.34 +chr10 116774295 116971212 1.34 +chr11 4745892 4948056 1.34 +chr11 27884866 28058344 1.34 +chr11 102869304 103066049 1.34 +chr11 133514734 133706144 1.34 +chr12 19280723 19281492 1.34 +chr12 78580511 78745903 1.34 +chr14 81191238 81376380 1.34 +chr15 71270959 71440808 1.34 +chr16 52428121 52626454 1.34 +chr16 52428138 52626428 1.34 +chr17 61192768 61362729 1.34 +chr19 6351112 6519252 1.34 +chr20 23543240 23638382 1.34 +chr20 40040735 40041178 1.34 +chrX 27025629 27210559 1.34 +chrY 18135971 18136728 1.34 +chr1 4476758 4608070 1.35 +chr1 83494525 83678458 1.35 +chr1 145434722 145622003 1.35 +chr1 157083023 157241521 1.35 +chr1 208589372 208727000 1.35 +chr1 230738768 230922038 1.35 +chr2 100117165 100292273 1.35 +chr2 147927761 148088945 1.35 +chr2 222493356 222674220 1.35 +chr3 158817691 158979197 1.35 +chr3 170366524 170543975 1.35 +chr3 187914183 187914670 1.35 +chr4 35206819 35400282 1.35 +chr4 52354879 52532792 1.35 +chr5 5879140 6020331 1.35 +chr5 36988555 37177051 1.35 +chr5 41662376 41829385 1.35 +chr5 93540085 93709247 1.35 +chr5 114441173 114603245 1.35 +chr5 178243610 178415257 1.35 +chr6 20630798 20759036 1.35 +chr6 46918436 47035989 1.35 +chr6 97824224 97971094 1.35 +chr6 164265539 164265787 1.35 +chr8 115760314 115915271 1.35 +chr9 7904520 8056890 1.35 +chr10 124083545 124268622 1.35 +chr11 55176672 55351006 1.35 +chr11 87911050 88091465 1.35 +chr11 102613735 102771566 1.35 +chr13 109210635 109357541 1.35 +chr14 31443192 31601773 1.35 +chr14 58323469 58470270 1.35 +chr14 60206359 60397740 1.35 +chr17 29512898 29667899 1.35 +chr18 6752916 6933612 1.35 +chr18 34297490 34473929 1.35 +chrX 65362816 65515922 1.35 +chr1 99670780 99835817 1.36 +chr1 142642787 142787183 1.36 +chr1 163383027 163567283 1.36 +chr1 218698700 218866017 1.36 +chr2 16384031 16384829 1.36 +chr2 21280646 21452847 1.36 +chr2 31383815 31583388 1.36 +chr2 204715244 204915179 1.36 +chr3 5501957 5696842 1.36 +chr3 6475280 6636703 1.36 +chr3 73463994 73612035 1.36 +chr3 100492740 100545594 1.36 +chr3 162044681 162207578 1.36 +chr3 174372436 174549467 1.36 +chr3 181855132 181855955 1.36 +chr4 4919317 5078151 1.36 +chr4 53060094 53207986 1.36 +chr4 82199002 82332981 1.36 +chr5 35388213 35513588 1.36 +chr6 20101886 20272058 1.36 +chr6 20101922 20272069 1.36 +chr6 79531824 79716854 1.36 +chr6 104973586 104973967 1.36 +chr7 41208853 41345448 1.36 +chr8 30945721 30946131 1.36 +chr8 68363210 68363812 1.36 +chr8 79633094 79802906 1.36 +chr8 133156783 133317994 1.36 +chr9 14419816 14586697 1.36 +chr10 27935308 28141733 1.36 +chr10 113860111 114056751 1.36 +chr11 35657471 35810959 1.36 +chr11 83021047 83205297 1.36 +chr11 101247718 101405229 1.36 +chr12 40716223 40716974 1.36 +chr13 20982693 21161824 1.36 +chr14 53557920 53721815 1.36 +chr14 80497781 80658071 1.36 +chr18 21431339 21583434 1.36 +chr18 21431373 21583433 1.36 +chr18 28257154 28415030 1.36 +chr18 50538203 50702047 1.36 +chr18 51445587 51648079 1.36 +chr18 74707626 74870951 1.36 +chr18 74707640 74870925 1.36 +chr19 22133724 22288871 1.36 +chr19 32891236 33099546 1.36 +chr19 36894580 36979850 1.36 +chr20 8156614 8328068 1.36 +chr20 19157422 19322167 1.36 +chrX 73387054 73556231 1.36 +chrX 103498150 103656697 1.36 +chrX 108414862 108517413 1.36 +chrX 146471609 146530342 1.36 +chrX 154334471 154482330 1.36 +chr1 18255166 18394021 1.37 +chr1 62131262 62274416 1.37 +chr1 94564158 94767971 1.37 +chr1 108330015 108512967 1.37 +chr1 225649911 225807728 1.37 +chr2 122827843 122978689 1.37 +chr3 34495510 34655794 1.37 +chr3 97022228 97176817 1.37 +chr3 108464746 108465227 1.37 +chr3 154955764 155113460 1.37 +chr3 177684955 177849617 1.37 +chr4 24079171 24250252 1.37 +chr4 89621629 89807633 1.37 +chr4 98622889 98832379 1.37 +chr4 98622912 98833448 1.37 +chr4 104276780 104433230 1.37 +chr4 175341179 175488565 1.37 +chr6 28599637 28767189 1.37 +chr6 65839991 65840775 1.37 +chr6 117915464 118016129 1.37 +chr8 49685131 49685914 1.37 +chr8 123525975 123697977 1.37 +chr9 106528618 106714989 1.37 +chr9 121736192 121918440 1.37 +chr10 17758041 17921747 1.37 +chr10 62853055 63041017 1.37 +chr10 82479195 82648679 1.37 +chr10 114874404 115058195 1.37 +chr11 96155681 96302306 1.37 +chr11 113740438 113878585 1.37 +chr12 129082307 129082773 1.37 +chr14 55649062 55827391 1.37 +chr14 97065728 97234127 1.37 +chr15 98650852 98814162 1.37 +chr16 33544449 33715230 1.37 +chr16 54425522 54559952 1.37 +chr18 11065439 11231074 1.37 +chr19 22570613 22729077 1.37 +chr20 6209827 6351381 1.37 +chr20 14556853 14704918 1.37 +chr20 57583626 57734574 1.37 +chr22 27103294 27317685 1.37 +chrX 42106039 42272650 1.37 +chrY 14466905 14467721 1.37 +chr1 59939363 60081655 1.38 +chr2 110173028 110348650 1.38 +chr2 119406291 119582870 1.38 +chr2 133513096 133688729 1.38 +chr2 169206769 169397061 1.38 +chr2 183145945 183311858 1.38 +chr2 228366700 228542036 1.38 +chr3 768376 958261 1.38 +chr3 66860333 67018889 1.38 +chr3 131457257 131616297 1.38 +chr3 137175189 137175931 1.38 +chr3 144846185 145031570 1.38 +chr3 193748164 193916444 1.38 +chr4 129322994 129526144 1.38 +chr4 143908621 143909369 1.38 +chr4 151531332 151532020 1.38 +chr5 173800904 173985872 1.38 +chr6 1241036 1397672 1.38 +chr6 75301349 75301837 1.38 +chr6 118676249 118787894 1.38 +chr6 131311502 131529403 1.38 +chr6 132102708 132208085 1.38 +chr6 132876727 132877160 1.38 +chr7 127107831 127266861 1.38 +chr8 58226066 58378511 1.38 +chr8 72600219 72767625 1.38 +chr8 90077495 90220312 1.38 +chr8 90077559 90220260 1.38 +chr8 119189820 119373386 1.38 +chr10 18842318 19038215 1.38 +chr12 37199957 37369589 1.38 +chr13 59662163 59888923 1.38 +chr13 107304438 107482925 1.38 +chr14 44775871 44957642 1.38 +chr15 93038000 93169754 1.38 +chr18 43069131 43274185 1.38 +chr20 21493032 21584641 1.38 +chr20 52655199 52655545 1.38 +chrX 7455566 7633822 1.38 +chrX 7633160 7633819 1.38 +chrX 44687783 44688510 1.38 +chrX 126290650 126444827 1.38 +chrY 12581761 12759800 1.38 +chr1 65371370 65460371 1.39 +chr1 143613581 143614341 1.39 +chr1 190142199 190319491 1.39 +chr1 207338657 207512534 1.39 +chr1 221352963 221534370 1.39 +chr1 241339986 241495285 1.39 +chr2 596981 744912 1.39 +chr2 41617751 41618169 1.39 +chr2 133791584 133959364 1.39 +chr2 135632580 135813995 1.39 +chr2 162682911 162846571 1.39 +chr2 240559993 240731366 1.39 +chr3 55606124 55797550 1.39 +chr3 68008697 68146379 1.39 +chr3 78395760 78564986 1.39 +chr3 99133076 99288212 1.39 +chr3 132613105 132784639 1.39 +chr3 170729328 170729895 1.39 +chr3 188693502 188694253 1.39 +chr4 31500911 31663999 1.39 +chr4 88909482 89097085 1.39 +chr4 101150836 101325274 1.39 +chr4 156373857 156552012 1.39 +chr5 27487185 27641956 1.39 +chr5 36275837 36380871 1.39 +chr5 102291745 102448736 1.39 +chr6 69938032 70037252 1.39 +chr6 146065393 146266178 1.39 +chr6 148426072 148609952 1.39 +chr7 16861638 17055792 1.39 +chr7 33718683 33897261 1.39 +chr7 46553602 46734204 1.39 +chr8 105218041 105402499 1.39 +chr9 78050479 78213759 1.39 +chr9 82003000 82148095 1.39 +chr10 47897443 48101918 1.39 +chr10 95774178 95924122 1.39 +chr10 123214100 123388480 1.39 +chr11 107487025 107646360 1.39 +chr11 107721404 107722054 1.39 +chr12 15977584 16171208 1.39 +chr12 20887518 20888220 1.39 +chr12 42413160 42579564 1.39 +chr12 64348953 64349751 1.39 +chr12 70624458 70765581 1.39 +chr12 88092437 88268617 1.39 +chr12 118059529 118238008 1.39 +chr13 52183624 52356899 1.39 +chr14 19570817 19730543 1.39 +chr14 24548302 24715386 1.39 +chr15 27464579 27613829 1.39 +chr20 28069297 28267569 1.39 +chrX 114736950 114737643 1.39 +chrY 6752267 6919657 1.39 +chr1 88594032 88769698 1.4 +chr1 176713719 176849348 1.4 +chr2 11122011 11294825 1.4 +chr2 190415420 190576094 1.4 +chr2 195231804 195390886 1.4 +chr2 205174891 205347264 1.4 +chr2 206824790 207023290 1.4 +chr2 217797626 218009975 1.4 +chr3 57370504 57371093 1.4 +chr3 137706751 137888006 1.4 +chr4 72239164 72405633 1.4 +chr4 81316240 81510343 1.4 +chr4 119888533 120067924 1.4 +chr4 164527243 164688085 1.4 +chr5 80045485 80216158 1.4 +chr5 82570102 82742889 1.4 +chr5 127603123 127759109 1.4 +chr6 71065292 71255915 1.4 +chr7 102949624 103109840 1.4 +chr7 107917096 108107231 1.4 +chr8 53938368 54105592 1.4 +chr8 120711400 120863829 1.4 +chr9 27417088 27590261 1.4 +chr9 102010405 102158130 1.4 +chr9 109317590 109318247 1.4 +chr10 21321389 21321781 1.4 +chr10 86226512 86390131 1.4 +chr10 119495422 119675862 1.4 +chr10 123333219 123506448 1.4 +chr10 129933095 130095088 1.4 +chr11 12368586 12369363 1.4 +chr11 114823737 114993848 1.4 +chr11 115161463 115257393 1.4 +chr12 25552960 25763321 1.4 +chr13 25898602 26065333 1.4 +chr13 26088616 26258539 1.4 +chr13 46548919 46699279 1.4 +chr13 50328528 50471832 1.4 +chr13 73124080 73302606 1.4 +chr15 85186245 85362775 1.4 +chr19 35629568 35766588 1.4 +chr19 42010298 42199663 1.4 +chr20 38174633 38311073 1.4 +chrX 2939526 3122894 1.4 +chrX 21448265 21601094 1.4 +chrX 56037806 56203576 1.4 +chrY 26625199 26798428 1.4 +chr1 34152124 34253960 1.41 +chr1 151011893 151012531 1.41 +chr1 191477741 191478202 1.41 +chr1 217542815 217743251 1.41 +chr1 246013826 246181016 1.41 +chr2 97472556 97702878 1.41 +chr3 24167259 24311442 1.41 +chr3 57297437 57406984 1.41 +chr3 127314118 127484692 1.41 +chr4 47264023 47432177 1.41 +chr4 54523234 54684548 1.41 +chr4 57434444 57434922 1.41 +chr4 74651758 74866776 1.41 +chr4 86385323 86558911 1.41 +chr4 111660518 111862003 1.41 +chr4 152766850 152907961 1.41 +chr6 7898162 8078796 1.41 +chr7 68251011 68418471 1.41 +chr7 119703966 119865943 1.41 +chr7 129658047 129658708 1.41 +chr8 9788890 9958420 1.41 +chr8 52776792 52920522 1.41 +chr8 81904402 82095683 1.41 +chr8 142612553 142790550 1.41 +chr10 4322070 4497152 1.41 +chr10 34543824 34750571 1.41 +chr11 11681667 11828279 1.41 +chr11 92720473 92905772 1.41 +chr12 76525022 76660149 1.41 +chr12 106096640 106265232 1.41 +chr15 74499839 74647142 1.41 +chr16 76434234 76606517 1.41 +chr18 168428 340445 1.41 +chr18 3863274 4046678 1.41 +chr18 66000802 66001509 1.41 +chrX 31427078 31590251 1.41 +chrX 132222083 132364969 1.41 +chr1 83633989 83797164 1.42 +chr1 175450745 175559022 1.42 +chr1 228957731 229156894 1.42 +chr2 6618691 6796407 1.42 +chr2 169796806 169937266 1.42 +chr2 169796837 169937261 1.42 +chr3 138231890 138416930 1.42 +chr3 190563330 190731931 1.42 +chr4 41340107 41527091 1.42 +chr4 88731532 88903393 1.42 +chr4 124508498 124677424 1.42 +chr6 28140043 28246315 1.42 +chr6 63764613 63922893 1.42 +chr6 152520069 152520603 1.42 +chr7 29450068 29624460 1.42 +chr7 69104887 69250742 1.42 +chr7 80906040 81083807 1.42 +chr7 102514284 102705988 1.42 +chr9 89829931 90026396 1.42 +chr10 42660894 42830194 1.42 +chr10 55596032 55596257 1.42 +chr11 5753721 5927294 1.42 +chr11 16626965 16801594 1.42 +chr11 93743479 93936234 1.42 +chr11 130972996 131149195 1.42 +chr12 64532787 64711495 1.42 +chr14 35613641 35778687 1.42 +chrY 15267652 15268179 1.42 +chr2 102713202 102893869 1.43 +chr2 210628113 210819973 1.43 +chr2 225556163 225739791 1.43 +chr2 235987180 236177484 1.43 +chr3 59701342 59875950 1.43 +chr3 81904890 82066186 1.43 +chr3 121651205 121651894 1.43 +chr3 125236945 125434205 1.43 +chr3 181453963 181647155 1.43 +chr4 77680481 77681311 1.43 +chr4 77782579 77943381 1.43 +chr4 139829369 139830163 1.43 +chr4 146558596 146724447 1.43 +chr4 166240548 166418926 1.43 +chr4 174136216 174137023 1.43 +chr5 69917500 70040712 1.43 +chr5 178243610 178455510 1.43 +chr6 24254617 24421543 1.43 +chr6 47870455 48002399 1.43 +chr6 116166727 116304219 1.43 +chr7 7176442 7292335 1.43 +chr7 56264968 56265717 1.43 +chr9 87247884 87248580 1.43 +chr9 87920679 88054773 1.43 +chr10 22399419 22575894 1.43 +chr11 123265105 123449312 1.43 +chr12 5163995 5164563 1.43 +chr12 97941659 98092858 1.43 +chr13 112025897 112208253 1.43 +chr14 50856299 51017418 1.43 +chr14 97835615 98011599 1.43 +chr15 79566763 79756514 1.43 +chr15 81302303 81488639 1.43 +chr15 99831524 100036168 1.43 +chr16 70783699 70973594 1.43 +chr17 28664934 28816067 1.43 +chr19 21345165 21504289 1.43 +chrX 6957788 7136832 1.43 +chr1 48002310 48202605 1.44 +chr1 182747653 182876763 1.44 +chr2 32483627 32649925 1.44 +chr2 173168767 173332941 1.44 +chr2 182686311 182851083 1.44 +chr2 198253860 198438524 1.44 +chr2 220237135 220431654 1.44 +chr4 55856185 55987940 1.44 +chr4 90029296 90200476 1.44 +chr4 131096084 131265054 1.44 +chr5 3426805 3559003 1.44 +chr5 79936331 80082792 1.44 +chr5 99850207 100001668 1.44 +chr6 40860381 41038343 1.44 +chr7 86933537 87062340 1.44 +chr7 131132439 131308575 1.44 +chr8 53078890 53252047 1.44 +chr8 82622728 82768419 1.44 +chr8 127563630 127764314 1.44 +chr9 97173317 97362893 1.44 +chr10 28985024 29145346 1.44 +chr10 33174826 33372767 1.44 +chr10 35838794 36004759 1.44 +chr10 108400155 108574810 1.44 +chr11 94849670 95025482 1.44 +chrX 102547184 102663453 1.44 +chr1 7059934 7186706 1.45 +chr1 181306539 181493746 1.45 +chr2 53784663 53945814 1.45 +chr2 119282256 119432338 1.45 +chr2 144588986 144779859 1.45 +chr3 115459232 115611158 1.45 +chr3 120110945 120111694 1.45 +chr3 140271723 140451428 1.45 +chr3 157113250 157286182 1.45 +chr4 176481822 176566236 1.45 +chr4 186650182 186822363 1.45 +chr5 60203778 60204237 1.45 +chr6 99896150 100068707 1.45 +chr6 119595033 119707570 1.45 +chr6 155961870 156122661 1.45 +chr7 64729936 64924073 1.45 +chr7 156044144 156202337 1.45 +chr8 65377850 65551264 1.45 +chr8 134871608 135047479 1.45 +chr9 85290912 85485840 1.45 +chr9 85859428 85948488 1.45 +chr9 89286203 89446273 1.45 +chr10 29841338 30014928 1.45 +chr10 92924800 93134214 1.45 +chr11 7401245 7402032 1.45 +chr11 111797499 111973902 1.45 +chr11 125962061 126160026 1.45 +chr12 75292746 75293210 1.45 +chr12 90372960 90550230 1.45 +chr14 56456968 56623373 1.45 +chr16 9634606 9827687 1.45 +chr17 61687769 61864883 1.45 +chr18 4926284 5083817 1.45 +chr18 10865355 11045527 1.45 +chr18 18072986 18222633 1.45 +chr19 20458451 20637715 1.45 +chr21 44672693 44846917 1.45 +chrX 47651895 47861422 1.45 +chrX 153861611 153862276 1.45 +chrY 2845480 3017508 1.45 +chr1 190332681 190333427 1.46 +chr2 134654316 134837314 1.46 +chr2 169789426 169953565 1.46 +chr3 31634839 31812066 1.46 +chr3 193863511 194023814 1.46 +chr4 99887462 100074907 1.46 +chr4 156624370 156778344 1.46 +chr5 73126031 73314954 1.46 +chr5 143298228 143480824 1.46 +chr5 145152023 145152440 1.46 +chr5 154417750 154418252 1.46 +chr6 73094225 73275800 1.46 +chr6 76566432 76685186 1.46 +chr6 97092117 97184863 1.46 +chr7 140977724 141092826 1.46 +chr8 5805153 6001336 1.46 +chr8 62411200 62572729 1.46 +chr9 73274833 73425564 1.46 +chr9 89446315 89614800 1.46 +chr10 38601694 38812575 1.46 +chr10 46999181 47159666 1.46 +chr10 131816154 131982049 1.46 +chr12 10845028 11018640 1.46 +chr12 62854244 62854997 1.46 +chr14 48879652 49050953 1.46 +chr14 100856556 101045895 1.46 +chr15 55906011 56096416 1.46 +chr16 8777563 8958332 1.46 +chr16 8957713 8958371 1.46 +chr18 17277421 17277775 1.46 +chr18 30699517 30859564 1.46 +chr20 10398757 10582167 1.46 +chrX 49444589 49644597 1.46 +chrY 7639125 7808753 1.46 +chr1 50594589 50802780 1.47 +chr1 63195526 63374556 1.47 +chr2 198360294 198526020 1.47 +chr3 32650193 32650799 1.47 +chr3 69383525 69541168 1.47 +chr3 170830820 170980842 1.47 +chr4 48524316 48683528 1.47 +chr6 126458972 126630444 1.47 +chr6 145034633 145199352 1.47 +chr6 168034078 168121846 1.47 +chr7 55009903 55010206 1.47 +chr7 103922036 104084401 1.47 +chr8 477653 658324 1.47 +chr9 86180561 86337441 1.47 +chr9 98020525 98190144 1.47 +chr10 8219578 8346526 1.47 +chr10 45156434 45157151 1.47 +chr10 100941261 101161313 1.47 +chr10 133629580 133778451 1.47 +chr11 14098699 14099413 1.47 +chr12 106899831 107062271 1.47 +chr12 117108431 117263859 1.47 +chr13 75289525 75439997 1.47 +chr13 94506903 94507480 1.47 +chr15 46837058 47011154 1.47 +chr15 53411855 53412638 1.47 +chr17 30523581 30685925 1.47 +chr18 9999822 10203371 1.47 +chr22 47740384 47838754 1.47 +chrX 64514128 64626801 1.47 +chr1 7641550 7843532 1.48 +chr1 46884600 47041895 1.48 +chr1 240002379 240203522 1.48 +chr2 70877489 70877986 1.48 +chr2 105850251 106014094 1.48 +chr2 151825968 151986544 1.48 +chr2 207687330 207861345 1.48 +chr3 8869825 9037105 1.48 +chr3 39667080 39836464 1.48 +chr3 115539059 115697740 1.48 +chr4 69694804 69912758 1.48 +chr4 78196278 78342827 1.48 +chr4 144367999 144533208 1.48 +chr4 144368003 144533204 1.48 +chr5 17110139 17284381 1.48 +chr5 70179512 70309338 1.48 +chr6 21512035 21698557 1.48 +chr6 56747408 56926686 1.48 +chr7 56639458 56803076 1.48 +chr7 97314295 97314786 1.48 +chr7 157549662 157756716 1.48 +chr9 93170720 93341568 1.48 +chr10 78781413 78996385 1.48 +chr10 85567652 85747565 1.48 +chr10 93738544 93911873 1.48 +chr12 45862358 46054280 1.48 +chr13 97268163 97467144 1.48 +chr15 46296351 46451066 1.48 +chr16 5440772 5644276 1.48 +chr18 46732299 46888453 1.48 +chr19 51083971 51239073 1.48 +chr21 32849426 33019513 1.48 +chrX 14761611 14921290 1.48 +chrX 37553903 37716432 1.48 +chrX 57816335 57960019 1.48 +chrX 100689549 100819590 1.48 +chrX 122218340 122344745 1.48 +chr1 87448603 87550379 1.49 +chr1 90457836 90458553 1.49 +chr1 173088739 173150600 1.49 +chr1 211480446 211481236 1.49 +chr1 219213887 219395917 1.49 +chr2 33026997 33209840 1.49 +chr2 78607472 78608236 1.49 +chr2 158164777 158343977 1.49 +chr2 234552685 234701715 1.49 +chr3 51558739 51559061 1.49 +chr3 112743047 112939558 1.49 +chr3 134620114 134771759 1.49 +chr4 106334083 106509296 1.49 +chr5 9436026 9606181 1.49 +chr5 10925654 10926306 1.49 +chr5 54562705 54764085 1.49 +chr6 170373241 170509779 1.49 +chr7 32275005 32459557 1.49 +chr7 107472400 107473057 1.49 +chr7 116327910 116481783 1.49 +chr9 4911574 5121406 1.49 +chr9 92185916 92362290 1.49 +chr9 107184004 107356920 1.49 +chr9 116607379 116766078 1.49 +chr11 22391987 22575753 1.49 +chr11 58901935 59116427 1.49 +chr11 93041021 93189057 1.49 +chr11 115593861 115757272 1.49 +chr12 9026914 9190867 1.49 +chr12 66740902 66931728 1.49 +chr13 31067305 31240545 1.49 +chr14 65027747 65175760 1.49 +chr14 65655482 65825533 1.49 +chr15 54901274 55123176 1.49 +chr16 26595131 26757807 1.49 +chr16 45067244 45231075 1.49 +chr17 5165932 5166614 1.49 +chr17 6065412 6217707 1.49 +chr17 20112029 20230301 1.49 +chr17 57872126 58038015 1.49 +chr18 8795684 8973445 1.49 +chr19 60914385 61088444 1.49 +chrX 18274436 18453317 1.49 +chr1 84728063 84892898 1.5 +chr1 163762157 163969574 1.5 +chr1 173739041 173886092 1.5 +chr2 44273710 44442507 1.5 +chr2 68007238 68163121 1.5 +chr2 195826912 195994198 1.5 +chr3 72455321 72633462 1.5 +chr3 106787070 106966498 1.5 +chr3 121148320 121326513 1.5 +chr5 62010355 62191843 1.5 +chr5 64142336 64300795 1.5 +chr5 73817311 73818090 1.5 +chr5 74847814 74970806 1.5 +chr5 79084142 79219538 1.5 +chr5 106747796 106920923 1.5 +chr5 108016191 108180023 1.5 +chr6 38990350 39129181 1.5 +chr6 86034671 86222195 1.5 +chr6 117034043 117169404 1.5 +chr7 157306039 157518326 1.5 +chr8 26206462 26412209 1.5 +chr9 3533199 3696631 1.5 +chr10 10830787 11020034 1.5 +chr10 97612977 97780354 1.5 +chr11 78034362 78206791 1.5 +chr12 49019643 49198741 1.5 +chr12 116586716 116748345 1.5 +chr12 129169567 129337248 1.5 +chr13 44121309 44301217 1.5 +chr15 32893109 33055196 1.5 +chr15 58181440 58324747 1.5 +chr17 69446939 69639765 1.5 +chr18 957589 958333 1.5 +chr18 41216566 41408713 1.5 +chr19 57561867 57749886 1.5 +chrX 106921602 107049634 1.5 +chrX 134430199 134632455 1.5 +chr1 146151766 146342686 1.51 +chr1 211426248 211601419 1.51 +chr2 39007277 39187672 1.51 +chr2 101278789 101279455 1.51 +chr2 127089207 127253875 1.51 +chr2 170471119 170690250 1.51 +chr2 200708454 200708868 1.51 +chr2 208715036 208876525 1.51 +chr2 222855303 223012639 1.51 +chr3 4492889 4493650 1.51 +chr3 154238282 154238738 1.51 +chr4 106900834 107064924 1.51 +chr4 184365850 184561840 1.51 +chr5 137785372 137785868 1.51 +chr5 159195225 159278790 1.51 +chr6 2316533 2413083 1.51 +chr6 15432597 15630682 1.51 +chr6 27152460 27241417 1.51 +chr6 45777376 45778135 1.51 +chr7 17756921 17931177 1.51 +chr7 55196125 55347935 1.51 +chr7 86906212 87054355 1.51 +chr7 133858782 134022392 1.51 +chr8 48686840 48844731 1.51 +chr8 59452746 59619686 1.51 +chr8 81493448 81648841 1.51 +chr8 104441333 104600167 1.51 +chr10 45425501 45633310 1.51 +chr11 2061361 2140325 1.51 +chr11 6747383 6920916 1.51 +chr11 95569628 95712661 1.51 +chr12 32552977 32729939 1.51 +chr13 23922512 24103118 1.51 +chr13 42155420 42301191 1.51 +chr15 31802458 31975417 1.51 +chr15 34700420 34860371 1.51 +chr16 22715683 22893642 1.51 +chr16 49618230 49794206 1.51 +chr18 45025819 45187979 1.51 +chr19 58164197 58336597 1.51 +chr20 47383084 47538301 1.51 +chrX 129672438 129773610 1.51 +chr1 77776731 77944517 1.52 +chr2 143499679 143667001 1.52 +chr2 143499703 143666971 1.52 +chr2 213160973 213347177 1.52 +chr2 216268805 216402525 1.52 +chr3 71622997 71816463 1.52 +chr3 121091572 121205097 1.52 +chr3 127928374 128108559 1.52 +chr3 135843691 136012579 1.52 +chr3 153771209 153915208 1.52 +chr5 34929543 35043272 1.52 +chr5 73872906 73873453 1.52 +chr5 86340653 86503917 1.52 +chr5 94527791 94663165 1.52 +chr5 124179322 124350943 1.52 +chr5 175139770 175266330 1.52 +chr6 110784447 110785221 1.52 +chr6 111684690 111819918 1.52 +chr6 137133725 137301607 1.52 +chr7 38185293 38185927 1.52 +chr7 45250029 45415202 1.52 +chr7 140911725 140979677 1.52 +chr8 124550520 124718428 1.52 +chr8 134125496 134255043 1.52 +chr9 5993718 6146499 1.52 +chr9 95248571 95419405 1.52 +chr10 14321741 14514319 1.52 +chr10 17006984 17194159 1.52 +chr10 25141002 25306820 1.52 +chr10 75725495 75896292 1.52 +chr11 531283 635941 1.52 +chr12 44739504 44903816 1.52 +chr13 27414257 27414749 1.52 +chr13 110508381 110681700 1.52 +chr14 30158279 30325251 1.52 +chr14 51766997 51964741 1.52 +chr14 66697392 66870415 1.52 +chr15 91550020 91714417 1.52 +chr16 17137869 17316758 1.52 +chr16 51466052 51670097 1.52 +chr16 81121446 81279554 1.52 +chr17 9466335 9611536 1.52 +chr17 15321328 15510169 1.52 +chr17 54149948 54280512 1.52 +chr19 42852412 43004473 1.52 +chrX 16028655 16242639 1.52 +chr1 34876929 35068982 1.53 +chr1 45348865 45523040 1.53 +chr1 92289815 92432460 1.53 +chr1 112311574 112481086 1.53 +chr1 117224641 117407078 1.53 +chr1 245125396 245280368 1.53 +chr2 157052263 157250015 1.53 +chr2 217224033 217224816 1.53 +chr3 133822599 133959231 1.53 +chr3 147710163 147710635 1.53 +chr4 76536770 76733068 1.53 +chr6 17735600 17880829 1.53 +chr6 112484860 112485640 1.53 +chr6 129664388 129861059 1.53 +chr6 150539782 150646794 1.53 +chr7 21588275 21784320 1.53 +chr7 23647193 23835807 1.53 +chr7 24489649 24659178 1.53 +chr8 39836613 39985589 1.53 +chr10 79687959 79844996 1.53 +chr10 120359649 120558865 1.53 +chr11 32892553 33037781 1.53 +chr12 101677282 101678070 1.53 +chr15 20363717 20542365 1.53 +chr18 70753437 70931323 1.53 +chr19 49045897 49206630 1.53 +chrX 45526480 45690299 1.53 +chr1 70717489 70814713 1.54 +chr1 87045853 87046483 1.54 +chr1 119501411 119708536 1.54 +chr1 165358469 165523632 1.54 +chr1 242491517 242647813 1.54 +chr2 36843825 37027743 1.54 +chr3 178148594 178149349 1.54 +chr4 4175379 4366937 1.54 +chr4 54045893 54204756 1.54 +chr5 95842138 96007440 1.54 +chr6 134784559 134979508 1.54 +chr6 143442949 143622404 1.54 +chr7 49785150 49947692 1.54 +chr8 95335790 95508534 1.54 +chr9 2136364 2296360 1.54 +chr10 31849338 32041147 1.54 +chr10 43864787 44024290 1.54 +chr10 63211719 63382032 1.54 +chr10 121490551 121670090 1.54 +chr11 18095206 18291271 1.54 +chr12 30740646 30939600 1.54 +chr14 93946930 93947549 1.54 +chr16 13833316 13983107 1.54 +chr17 2312031 2492162 1.54 +chr18 22104629 22273853 1.54 +chr18 32335236 32502050 1.54 +chr18 55962809 56151527 1.54 +chr20 61511174 61703926 1.54 +chr21 38078283 38078909 1.54 +chr22 26497322 26722737 1.54 +chr22 44414990 44543853 1.54 +chrX 71369082 71564002 1.54 +chrX 109597530 109598132 1.54 +chr1 66677291 66805722 1.55 +chr1 101075385 101207106 1.55 +chr1 116580533 116707521 1.55 +chr2 45253782 45449719 1.55 +chr2 68860263 69024185 1.55 +chr2 68860398 69024204 1.55 +chr2 152428325 152613810 1.55 +chr2 173995308 174174745 1.55 +chr2 190280619 190450126 1.55 +chr2 202418564 202419200 1.55 +chr3 99856407 100037503 1.55 +chr3 184141428 184305417 1.55 +chr5 170996378 171159953 1.55 +chr5 170996460 171159864 1.55 +chr8 48736329 48914366 1.55 +chr8 66009709 66192073 1.55 +chr9 74281674 74441354 1.55 +chr9 123651685 123824659 1.55 +chr10 48967298 49121505 1.55 +chr10 52151517 52328336 1.55 +chr10 103281369 103282021 1.55 +chr10 126605954 126760588 1.55 +chr11 2261639 2430734 1.55 +chr11 44879168 45046134 1.55 +chr11 68072319 68278585 1.55 +chr11 122124382 122299430 1.55 +chr12 14521944 14697706 1.55 +chr13 32790732 32947408 1.55 +chr15 70693934 70907539 1.55 +chr15 92324496 92491964 1.55 +chr17 43957711 44125703 1.55 +chr18 17227954 17228530 1.55 +chr18 52656262 52867700 1.55 +chr19 12006544 12159478 1.55 +chr20 36675397 36843266 1.55 +chrX 10267319 10437613 1.55 +chr1 42987601 43152105 1.56 +chr1 44062032 44062761 1.56 +chr1 196605835 196750596 1.56 +chr2 8025181 8218125 1.56 +chr2 24702643 24890535 1.56 +chr2 108344392 108345037 1.56 +chr3 42599838 42784967 1.56 +chr3 153502270 153697655 1.56 +chr3 153502331 153697658 1.56 +chr4 37192091 37353152 1.56 +chr4 146929902 147089763 1.56 +chr5 112488884 112695479 1.56 +chr5 157384965 157573351 1.56 +chr6 75947919 76073439 1.56 +chr7 21173940 21341345 1.56 +chr7 130078078 130270796 1.56 +chr8 8029024 8179923 1.56 +chr8 30816370 30990446 1.56 +chr8 90947388 91066955 1.56 +chr8 130562466 130702085 1.56 +chr9 15219945 15401987 1.56 +chr9 94281900 94449618 1.56 +chr9 124769226 124954691 1.56 +chr9 137112108 137333436 1.56 +chr10 13058510 13254677 1.56 +chr10 23123798 23319758 1.56 +chr10 126007431 126190348 1.56 +chr11 70596609 70789798 1.56 +chr12 3099629 3265088 1.56 +chr12 7918143 7919008 1.56 +chr12 7918177 7918979 1.56 +chr12 10122402 10291075 1.56 +chr13 22887527 23065746 1.56 +chr13 47897121 47897916 1.56 +chr14 53262095 53262810 1.56 +chr14 75653348 75833248 1.56 +chr14 75653368 75833225 1.56 +chr15 43413012 43613760 1.56 +chr15 48942176 49085328 1.56 +chr17 25071274 25261602 1.56 +chr18 9449205 9633157 1.56 +chr18 27366553 27569454 1.56 +chr22 42473589 42651043 1.56 +chrX 100064566 100164684 1.56 +chrY 13952205 14124161 1.56 +chr2 43518282 43698576 1.57 +chr2 71104398 71285149 1.57 +chr2 161596505 161786674 1.57 +chr2 239981214 240142437 1.57 +chr3 23574974 23575716 1.57 +chr3 151508467 151678827 1.57 +chr3 184802897 185002162 1.57 +chr6 52995970 53088442 1.57 +chr6 167790712 167979638 1.57 +chr7 137919306 138093868 1.57 +chr9 20172465 20351121 1.57 +chr10 94734535 94907303 1.57 +chr10 112963144 113136917 1.57 +chr11 94000385 94161720 1.57 +chr11 110352104 110510107 1.57 +chr11 129326505 129514541 1.57 +chr11 129513858 129514567 1.57 +chr12 38991144 39141476 1.57 +chr12 129648137 129839417 1.57 +chr13 31703968 31870598 1.57 +chr13 32842471 33025319 1.57 +chr13 44939386 45124346 1.57 +chr14 96318570 96319342 1.57 +chr14 100296015 100460986 1.57 +chr15 49615114 49801299 1.57 +chr16 46446494 46605357 1.57 +chr17 56657760 56847011 1.57 +chr18 11684280 11857259 1.57 +chr18 42350383 42350954 1.57 +chr18 46462730 46633371 1.57 +chr19 9041365 9197603 1.57 +chr20 1128792 1320966 1.57 +chr20 13564713 13774030 1.57 +chr1 37166110 37304834 1.58 +chr1 85106641 85107391 1.58 +chr1 113304939 113473675 1.58 +chr1 167621290 167722803 1.58 +chr2 27426604 27622433 1.58 +chr2 52630527 52788518 1.58 +chr2 109111217 109319695 1.58 +chr2 203317371 203318141 1.58 +chr3 9054010 9236123 1.58 +chr3 12590855 12735917 1.58 +chr3 179511139 179653380 1.58 +chr3 179511200 179653300 1.58 +chr4 152963195 153151424 1.58 +chr5 6515876 6708632 1.58 +chr5 58128952 58280997 1.58 +chr5 65059220 65239073 1.58 +chr5 179870172 180060082 1.58 +chr6 14536149 14673791 1.58 +chr6 52294378 52456670 1.58 +chr6 108026474 108214952 1.58 +chr8 10795780 10796421 1.58 +chr8 80726279 80893793 1.58 +chr9 96434170 96615981 1.58 +chr12 95366392 95367124 1.58 +chr13 28837843 29008462 1.58 +chr13 39497391 39498139 1.58 +chr14 20383517 20541450 1.58 +chr15 81186971 81351624 1.58 +chr17 36236729 36373682 1.58 +chr19 59041063 59235473 1.58 +chrX 12582350 12582807 1.58 +chrX 77111573 77195698 1.58 +chrX 108665629 108804379 1.58 +chrX 119444059 119591078 1.58 +chr1 29538364 29679865 1.59 +chr1 53077428 53193598 1.59 +chr1 181839038 182030843 1.59 +chr2 60501800 60696137 1.59 +chr2 86130524 86330997 1.59 +chr2 168899606 169089871 1.59 +chr2 172117570 172284258 1.59 +chr2 230386971 230387771 1.59 +chr3 192498455 192698313 1.59 +chr4 113577898 113578539 1.59 +chr5 72794439 72920149 1.59 +chr6 3292077 3471147 1.59 +chr6 10551747 10751130 1.59 +chr6 151442675 151561999 1.59 +chr9 97394031 97558497 1.59 +chr10 111497648 111663117 1.59 +chr11 34935748 34935995 1.59 +chr12 25269925 25447810 1.59 +chr12 94119164 94318610 1.59 +chr12 108178811 108360902 1.59 +chr13 48903929 49068849 1.59 +chr15 75185326 75348634 1.59 +chr17 3241556 3398985 1.59 +chr18 46310160 46473840 1.59 +chr18 58572418 58756466 1.59 +chr19 9688758 9689294 1.59 +chr20 45970037 46146093 1.59 +chr22 24560716 24734198 1.59 +chr1 9083937 9084661 1.6 +chr1 20428363 20560773 1.6 +chr1 58653712 58843461 1.6 +chr1 243222295 243431391 1.6 +chr2 223180236 223349542 1.6 +chr2 238922713 239091229 1.6 +chr3 14886339 14886921 1.6 +chr3 150340151 150506775 1.6 +chr4 71879366 72033748 1.6 +chr4 84626781 84801949 1.6 +chr5 118271452 118444329 1.6 +chr7 50442692 50617642 1.6 +chr7 139186224 139329952 1.6 +chr7 154264109 154414449 1.6 +chr8 57211535 57353615 1.6 +chr9 4819733 4991796 1.6 +chr9 101231164 101410218 1.6 +chr10 126760603 126954956 1.6 +chr11 8144704 8368299 1.6 +chr11 8634965 8795687 1.6 +chr11 124585420 124761526 1.6 +chr12 15496644 15658323 1.6 +chr13 19136604 19137297 1.6 +chr15 83685716 83883618 1.6 +chr16 919727 1097284 1.6 +chr17 10055064 10253595 1.6 +chr18 54867340 55049496 1.6 +chr18 54867345 55049546 1.6 +chr20 49489010 49588240 1.6 +chr21 29627040 29804741 1.6 +chrX 13814937 14015038 1.6 +chrX 19155814 19156191 1.6 +chrX 78252510 78419277 1.6 +chrX 153783896 153962641 1.6 +chr1 23210578 23311475 1.61 +chr1 67239552 67370740 1.61 +chr1 221890096 221891232 1.61 +chr2 178079912 178252288 1.61 +chr3 5027953 5180913 1.61 +chr3 15780432 15940568 1.61 +chr3 54007394 54188434 1.61 +chr3 124853438 125034577 1.61 +chr3 143078415 143237513 1.61 +chr3 144154881 144312253 1.61 +chr3 158164991 158323275 1.61 +chr4 47677803 47874056 1.61 +chr5 78198048 78198544 1.61 +chr5 79172773 79351477 1.61 +chr5 115966010 116138288 1.61 +chr6 5028697 5217232 1.61 +chr6 25482580 25677388 1.61 +chr6 26302873 26459898 1.61 +chr6 39977178 40078522 1.61 +chr6 144431624 144627356 1.61 +chr8 99215537 99366374 1.61 +chr9 99583829 99744600 1.61 +chr9 100604522 100754437 1.61 +chr9 114501030 114687327 1.61 +chr10 105830098 106023281 1.61 +chr12 5926585 5926958 1.61 +chr12 25394087 25557266 1.61 +chr12 94481227 94629374 1.61 +chr13 45700578 45869746 1.61 +chr14 69594037 69800258 1.61 +chr15 28255532 28543404 1.61 +chr15 67981331 68181027 1.61 +chr18 41984362 41984869 1.61 +chr19 56315908 56503568 1.61 +chr20 54786438 54948801 1.61 +chr21 25481427 25638118 1.61 +chr1 3214534 3355042 1.62 +chr2 30665578 30833166 1.62 +chr2 62595499 62742752 1.62 +chr2 171080964 171261141 1.62 +chr3 17809328 17960713 1.62 +chr3 142172454 142173036 1.62 +chr3 173181145 173347945 1.62 +chr3 173688016 173855765 1.62 +chr5 88149138 88309856 1.62 +chr5 169138790 169139510 1.62 +chr6 7427280 7575675 1.62 +chr8 126798041 126798844 1.62 +chr9 99080853 99246269 1.62 +chr9 100883044 101058940 1.62 +chr11 10289439 10454001 1.62 +chr11 18541571 18742470 1.62 +chr11 73706065 73867546 1.62 +chr11 101724679 101939102 1.62 +chr11 102010679 102011116 1.62 +chr11 124487806 124646248 1.62 +chr14 91389787 91561645 1.62 +chr17 28282582 28423748 1.62 +chr17 62106614 62246742 1.62 +chr18 53397307 53398022 1.62 +chrX 40755147 40944392 1.62 +chrX 70619662 70620234 1.62 +chrX 117482391 117604994 1.62 +chr1 110897345 111090422 1.63 +chr1 203522749 203523529 1.63 +chr2 47693987 47694786 1.63 +chr3 36866879 37036085 1.63 +chr4 7937465 8151463 1.63 +chr6 37889307 38009607 1.63 +chr6 130232415 130405154 1.63 +chr6 155289807 155290365 1.63 +chr7 19997089 20186026 1.63 +chr7 134684519 134842787 1.63 +chr9 19310518 19506748 1.63 +chr10 73834334 73834924 1.63 +chr10 74954877 74955199 1.63 +chr11 3573936 3757991 1.63 +chr11 9666315 9667022 1.63 +chr11 108226260 108393306 1.63 +chr12 51027859 51174187 1.63 +chr13 41262233 41262987 1.63 +chr14 34839207 35022652 1.63 +chr15 64985158 65179621 1.63 +chr16 4702956 4851459 1.63 +chr17 51993140 52172780 1.63 +chr19 5713441 5893445 1.63 +chr20 55914021 56121669 1.63 +chr22 30835441 30932660 1.63 +chrX 24176952 24341732 1.63 +chr1 93768326 93910259 1.64 +chr1 172077256 172222546 1.64 +chr1 180226264 180383894 1.64 +chr2 73670508 73888473 1.64 +chr2 98011957 98199214 1.64 +chr2 191161999 191354926 1.64 +chr3 10613350 10788431 1.64 +chr5 77156061 77320968 1.64 +chr6 32404070 32498505 1.64 +chr7 36168129 36337418 1.64 +chr7 76963834 77117182 1.64 +chr8 61676214 61851180 1.64 +chr8 141395906 141551776 1.64 +chr9 96375866 96551155 1.64 +chr9 111106642 111321237 1.64 +chr9 134611262 134820433 1.64 +chr10 101701720 101872241 1.64 +chr10 102895114 103074769 1.64 +chr11 69162501 69323924 1.64 +chr12 26494860 26676684 1.64 +chr12 56286203 56497599 1.64 +chr12 127784893 127938159 1.64 +chr13 76523008 76666874 1.64 +chr15 96860070 97061802 1.64 +chr15 97225257 97225937 1.64 +chr15 98388085 98574555 1.64 +chr16 12336699 12522785 1.64 +chr16 15983938 16191749 1.64 +chr16 48210521 48394084 1.64 +chr18 49786579 49971825 1.64 +chr20 2441363 2600552 1.64 +chr20 4636907 4766293 1.64 +chr20 20143748 20240111 1.64 +chrX 2789811 2790584 1.64 +chrX 46494737 46616168 1.64 +chr1 101220701 101413217 1.65 +chr1 157953509 158150295 1.65 +chr2 70374526 70514351 1.65 +chr4 38175739 38357481 1.65 +chr4 38356729 38357541 1.65 +chr4 148951122 148951822 1.65 +chr5 67498164 67677068 1.65 +chr5 112074297 112242499 1.65 +chr6 106038433 106153492 1.65 +chr7 101175494 101390682 1.65 +chr7 155193735 155348314 1.65 +chr8 21698622 21872574 1.65 +chr8 27277284 27462575 1.65 +chr8 102544799 102710032 1.65 +chr9 128256298 128420866 1.65 +chr10 70288208 70459535 1.65 +chr10 103204733 103388522 1.65 +chr12 91204898 91367961 1.65 +chr12 101481982 101482369 1.65 +chr13 19738155 19920905 1.65 +chr13 99118130 99118732 1.65 +chr15 41409192 41591912 1.65 +chr15 42593672 42756125 1.65 +chr15 57657503 57658171 1.65 +chr15 87338954 87527174 1.65 +chr16 10432694 10596527 1.65 +chr16 24637564 24856451 1.65 +chr17 33794727 33977073 1.65 +chr17 58063444 58237753 1.65 +chr17 68199396 68366917 1.65 +chr17 72107950 72271167 1.65 +chr20 56817450 56973377 1.65 +chr21 39567363 39702870 1.65 +chrX 2466238 2622009 1.65 +chrX 122792466 122990233 1.65 +chrY 2466238 2622009 1.65 +chr1 114409026 114555080 1.66 +chr2 42333837 42518734 1.66 +chr2 61421368 61422149 1.66 +chr2 95675764 95676554 1.66 +chr2 238251721 238252119 1.66 +chr3 45989826 46150224 1.66 +chr3 153173960 153322734 1.66 +chr5 42897645 42898147 1.66 +chr5 61659616 61832898 1.66 +chr5 67322915 67501695 1.66 +chr5 136892127 137064639 1.66 +chr5 157180105 157349297 1.66 +chr6 108524960 108647021 1.66 +chr7 27976454 28166805 1.66 +chr7 43425453 43594825 1.66 +chr8 67179019 67344055 1.66 +chr8 135577514 135739934 1.66 +chr10 26797315 26941968 1.66 +chr10 69441097 69601751 1.66 +chr10 90596444 90783689 1.66 +chr10 91315272 91463649 1.66 +chr12 127840870 127841450 1.66 +chr13 26843897 27009985 1.66 +chr13 51243111 51243933 1.66 +chr14 77328585 77513223 1.66 +chr17 42087499 42268621 1.66 +chr19 53269630 53431701 1.66 +chr20 3475997 3580840 1.66 +chr20 51561551 51690342 1.66 +chr21 38601021 38740380 1.66 +chr1 156056126 156056893 1.67 +chr1 222762761 222926484 1.67 +chr2 54748826 54918968 1.67 +chr2 175041508 175231393 1.67 +chr3 12420266 12590814 1.67 +chr3 113764691 113980300 1.67 +chr3 196585074 196768555 1.67 +chr6 6119858 6280511 1.67 +chr7 37407807 37597713 1.67 +chr7 76490702 76687246 1.67 +chr7 129214418 129214818 1.67 +chr8 29818492 30006914 1.67 +chr8 125633463 125840510 1.67 +chr9 131818912 132006030 1.67 +chr10 63780499 63955335 1.67 +chr10 99620303 99781492 1.67 +chr11 47289757 47477534 1.67 +chr11 85137550 85342098 1.67 +chr12 54839394 55007887 1.67 +chr13 47794572 47964053 1.67 +chr14 54376859 54570373 1.67 +chr15 40544105 40544306 1.67 +chr15 41680206 41851733 1.67 +chr16 70269439 70453945 1.67 +chr17 16526746 16718786 1.67 +chr17 19211694 19360796 1.67 +chr17 42875381 43043105 1.67 +chr17 57301021 57473353 1.67 +chr1 200974587 201176406 1.68 +chr1 233871784 234056082 1.68 +chr2 23670212 23881957 1.68 +chr2 23670217 23881959 1.68 +chr2 28412635 28572084 1.68 +chr2 178966383 179139203 1.68 +chr3 49222599 49393280 1.68 +chr3 123838303 124037428 1.68 +chr3 195537808 195699146 1.68 +chr5 171714383 171955264 1.68 +chr6 34918606 35051428 1.68 +chr6 88286044 88456239 1.68 +chr6 135475561 135663726 1.68 +chr6 149545624 149673890 1.68 +chr7 7947777 8126039 1.68 +chr7 105583898 105763396 1.68 +chr8 101279091 101431725 1.68 +chr11 33907530 34049097 1.68 +chr11 116204416 116368601 1.68 +chr12 102490337 102647688 1.68 +chr12 115487351 115487968 1.68 +chr12 131855464 132036472 1.68 +chr14 63342212 63500585 1.68 +chr14 99081660 99247736 1.68 +chr16 68164927 68165567 1.68 +chr16 70619026 70788093 1.68 +chr18 58408938 58578499 1.68 +chr20 34168760 34346661 1.68 +chr20 35834739 35931146 1.68 +chr20 39296098 39390691 1.68 +chr20 54455167 54566131 1.68 +chrX 12952343 13058902 1.68 +chrX 41108001 41241482 1.68 +chrX 54079685 54157798 1.68 +chr1 12351219 12462984 1.69 +chr1 16660236 16806953 1.69 +chr1 152333539 152334265 1.69 +chr1 177942471 178115476 1.69 +chr1 234121480 234285489 1.69 +chr2 37650188 37838486 1.69 +chr2 201556033 201729254 1.69 +chr2 213713161 213874167 1.69 +chr3 4490514 4652956 1.69 +chr3 47858146 48047371 1.69 +chr3 102991643 103144247 1.69 +chr3 109236225 109385450 1.69 +chr3 153144934 153314608 1.69 +chr3 170975630 171133894 1.69 +chr3 186612234 186768120 1.69 +chr5 142505106 142742827 1.69 +chr5 148642237 148846264 1.69 +chr6 107020676 107112028 1.69 +chr9 100705049 100705492 1.69 +chr9 135531598 135710682 1.69 +chr10 14952575 15140784 1.69 +chr10 89694476 89694968 1.69 +chr11 3854187 4005769 1.69 +chr12 100590406 100752075 1.69 +chr14 103217347 103382885 1.69 +chr15 39583755 39584254 1.69 +chr16 12001822 12176016 1.69 +chr16 29408714 29609851 1.69 +chr18 2479744 2665791 1.69 +chr20 24936530 25109994 1.69 +chrX 9560995 9684171 1.69 +chrX 135361340 135484299 1.69 +chr1 44428131 44608760 1.7 +chr1 51656028 51797618 1.7 +chr2 26821205 26990177 1.7 +chr2 45893818 46085030 1.7 +chr2 65586797 65745895 1.7 +chr2 197681931 197850356 1.7 +chr3 103202967 103378300 1.7 +chr3 114747418 114895013 1.7 +chr3 172293476 172503074 1.7 +chr4 10207878 10373914 1.7 +chr4 37870461 37871030 1.7 +chr5 141852268 141998278 1.7 +chr6 137962604 138142635 1.7 +chr7 30381583 30382165 1.7 +chr9 21157685 21158452 1.7 +chr9 130501524 130668171 1.7 +chr10 104652460 104813472 1.7 +chr11 1585438 1785273 1.7 +chr13 102258686 102424760 1.7 +chr14 68520131 68700019 1.7 +chr15 70142508 70313623 1.7 +chr15 80867047 81024206 1.7 +chr16 15411955 15573494 1.7 +chr16 67318530 67478802 1.7 +chr17 34978799 34979103 1.7 +chr19 51665825 51835148 1.7 +chr1 53509194 53688167 1.71 +chr1 53802593 53986331 1.71 +chr1 209767120 209981362 1.71 +chr2 96435158 96597950 1.71 +chr2 96435213 96597933 1.71 +chr2 237040133 237220803 1.71 +chr3 52658450 52743025 1.71 +chr5 131817042 131977001 1.71 +chr6 6217514 6382300 1.71 +chr6 160224184 160354066 1.71 +chr7 65934746 66092390 1.71 +chr7 128410279 128568571 1.71 +chr8 11802619 11803057 1.71 +chr8 23107287 23286627 1.71 +chr8 95978837 96177980 1.71 +chr9 91887431 91888260 1.71 +chr10 1206029 1370732 1.71 +chr10 30925759 30926104 1.71 +chr10 89597006 89798429 1.71 +chr11 119363929 119364258 1.71 +chr13 112739429 112901590 1.71 +chr15 50535080 50684151 1.71 +chr15 82907017 82907548 1.71 +chr17 32189512 32371445 1.71 +chr17 33321866 33322229 1.71 +chr20 31582329 31675113 1.71 +chrX 118438914 118566974 1.71 +chrX 129036286 129145716 1.71 +chr1 10075630 10169954 1.72 +chr1 112837726 112960659 1.72 +chr1 160195323 160195745 1.72 +chr1 166251815 166352458 1.72 +chr1 227623770 227744106 1.72 +chr2 196538333 196753496 1.72 +chr3 123486790 123658097 1.72 +chr5 10644878 10824706 1.72 +chr5 55546666 55547446 1.72 +chr5 138219217 138414953 1.72 +chr6 88811053 88923153 1.72 +chr9 115619446 115806020 1.72 +chr11 85720098 85720574 1.72 +chr12 121815155 121982522 1.72 +chr12 122658602 122876941 1.72 +chr17 23307568 23308267 1.72 +chr17 76183573 76354819 1.72 +chr20 17799943 17961951 1.72 +chr20 54471430 54472070 1.72 +chr1 15425500 15560568 1.73 +chr1 219787706 219961446 1.73 +chr2 9588319 9772865 1.73 +chr2 217062057 217226494 1.73 +chr3 128599495 128752041 1.73 +chr5 149219507 149349990 1.73 +chr6 12406236 12579760 1.73 +chr7 28459809 28603601 1.73 +chr9 70318675 70488655 1.73 +chr10 90435254 90615085 1.73 +chr11 47864340 48047114 1.73 +chr11 76074455 76232396 1.73 +chr12 52417126 52574777 1.73 +chr12 93046177 93206723 1.73 +chr12 112653797 112826527 1.73 +chr14 64199412 64200177 1.73 +chr17 41439719 41565972 1.73 +chr17 53510745 53677092 1.73 +chr18 72266630 72448118 1.73 +chr19 56161423 56314474 1.73 +chr19 63516095 63516696 1.73 +chrX 148437108 148548633 1.73 +chr1 19127387 19266745 1.74 +chr1 24351995 24532035 1.74 +chr1 24599009 24781759 1.74 +chr1 144149973 144150623 1.74 +chr1 202748911 202933682 1.74 +chr2 46105110 46297018 1.74 +chr3 51806508 51927169 1.74 +chr3 188090387 188091023 1.74 +chr4 8153524 8313477 1.74 +chr6 160501428 160502221 1.74 +chr7 98844510 99008804 1.74 +chr7 148462610 148626409 1.74 +chr8 28927769 29094558 1.74 +chr9 111914216 112092130 1.74 +chr9 132470436 132643733 1.74 +chr11 44421785 44422195 1.74 +chr11 82426839 82581132 1.74 +chr11 111095171 111276128 1.74 +chr12 11528806 11739539 1.74 +chr13 27824701 27970943 1.74 +chr14 74588657 74589380 1.74 +chr14 95199640 95365712 1.74 +chr16 14529557 14723144 1.74 +chr17 6893949 6894703 1.74 +chr17 6894000 7063252 1.74 +chr20 29778833 29779350 1.74 +chr20 31093124 31211874 1.74 +chr22 21400636 21552718 1.74 +chr1 918164 1145847 1.75 +chr1 21904425 22097047 1.75 +chr1 39307926 39450770 1.75 +chr1 198548931 198714422 1.75 +chr3 15628045 15764879 1.75 +chr3 187299633 187451862 1.75 +chr7 6234052 6234587 1.75 +chr8 8489834 8677590 1.75 +chr9 70739219 70739964 1.75 +chr9 125452145 125649123 1.75 +chr12 858648 1084584 1.75 +chr14 72043432 72226400 1.75 +chr15 60932064 61119675 1.75 +chr16 21936647 21937453 1.75 +chr17 8408778 8574309 1.75 +chr20 32566759 32687747 1.75 +chr21 36765659 36936719 1.75 +chr1 5866184 6042798 1.76 +chr1 17140091 17356497 1.76 +chr1 109536893 109537533 1.76 +chr1 151954735 152104802 1.76 +chr1 180895642 181046126 1.76 +chr2 101950474 102136502 1.76 +chr3 58149241 58327168 1.76 +chr4 185872338 186040225 1.76 +chr6 156885777 157054875 1.76 +chr6 159404565 159586466 1.76 +chr10 6277674 6476639 1.76 +chr10 72645818 72803738 1.76 +chr11 56867267 57054678 1.76 +chr11 72378220 72573798 1.76 +chr11 75107934 75273492 1.76 +chr12 1719982 1855189 1.76 +chr12 4174840 4341675 1.76 +chr12 31916950 32086540 1.76 +chr12 92994994 93158935 1.76 +chr15 65164181 65330882 1.76 +chr16 21304276 21473433 1.76 +chr16 65152761 65333416 1.76 +chr16 69164308 69365092 1.76 +chr17 74607599 74792971 1.76 +chr19 34956752 35094150 1.76 +chr20 29939223 30102016 1.76 +chr20 46989159 47129952 1.76 +chr20 48846717 48952485 1.76 +chr1 11483207 11618991 1.77 +chr1 29481197 29601498 1.77 +chr1 46259106 46477518 1.77 +chr1 148424844 148425333 1.77 +chr1 202668576 202832483 1.77 +chr10 98421164 98525673 1.77 +chr11 75965024 76130350 1.77 +chr12 56117999 56288135 1.77 +chr12 109799653 109949314 1.77 +chr15 65872293 66060824 1.77 +chr17 54693901 54872929 1.77 +chr22 36120293 36212700 1.77 +chr1 28834716 28835372 1.78 +chr2 25697038 25870821 1.78 +chr2 60927184 61134244 1.78 +chr3 9803852 9993473 1.78 +chr4 3367127 3557717 1.78 +chr6 135533146 135718367 1.78 +chr12 3999892 4189702 1.78 +chr17 25520145 25693265 1.78 +chr17 38137289 38327610 1.78 +chr17 41095612 41096064 1.78 +chr17 52616475 52760769 1.78 +chr17 63619442 63836050 1.78 +chr20 34886009 34886445 1.78 +chr20 42799956 42800600 1.78 +chr1 36604265 36743128 1.79 +chr2 29261117 29453743 1.79 +chr2 85513681 85689354 1.79 +chr5 131642147 131731304 1.79 +chr7 25786649 25787407 1.79 +chr7 65647563 65782464 1.79 +chr9 122851372 123029736 1.79 +chr10 76881864 76882526 1.79 +chr10 82019058 82019379 1.79 +chr11 1687193 1687713 1.79 +chr11 65876904 66035373 1.79 +chr11 117022863 117194572 1.79 +chr12 124056104 124234341 1.79 +chr17 38927004 39091531 1.79 +chr19 51988000 52169534 1.79 +chr1 200418761 200596067 1.8 +chr11 68509567 68701448 1.8 +chr14 22435305 22627125 1.8 +chr15 63070265 63255108 1.8 +chr17 7438479 7599079 1.8 +chr17 27416322 27596064 1.8 +chr17 39379048 39379525 1.8 +chr20 41668572 41668909 1.8 +chr20 59894661 60016627 1.8 +chr22 29897783 30032831 1.8 +chr1 159890346 159891087 1.81 +chr2 64478283 64684358 1.81 +chr5 149407439 149494253 1.81 +chr7 71970754 71971176 1.81 +chr10 125057762 125058197 1.81 +chr11 118726259 118898358 1.81 +chr13 98801820 98995341 1.81 +chr14 67822628 67995564 1.81 +chr15 88048919 88207829 1.81 +chr15 88048967 88207890 1.81 +chr16 30954589 31190651 1.81 +chr17 35682130 35682517 1.81 +chr17 37192625 37370848 1.81 +chr17 55350495 55498032 1.81 +chr17 59626448 59627242 1.81 +chr19 6523485 6701760 1.81 +chr20 275014 424783 1.81 +chr20 52147821 52257711 1.81 +chr22 38345454 38559089 1.81 +chr22 39986642 39987174 1.81 +chr1 89836304 90021445 1.82 +chr1 151954692 152104830 1.82 +chr1 179093657 179296252 1.82 +chr1 225105730 225267337 1.82 +chr2 8605936 8781369 1.82 +chr2 111491810 111673234 1.82 +chr3 10137523 10347462 1.82 +chr5 126045881 126232819 1.82 +chr7 5704715 5845078 1.82 +chr9 127232116 127408009 1.82 +chr11 63843056 64009168 1.82 +chr12 46666796 46667275 1.82 +chr12 119399640 119534882 1.82 +chr14 74297073 74478914 1.82 +chr16 67553299 67731519 1.82 +chr17 38221350 38425990 1.82 +chr20 44219240 44398282 1.82 +chr20 47068467 47243628 1.82 +chr1 159338985 159509999 1.83 +chr1 202583864 202785074 1.83 +chr3 120926480 121091560 1.83 +chr6 41997884 42172187 1.83 +chr8 61829586 62011027 1.83 +chr10 5677116 5844781 1.83 +chr12 111697547 111800218 1.83 +chr15 56535071 56693130 1.83 +chr16 49081891 49252657 1.83 +chr17 35065379 35227093 1.83 +chr17 69779007 69974731 1.83 +chr17 72330424 72523882 1.83 +chr18 19024453 19181574 1.83 +chrX 152359276 152538530 1.83 +chrX 152359331 152538528 1.83 +chr1 10986318 11098993 1.84 +chr2 43281993 43437205 1.84 +chr2 113063927 113232323 1.84 +chr3 13474877 13656208 1.84 +chr3 125910485 126076506 1.84 +chr3 129803278 129803833 1.84 +chr5 176797689 176798343 1.84 +chr6 11482211 11482901 1.84 +chr6 42887766 42955906 1.84 +chr7 44476301 44657269 1.84 +chr7 73732805 73733598 1.84 +chr7 74971406 75153739 1.84 +chr11 34986532 35140537 1.84 +chr11 127930601 128090813 1.84 +chr15 64818125 64818855 1.84 +chr16 83357071 83531668 1.84 +chr17 45862395 46041800 1.84 +chr1 32874030 33046828 1.85 +chr3 48357048 48512499 1.85 +chr6 6647138 6805876 1.85 +chr8 38289058 38488992 1.85 +chr11 120790891 120960982 1.85 +chr13 29914797 30048599 1.85 +chr14 49608076 49608742 1.85 +chr14 61029917 61208384 1.85 +chr17 4001094 4169898 1.85 +chr17 46171561 46380473 1.85 +chr20 43889716 43975657 1.85 +chr20 45646214 45732628 1.85 +chr22 33798821 33828353 1.85 +chr5 176391379 176542031 1.86 +chr6 36906248 37077769 1.86 +chr7 1611841 1807604 1.86 +chr9 133967495 133968038 1.86 +chr11 62270871 62271625 1.86 +chr12 120268349 120440851 1.86 +chr16 4441449 4598878 1.86 +chr16 66952765 67134295 1.86 +chr18 75538799 75701174 1.86 +chr20 32907707 33049748 1.86 +chr20 60992889 61121982 1.86 +chr21 43550744 43727862 1.86 +chr1 11765728 11932249 1.87 +chr1 31946121 32129493 1.87 +chr1 38263874 38379421 1.87 +chr9 132841764 133018395 1.87 +chr10 88054888 88222653 1.87 +chr11 60851894 61061395 1.87 +chr15 38766152 38944382 1.87 +chr15 71764680 71904803 1.87 +chr19 44539203 44539916 1.87 +chr20 41804350 41804986 1.87 +chr1 32275078 32440655 1.88 +chr1 205087972 205280504 1.88 +chr3 39135103 39311642 1.88 +chr6 158244856 158440772 1.88 +chr10 97210760 97374404 1.88 +chr11 66710945 66872614 1.88 +chr12 12759639 12961302 1.88 +chr12 103550316 103714143 1.88 +chr14 23511789 23749068 1.88 +chr14 102319935 102479461 1.88 +chr16 11795525 11796320 1.88 +chr19 52914892 53095367 1.88 +chr20 33649536 33824840 1.88 +chr20 48091894 48092464 1.88 +chr2 74335462 74515970 1.89 +chr5 71659788 71823446 1.89 +chr7 1080569 1081263 1.89 +chr7 99507554 99508034 1.89 +chr7 149650854 149838701 1.89 +chr9 126127930 126290748 1.89 +chr10 49230238 49437120 1.89 +chr12 51493439 51493869 1.89 +chr12 51493494 51676901 1.89 +chr12 130834289 130834783 1.89 +chr13 113770458 113954546 1.89 +chr15 76071170 76279435 1.89 +chr15 78200185 78353829 1.89 +chr15 89307983 89308593 1.89 +chr16 79834964 80033098 1.89 +chr17 7563870 7733801 1.89 +chr19 13100938 13216789 1.89 +chr19 16712801 16899931 1.89 +chr19 51023195 51151414 1.89 +chr22 22510132 22789754 1.89 +chr1 38051102 38051746 1.9 +chr1 45524166 45524938 1.9 +chr1 226955326 227082922 1.9 +chr2 219018312 219018782 1.9 +chr2 241620162 241796109 1.9 +chr6 33467588 33615999 1.9 +chr8 128735637 128884770 1.9 +chr10 89705200 89888619 1.9 +chr11 61394968 61395556 1.9 +chr15 64628666 64838826 1.9 +chr17 7251966 7438343 1.9 +chr19 4806367 5009897 1.9 +chr1 153110698 153292792 1.91 +chr9 6566990 6567805 1.91 +chr15 63873345 64073088 1.91 +chr17 801989 1008120 1.91 +chr1 153411319 153581460 1.92 +chr3 197391879 197392198 1.92 +chr8 56945439 57129714 1.92 +chr9 113691304 113877333 1.92 +chr10 135071590 135072261 1.92 +chr11 62229957 62403825 1.92 +chr11 64564018 64751914 1.92 +chr14 94901754 95067534 1.92 +chr22 35478918 35479413 1.92 +chr14 76653370 76654077 1.93 +chr14 92565537 92566181 1.93 +chr16 30539964 30747822 1.93 +chr17 23823737 23824010 1.93 +chr19 46479944 46480344 1.93 +chr19 53951044 54181980 1.93 +chr1 109816111 109932200 1.94 +chr7 72073705 72269192 1.94 +chr10 11805252 12011792 1.94 +chr14 68409982 68576217 1.94 +chr8 103802271 103959652 1.95 +chr9 129242753 129460262 1.95 +chr16 73977120 73977825 1.95 +chr17 1705906 1882486 1.95 +chr22 28379404 28551049 1.95 +chr22 39482928 39604225 1.95 +chr4 1595392 1784503 1.96 +chr17 35754749 35920940 1.96 +chr2 10189987 10393154 1.97 +chr12 119062359 119247085 1.97 +chr16 84379908 84540021 1.97 +chr17 33964845 34251503 1.97 +chr17 71063165 71247066 1.97 +chr19 10248852 10410584 1.97 +chr19 43766631 43893238 1.97 +chr1 40830627 40958746 1.98 +chr1 154215584 154399117 1.98 +chr10 71487253 71487667 1.98 +chr16 66199799 66393264 1.98 +chr17 17338154 17521462 1.98 +chr19 49912622 50122984 1.98 +chr21 42507308 42689289 1.98 +chr6 32031967 32200774 1.99 +chr12 51745498 51903199 1.99 +chr17 78374103 78374783 1.99 +chr2 96022598 96200121 2 +chr6 43965943 44149957 2 +chr7 2777815 2778183 2 +chr12 47562010 47745122 2 +chr15 38393810 38394432 2 +chr15 72737468 72953319 2 +chr22 39890713 40012950 2 +chr22 40022423 40142026 2 +chr12 55727385 55912942 2.01 +chr12 55727416 55727833 2.01 +chr16 55369598 55559295 2.01 +chr19 10699196 10699522 2.01 +chr2 231787962 231788629 2.02 +chr22 41529087 41529734 2.02 +chr1 27575366 27652624 2.03 +chr17 43957677 44125703 2.03 +chr18 54567090 54747580 2.03 +chr22 39701633 39877788 2.03 +chr1 204651282 204822719 2.04 +chr3 49843918 50006856 2.04 +chr19 13815786 13816238 2.04 +chr19 55588771 55753067 2.04 +chr1 9135174 9264249 2.05 +chr11 63543757 63721176 2.05 +chr19 2595842 2773344 2.05 +chr1 27728323 27853587 2.06 +chr16 80148499 80307562 2.06 +chr1 27231523 27352489 2.08 +chr7 99873610 100038722 2.08 +chr16 2067723 2267343 2.08 +chr20 30526698 30684121 2.08 +chr11 46216229 46383009 2.09 +chr16 431054 431761 2.1 +chr17 17940393 18114679 2.1 +chr19 12627174 12774845 2.1 +chr17 73636230 73813841 2.11 +chr11 71992715 72182751 2.13 +chr19 7474977 7636673 2.13 +chr19 54525652 54718270 2.15 +chr19 2928780 3109070 2.17 +chr1 25064047 25159219 2.18 +chr20 62108686 62284676 2.38 \ No newline at end of file diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 tool_conf.xml.sample --- a/tool_conf.xml.sample +++ b/tool_conf.xml.sample @@ -175,6 +175,7 @@ <section name="Regional Variation" id="regVar"><tool file="regVariation/windowSplitter.xml" /><tool file="regVariation/featureCounter.xml" /> + <tool file="regVariation/WeightedAverage.xml" /><tool file="regVariation/quality_filter.xml" /><tool file="regVariation/maf_cpg_filter.xml" /><tool file="regVariation/getIndels_2way.xml" /> @@ -194,8 +195,10 @@ </section><section name="Multiple regression" id="multReg"><tool file="regVariation/linear_regression.xml" /> + <tool file="regVariation/logistic_regression_vif.xml" /><tool file="regVariation/best_regression_subsets.xml" /><tool file="regVariation/rcve.xml" /> + <tool file="regVariation/partialR_square.xml" /></section><section name="Multivariate Analysis" id="multVar"><tool file="multivariate_stats/pca.xml" /> diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 tools/regVariation/WeightedAverage.py --- /dev/null +++ b/tools/regVariation/WeightedAverage.py @@ -0,0 +1,103 @@ +#!/usr/bin/env python +""" + +usage: %prog bed_file_1 bed_file_2 out_file + -1, --cols1=N,N,N,N: Columns for chr, start, end, strand in first file + -2, --cols2=N,N,N,N,N: Columns for chr, start, end, strand, name/value in second file +""" +from galaxy import eggs + +import collections +import sys, string +#import numpy +from galaxy import eggs +import pkg_resources +pkg_resources.require( "bx-python" ) +import sys, traceback, fileinput +from warnings import warn +from galaxy.tools.util.galaxyops import * +from bx.cookbook import doc_optparse + + +#export PYTHONPATH=~/galaxy/lib/ +#running command python WeightedAverage.py interval_interpolate.bed value_interpolate.bed interpolate_result.bed + +def stop_err(msg): + sys.stderr.write(msg) + sys.exit() + + +def FindRate(chromosome,start_stop,dictType): + OverlapList=[] + for tempO in dictType[chromosome]: + DatabaseInterval=[tempO[0],tempO[1]] + Overlap=GetOverlap(start_stop,DatabaseInterval) + if Overlap>0: + OverlapList.append([Overlap,tempO[2]]) + + if len(OverlapList)>0: + + SumRecomb=0 + SumOverlap=0 + for member in OverlapList: + SumRecomb+=member[0]*member[1] + SumOverlap+=member[0] + averageRate=SumRecomb/SumOverlap + + return averageRate + + else: + return 'NA' + + + +def GetOverlap(a,b): + return min(a[1],b[1])-max(a[0],b[0]) + +options, args = doc_optparse.parse( __doc__ ) + +try: + chr_col_1, start_col_1, end_col_1, strand_col1 = parse_cols_arg( options.cols1 ) + chr_col_2, start_col_2, end_col_2, strand_col2, name_col_2 = parse_cols_arg( options.cols2 ) + input1, input2, input3 = args +except Exception, eee: + print eee + stop_err( "Data issue: click the pencil icon in the history item to correct the metadata attributes." ) + + + +fd2=open(input2) +lines2=fd2.readlines() +RecombChrDict=collections.defaultdict(list) + +skipped=0 +for line in lines2: + temp=line.strip().split() + try: + assert float(temp[int(name_col_2)]) + except: + skipped+=1 + continue + tempIndex=[int(temp[int(start_col_2)]),int(temp[int(end_col_2)]),float(temp[int(name_col_2)])] + RecombChrDict[temp[int(chr_col_2)]].append(tempIndex) + +print "Skipped %d features with invalid values" %(skipped) + +fd1=open(input1) +lines=fd1.readlines() +finalProduct='' +for line in lines: + temp=line.strip().split('\t') + chromosome=temp[int(chr_col_1)] + start=int(temp[int(start_col_1)]) + stop=int(temp[int(end_col_1)]) + start_stop=[start,stop] + RecombRate=FindRate(chromosome,start_stop,RecombChrDict) + try: + RecombRate="%.4f" %(float(RecombRate)) + except: + RecombRate=RecombRate + finalProduct+=line.strip()+'\t'+str(RecombRate)+'\n' +fdd=open(input3,'w') +fdd.writelines(finalProduct) +fdd.close() diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 tools/regVariation/WeightedAverage.xml --- /dev/null +++ b/tools/regVariation/WeightedAverage.xml @@ -0,0 +1,71 @@ +<tool id="wtavg" name="Assign weighted-average" version="1.0.0"> + <description> of the values of features overlapping an interval </description> + <command interpreter="python">WeightedAverage.py $genomic_interval $genomic_feature $out_file1 -1 ${genomic_interval.metadata.chromCol},${genomic_interval.metadata.startCol},${genomic_interval.metadata.endCol},${genomic_interval.metadata.strandCol} -2 ${genomic_feature.metadata.chromCol},${genomic_feature.metadata.startCol},${genomic_feature.metadata.endCol},${genomic_feature.metadata.strandCol},${genomic_feature.metadata.nameCol}</command> + + <inputs> + <param format="interval" name="genomic_interval" type="data" label="Genomic intervals (first dataset)" help="Dataset missing? See Note below."/> + <param format="interval" name="genomic_feature" label="Genomic features (second dataset)" type="data" help="Make sure the value column is specified. See Note below." /> + </inputs> + <outputs> + <data format="input" name="out_file1" metadata_source="genomic_interval" /> + </outputs> + <tests> + <!-- Test data with valid values --> + <test> + <param name="genomic_interval" value="interval_interpolate.bed"/> + <param name="genomic_feature" value="value_interpolate.bed"/> + <output name="out_file1" file="interpolate_result.bed"/> + </test> + + </tests> + <help> + + +.. class:: infomark + +**What it does** + +For each interval in your first dataset, this tool calculates the weighted average value of the overlapping features in your second dataset. + +- When a genomic interval partially or totally overlaps a single genomic feature, the value of that genomic feature is assigned to the genomic interval. +- When a genomic interval partially or totally overlaps with more than one genomic features, the average of the values of the overlapping genomic features weighted by the corresponding number of overlapping bases is assigned to the genomic interval. +- When a genomic interval does not overlap with any genomic feature, 'NA' will be assigned as it's value. + +----- + +.. class:: warningmark + +**Note** + +The input datasets should be in **bed** or **interval** format. Please use "edit attributes"/pencil icon to specify the column containing the values for the features in the second dataset as **name/identifier** column. + +The output will contain all the columns in the first input plus a new column containing the assigned value for each interval. + +----- + +**Example** + +- Suppose our first dataset contains the following **genomic intervals**:: + + chr start stop + chr1 1000 2000 + chr1 3000 5000 + chr1 8000 9000 + +- and our second dataset contains the following **genomic features** each having an associated value (in fourth column) :: + + chr start stop name + chr1 900 1200 0.5 + chr1 2900 3100 0.2 + chr1 4800 5100 0.8 + +- For each **genomic interval** in our first dataset, this tool calculates the weighted average value of the overlapping **genomic features** in our second dataset :: + + chr1 1000 2000 0.5 + chr1 3000 5000 0.6 + chr1 8000 9000 NA + + + +</help> +</tool> \ No newline at end of file diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 tools/regVariation/logistic_regression_vif.py --- /dev/null +++ b/tools/regVariation/logistic_regression_vif.py @@ -0,0 +1,185 @@ +#!/usr/bin/env python + +from galaxy import eggs +import sys, string +from rpy import * +import numpy + +#export PYTHONPATH=~/galaxy/lib/ + +def stop_err(msg): + sys.stderr.write(msg) + sys.exit() +#infile = 'logreg_inp.tab' +#y_col=3 +#x_cols=[1,2,3] +#outfile='logreg_out.txt' +#python logistic_regression_vif.py logreg_inp.tab 4 1,2,3 logreg_out2.tabular # running test +infile = sys.argv[1] +y_col = int(sys.argv[2])-1 +x_cols = sys.argv[3].split(',') +outfile = sys.argv[4] + + +print "Predictor columns: %s; Response column: %d" %(x_cols,y_col+1) +fout = open(outfile,'w') +elems = [] +for i, line in enumerate( file ( infile )): + line = line.rstrip('\r\n') + if len( line )>0 and not line.startswith( '#' ): + elems = line.split( '\t' ) + break + if i == 30: + break # Hopefully we'll never get here... + +if len( elems )<1: + stop_err( "The data in your input dataset is either missing or not formatted properly." ) + +y_vals = [] +x_vals = [] + +for k,col in enumerate(x_cols): + x_cols[k] = int(col)-1 + x_vals.append([]) + +NA = 'NA' +for ind,line in enumerate( file( infile )): + if line and not line.startswith( '#' ): + try: + fields = line.split("\t") + try: + yval = float(fields[y_col]) + except: + yval = r('NA') + y_vals.append(yval) + for k,col in enumerate(x_cols): + try: + xval = float(fields[col]) + except: + xval = r('NA') + x_vals[k].append(xval) + except: + pass + +x_vals1 = numpy.asarray(x_vals).transpose() + +check1=0 +check0=0 +for i in y_vals: + if i == 1: + check1=1 + if i == 0: + check0=1 +if check1==0 or check0==0: + sys.exit("Warning: logistic regression must have at least two classes") + +for i in y_vals: + if i not in [1,0,r('NA')]: + print >>fout, str(i) + sys.exit("Warning: the current version of this tool can run only with two classes and need to be labeled as 0 and 1.") + + +dat= r.list(x=array(x_vals1), y=y_vals) +novif=0 +set_default_mode(NO_CONVERSION) +try: + linear_model = r.glm(r("y ~ x"), data = r.na_exclude(dat),family="binomial") + #r('library(car)') + #r.assign('dat',dat) + #r.assign('ncols',len(x_cols)) + #r.vif(r('glm(dat$y ~ ., data = na.exclude(data.frame(as.matrix(dat$x,ncol=ncols))->datx),family="binomial")')).as_py() + +except RException, rex: + stop_err("Error performing logistic regression on the input data.\nEither the response column or one of the predictor columns contain only non-numeric or invalid values.") +if len(x_cols)>1: + try: + + r('library(car)') + r.assign('dat',dat) + r.assign('ncols',len(x_cols)) + vif=r.vif(r('glm(dat$y ~ ., data = na.exclude(data.frame(as.matrix(dat$x,ncol=ncols))->datx),family="binomial")')) + except RException, rex: + print rex +else: + novif=1 + +set_default_mode(BASIC_CONVERSION) + +coeffs=linear_model.as_py()['coefficients'] +null_deviance=linear_model.as_py()['null.deviance'] +residual_deviance=linear_model.as_py()['deviance'] +yintercept= coeffs['(Intercept)'] +summary = r.summary(linear_model) +co = summary.get('coefficients', 'NA') +""" +if len(co) != len(x_vals)+1: + stop_err("Stopped performing logistic regression on the input data, since one of the predictor columns contains only non-numeric or invalid values.") +""" + +try: + yintercept = r.round(float(yintercept), digits=10) + pvaly = r.round(float(co[0][3]), digits=10) +except: + pass +print >>fout, "response column\tc%d" %(y_col+1) +tempP=[] +for i in x_cols: + tempP.append('c'+str(i+1)) +tempP=','.join(tempP) +print >>fout, "predictor column(s)\t%s" %(tempP) +print >>fout, "Y-intercept\t%s" %(yintercept) +print >>fout, "p-value (Y-intercept)\t%s" %(pvaly) + +if len(x_vals) == 1: #Simple linear regression case with 1 predictor variable + try: + slope = r.round(float(coeffs['x']), digits=10) + except: + slope = 'NA' + try: + pval = r.round(float(co[1][3]), digits=10) + except: + pval = 'NA' + print >>fout, "Slope (c%d)\t%s" %(x_cols[0]+1,slope) + print >>fout, "p-value (c%d)\t%s" %(x_cols[0]+1,pval) +else: #Multiple regression case with >1 predictors + ind=1 + while ind < len(coeffs.keys()): + try: + slope = r.round(float(coeffs['x'+str(ind)]), digits=10) + except: + slope = 'NA' + print >>fout, "Slope (c%d)\t%s" %(x_cols[ind-1]+1,slope) + try: + pval = r.round(float(co[ind][3]), digits=10) + except: + pval = 'NA' + print >>fout, "p-value (c%d)\t%s" %(x_cols[ind-1]+1,pval) + ind+=1 + +rsq = summary.get('r.squared','NA') + + +try: + rsq= r.round(float((null_deviance-residual_deviance)/null_deviance), digits=5) + null_deviance= r.round(float(null_deviance), digits=5) + residual_deviance= r.round(float(residual_deviance), digits=5) + + #rsq = r.round(float(rsq), digits=5) + +except: + pass + +print >>fout, "Null deviance\t%s" %(null_deviance) +print >>fout, "Residual deviance\t%s" %(residual_deviance) +print >>fout, "pseudo R-squared\t%s" %(rsq) +print >>fout, "\n" +print >>fout, 'vif' + +if novif==0: + py_vif=vif.as_py() + count=0 + for i in sorted(py_vif.keys()): + print >>fout,'c'+str(x_cols[count]+1) ,str(py_vif[i]) + count+=1 +elif novif==1: + print >>fout, "vif can calculate only when model have more than 1 predictor" diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 tools/regVariation/logistic_regression_vif.xml --- /dev/null +++ b/tools/regVariation/logistic_regression_vif.xml @@ -0,0 +1,74 @@ +<tool id="LogisticRegression" name="Perform Logistic Regression with vif" version="1.0.1"> + <description></description> + <command interpreter="python"> + logistic_regression_vif.py + $input1 + $response_col + $predictor_cols + $out_file1 + 1>/dev/null + </command> + <inputs> + <param format="tabular" name="input1" type="data" label="Select data" help="Dataset missing? See TIP below."/> + <param name="response_col" label="Response column (Y)" type="data_column" data_ref="input1" numerical="True"/> + <param name="predictor_cols" label="Predictor columns (X)" type="data_column" data_ref="input1" numerical="True" multiple="true" > + <validator type="no_options" message="Please select at least one column."/> + </param> + </inputs> + <outputs> + <data format="input" name="out_file1" metadata_source="input1" /> + + </outputs> + <requirements> + <requirement type="python-module">rpy</requirement> + </requirements> + <tests> + <test> + <param name="input1" value="logreg_inp.tabular"/> + <param name="response_col" value="4"/> + <param name="predictor_cols" value="1,2,3"/> + <output name="out_file1" file="logreg_out2.tabular"/> + + </test> + </tests> + <help> + + +.. class:: infomark + +**TIP:** If your data is not TAB delimited, use *Edit Datasets->Convert characters* + +----- + +.. class:: infomark + +**What it does** + +This tool uses the **'glm'** function from R statistical package to perform logistic regression on the input data. It outputs one file containing the summary statistics of the performed regression. Also, it calculates VIF(Variance Inflation Factor) with **'vif'** function from library (car) in R. + + +*R Development Core Team (2010). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org.* + +----- + +.. class:: warningmark + +**Note** + +- This tool currently treats all predictor variables as continuous numeric variables and response variable as categorical variable. Currently, the response variable can have only two classes, namely 0 and 1. The program will take 0 as base class. + +- Rows containing non-numeric (or missing) data in any of the chosen columns will be skipped from the analysis. + +- The summary statistics in the output are described below: + +- Pseudo R-squared: the proportion of model improvement from null model +- p-value: p-value for the z-test of the null hypothesis that the corresponding slope is equal to zero against the two-sided alternative. +- Coefficient indicates log ratio of (probability to be class 1 / probability to be class 0) + +- This tool also provides **Variance Inflation Factor or VIF** which quantifies the level of multicollinearity. The tool will automatic generate VIF if the model has more than one predictor. The higher the VIF, the higher is the multicollinearity. Multicollinearity will inflate standard error and reduce level of significance of the predictor. In the worst case, it can reverse direction of slope for highly correlated predictors if one of them is significant. A general thumb-rule is to use those predictors having VIF lower than 10 or 5. +- **vif** is calculated by + - First, regressing each predictor over all other predictors, and recording R-squared for each regression. + - Second, computing vif as 1/(1- R_squared) + + </help> +</tool> diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 tools/regVariation/partialR_square.py --- /dev/null +++ b/tools/regVariation/partialR_square.py @@ -0,0 +1,146 @@ +#!/usr/bin/env python + +from galaxy import eggs + +import sys, string +from rpy import * +import numpy + +#export PYTHONPATH=~/galaxy/lib/ +#running command python partialR_square.py reg_inp.tab 4 1,2,3 partialR_result.tabular + +def stop_err(msg): + sys.stderr.write(msg) + sys.exit() + +def sscombs(s): + if len(s) == 1: + return [s] + else: + ssc = sscombs(s[1:]) + return [s[0]] + [s[0]+comb for comb in ssc] + ssc + + +infile = sys.argv[1] +y_col = int(sys.argv[2])-1 +x_cols = sys.argv[3].split(',') +outfile = sys.argv[4] + +print "Predictor columns: %s; Response column: %d" %(x_cols,y_col+1) +fout = open(outfile,'w') + +for i, line in enumerate( file ( infile )): + line = line.rstrip('\r\n') + if len( line )>0 and not line.startswith( '#' ): + elems = line.split( '\t' ) + break + if i == 30: + break # Hopefully we'll never get here... + +if len( elems )<1: + stop_err( "The data in your input dataset is either missing or not formatted properly." ) + +y_vals = [] +x_vals = [] + +for k,col in enumerate(x_cols): + x_cols[k] = int(col)-1 + x_vals.append([]) + """ + try: + float( elems[x_cols[k]] ) + except: + try: + msg = "This operation cannot be performed on non-numeric column %d containing value '%s'." %( col, elems[x_cols[k]] ) + except: + msg = "This operation cannot be performed on non-numeric data." + stop_err( msg ) + """ +NA = 'NA' +for ind,line in enumerate( file( infile )): + if line and not line.startswith( '#' ): + try: + fields = line.split("\t") + try: + yval = float(fields[y_col]) + except Exception, ey: + yval = r('NA') + #print >>sys.stderr, "ey = %s" %ey + y_vals.append(yval) + for k,col in enumerate(x_cols): + try: + xval = float(fields[col]) + except Exception, ex: + xval = r('NA') + #print >>sys.stderr, "ex = %s" %ex + x_vals[k].append(xval) + except: + pass + +x_vals1 = numpy.asarray(x_vals).transpose() +dat= r.list(x=array(x_vals1), y=y_vals) + +set_default_mode(NO_CONVERSION) +try: + full = r.lm(r("y ~ x"), data= r.na_exclude(dat)) #full model includes all the predictor variables specified by the user +except RException, rex: + stop_err("Error performing linear regression on the input data.\nEither the response column or one of the predictor columns contain no numeric values.") +set_default_mode(BASIC_CONVERSION) + +summary = r.summary(full) +fullr2 = summary.get('r.squared','NA') + +if fullr2 == 'NA': + stop_error("Error in linear regression") + +if len(x_vals) < 10: + s = "" + for ch in range(len(x_vals)): + s += str(ch) +else: + stop_err("This tool only works with less than 10 predictors.") + +print >>fout, "#Model\tR-sq\tpartial_R_Terms\tpartial_R_Value" +all_combos = sorted(sscombs(s), key=len) +all_combos.reverse() +for j,cols in enumerate(all_combos): + #if len(cols) == len(s): #Same as the full model above + # continue + if len(cols) == 1: + x_vals1 = x_vals[int(cols)] + else: + x_v = [] + for col in cols: + x_v.append(x_vals[int(col)]) + x_vals1 = numpy.asarray(x_v).transpose() + dat= r.list(x=array(x_vals1), y=y_vals) + set_default_mode(NO_CONVERSION) + red = r.lm(r("y ~ x"), data= dat) #Reduced model + set_default_mode(BASIC_CONVERSION) + summary = r.summary(red) + redr2 = summary.get('r.squared','NA') + try: + partial_R = (float(fullr2)-float(redr2))/(1-float(redr2)) + except: + partial_R = 'NA' + col_str = "" + for col in cols: + col_str = col_str + str(int(x_cols[int(col)]) + 1) + " " + col_str.strip() + partial_R_col_str = "" + for col in s: + if col not in cols: + partial_R_col_str = partial_R_col_str + str(int(x_cols[int(col)]) + 1) + " " + partial_R_col_str.strip() + if len(cols) == len(s): #full model + partial_R_col_str = "-" + partial_R = "-" + try: + redr2 = "%.4f" %(float(redr2)) + except: + pass + try: + partial_R = "%.4f" %(float(partial_R)) + except: + pass + print >>fout, "%s\t%s\t%s\t%s" %(col_str,redr2,partial_R_col_str,partial_R) diff -r e13fb2823aa4cbd835de79c89f53feca37ab7969 -r 0622b489a0f17f4b328473a2930f335d8d879fa1 tools/regVariation/partialR_square.xml --- /dev/null +++ b/tools/regVariation/partialR_square.xml @@ -0,0 +1,68 @@ +<tool id="partialRsq" name="Compute partial R square" version="1.0.0"> + <description></description> + <command interpreter="python"> + partialR_square.py + $input1 + $response_col + $predictor_cols + $out_file1 + 1>/dev/null + </command> + <inputs> + <param format="tabular" name="input1" type="data" label="Select data" help="Dataset missing? See TIP below."/> + <param name="response_col" label="Response column (Y)" type="data_column" data_ref="input1" /> + <param name="predictor_cols" label="Predictor columns (X)" type="data_column" data_ref="input1" multiple="true"> + <validator type="no_options" message="Please select at least one column."/> + </param> + </inputs> + <outputs> + <data format="input" name="out_file1" metadata_source="input1" /> + </outputs> + <requirements> + <requirement type="python-module">rpy</requirement> + </requirements> + <tests> + <!-- Test data with vlid values --> + <test> + <param name="input1" value="regr_inp.tabular"/> + <param name="response_col" value="3"/> + <param name="predictor_cols" value="1,2"/> + <output name="out_file1" file="partialR_result.tabular"/> + </test> + + </tests> + <help> + +.. class:: infomark + +**TIP:** If your data is not TAB delimited, use *Edit Datasets->Convert characters* + +----- + +.. class:: infomark + +**What it does** + +This tool computes the Partial R squared for all possible variable subsets using the following formula: + +**Partial R squared = [SSE(without i: 1,2,...,p-1) - SSE (full: 1,2,..,i..,p-1) / SSE(without i: 1,2,...,p-1)]**, which denotes the case where the 'i'th predictor is dropped. + + + +In general, **Partial R squared = [SSE(without i: 1,2,...,p-1) - SSE (full: 1,2,..,i..,p-1) / SSE(without i: 1,2,...,p-1)]**, where, + +- SSE (full: 1,2,..,i..,p-1) = Sum of Squares left out by the full set of predictors SSE(X1, X2 … Xp) +- SSE (full: 1,2,..,i..,p-1) = Sum of Squares left out by the set of predictors excluding; for example, if we omit the first predictor, it will be SSE(X2 … Xp). + + +The 4 columns in the output are described below: + +- Column 1 (Model): denotes the variables present in the model +- Column 2 (R-sq): denotes the R-squared value corresponding to the model in Column 1 +- Column 3 (Partial R squared_Terms): denotes the variable/s for which Partial R squared is computed. These are the variables that are absent in the reduced model in Column 1. A '-' in this column indicates that the model in Column 1 is the Full model. +- Column 4 (Partial R squared): denotes the Partial R squared value corresponding to the variable/s in Column 3. A '-' in this column indicates that the model in Column 1 is the Full model. + +*R Development Core Team (2010). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org.* + + </help> +</tool> Repository URL: https://bitbucket.org/galaxy/galaxy-central/ -- This is a commit notification from bitbucket.org. You are receiving this because you have the service enabled, addressing the recipient of this email.