3月3日作業履歴
100Kでアジレントのプローブを選択予定、プローブがそれぞれ100K以上離れている(多分NNNNNの部分)
他も作ってみる。
[minako1@minako agilent]# grep "[[:space:]]chr20[[:space:]]" agilentCgh244a.txt | perl 20100301.pl - 100000 > chr20_100k_dist.v2.bed
[minako1@minako agilent]# grep "[[:space:]]chrX[[:space:]]" agilentCgh244a.txt | perl 20100301.pl - 100000 > chrX_100k_dist.v2.bed
catで連結
[minako1@minako chr_100k_dist.v2]# cat chr1_100k_dist.v2.bed chr2_100k_dist.v2.bed chr3_100k_dist.v2.bed chr4_100k_dist.v2.bed chr5_100k_dist.v2.bed chr6_100k_dist.v2.bed chr7_100k_dist.v2.bed chr8_100k_dist.v2.bed chr9_100k_dist.v2.bed chr10_100k_dist.v2.bed chr11_100k_dist.v2.bed chr12_100k_dist.v2.bed chr13_100k_dist.v2.bed chr14_100k_dist.v2.bed chr15_100k_dist.v2.bed chr16_100k_dist.v2.bed chr17_100k_dist.v2.bed chr18_100k_dist.v2.bed chr19_100k_dist.v2.bed chr20_100k_dist.v2.bed chrX_100k_dist.v2.bed > allchr_100k_dist.v2.bed
.bedのままでratのgenome graphにアップロード。
