rat chr20のn...の数いくつ？
[Minako@oyoyo ~/test]# perl probe.distance.v2.a.pl chr20_nnncut.gff
55263898 55263898 1
55263901 55263903 3
55263905 55263906 2
55263908 55263908 1
55263911 55263913 3
55263916 55263919 4
55263922 55263922 1
55263925 55263925 1
55263927 55263938 12
55263940 55263941 2
55263943 55263957 15
55263959 55263976 18
55263980 55263980 1
55263982 55263989 8
start end n個数
これをucscに前後1つづつ伸ばして入れてみると、たしかに。
[Minako@oyoyo ~/test]# perl probe.distance.v2.a.pl chr20_nnncut.gff | perl /home/Minako/test/agilent/bar2.pl - 10 | head
0 11534 #1-9個のnの個数
10 1013
20 302
30 152
40 105
50 1424
60 49
70 57
80 38
90 53
[Minako@oyoyo ~/test]# perl probe.distance.v2.a.pl chr20_nnncut.gff | perl /home/Minako/test/agilent/bar2.pl - 10 | tail
95380 1
100000 2
104320 1
135420 1
146000 1
162000 1
186000 1
188000 1
212000 1
386000 1 #結構長いnnnn...
[Minako@oyoyo ~/test]# perl probe.distance.v2.a.pl chr20_nnncut.gff | perl /home/Minako/test/agilent/bar2.pl - 1 | head -20
1 5697 #n1つが最も多い
2 1805
3 1291
4 740
5 623
6 497
7 384
8 269
9 228
10 192
11 134
12 120
13 111
14 98
15 91
16 69
17 69
18 76
19 53
20 36
[Minako@oyoyo ~/test]# perl probe.distance.v2.a.pl chr20_nnncut.gff | perl /home/Minako/test/agilent/bar2.pl - 10000
0 16137
10000 57
20000 14
30000 7
40000 2
50000 12
60000 3
70000 1
80000 1
90000 1
100000 3
130000 1
140000 1
160000 1
180000 2
210000 1
380000 1
[Minako@oyoyo ~/test]# perl probe.distance.v2.a.pl chr20_nnncut.gff | perl /home/Minako/test/agilent/bar2.pl - 10000 > chr20_nnn_bin10k.txt
なんか当たり前のような、、こんなのつくってどうするんだ？
