在展示染色體信息時,如果想要重點展示其中某一段區(qū)域的信息,可以借助zooms
來實現(xiàn)。zooms
起到一個伸縮的功能,將原本的區(qū)域放大或者縮小。
配置文件寫法如下:
zooms
由多個zoom
構(gòu)成,每個zoom
必須具備以下幾個參數(shù)
chr
start
end
-
scale
其中chr
, start
, end
這三個參數(shù)指定需要縮放的染色體區(qū)域,scale
指定需要縮放的倍數(shù)。
看一個具體的例子:
在上面的示意圖中,對染色體多處區(qū)域進行了縮放。其中對1號染色體上的部分區(qū)域進行了放大,對2號染色體上的部分區(qū)域進行了縮小。
1號染色體100-120Mb區(qū)域放大2倍,接下來的區(qū)域依次類推,分別放大了3倍, 4倍,5倍;對應(yīng)配置文件的寫法為
2號染色體100-120Mb 區(qū)域縮小為原來的0.5倍,接下來的區(qū)域依次類推,分別縮小為原來0.25, 0.1,0.25, 0.5;對應(yīng)配置文件的寫法為
在上述的例子中,需要縮放的區(qū)域之間沒有重疊,當(dāng)這些區(qū)域之間以后重疊時,重疊區(qū)的縮放比例以最高的為準(zhǔn)。比如下面的縮放比例
實際處理時會變成如下的縮放比例
對于沒有爭議的區(qū)域,就按照指定的比例縮放,對于指定了多個縮放比例的overlap
區(qū),以最大的數(shù)值作為最終的縮放比例。
在zoom
相關(guān)參數(shù)中,還有兩個很特別的參數(shù)
smooth_distance
-
smooth_steps
寫法如下:
這兩個參數(shù)的作用是控制鄰近區(qū)域的縮放情況。不指定這兩個參數(shù)時,只有目標(biāo)區(qū)域會進行縮放;指定這兩個參數(shù)后,會對上下游區(qū)域進行縮放。
smooth_distance
控制需要縮放的上下游區(qū)域的距離,2r
代表是目標(biāo)區(qū)域的兩倍,由于目標(biāo)區(qū)域為120u到125u
, 長度為5M,需要縮放的上下游區(qū)域就是2 X 5 = 10M,這些區(qū)域也需要方法scale
倍數(shù); smooth_steps
控制上下游區(qū)域縮放的步數(shù),上下游區(qū)域的縮放比例是一個線性增長的關(guān)系,在下面的示意圖中,可以看到,1號染色體目的區(qū)域上游的10M區(qū)域110-120
, 也放大為原來的10倍,其縮放不是整個區(qū)域直接放大10倍,而是劃分成了10個子的區(qū)域,每個區(qū)域縮放比例逐漸遞增,有1.8遞增到9.2,總和加起來是原先的10倍。
對于zooms
而言,其參數(shù)一共只有6個,具體的含義上文也給很給出了詳細的解釋。但是對于染色體的縮放,除了zooms
之外,還有其他的參數(shù)也可以實現(xiàn)
1. 對整個染色體進行縮放
通過chromosomes_scale
參數(shù)來實現(xiàn),寫法如下
chromosomes_scale = hs1:0.2;hs2:0.2;hs3:0.2;hs8:5;hs9:5;hs10:5
hs1:0.2
表示將1號染色體縮小為原來的0.2倍,對于每條染色體的縮放比例,采用分號分隔。
2. 對染色體部分區(qū)域進行縮放
寫法如下
chromosomes = hs1[a]:0-20;hs2[b]:0-20;hs1[c]:20-40;hs2[d]:20-40;hs1[e]:40-60;hs2[f]:40-60
chromosomes_scale = a:0.5;b:0.5;e:5;f:5
以hs1
為例,在chromosomes
中,指定了a
, c
, e
3個區(qū)間,在chromosomes_scale
中,對這3個區(qū)間進行了縮放。