区别intrachr和outrachr
有的时候我们想把染色体内部的连线与不同染色体间的连线区别开来,在这个例子中我们将使用rule来帮助实现这一目的. 先看看效果图:
CIRCOS.CONF 完整
···
<links>
z = 0
radius = 0.90r
crest = 1
color = vvlgrey
bezier_radius = 0r
bezier_radius_purity = 0.50
thickness = 2
<link>
file = data/5/segdup.txt
<rules>
<rule>
condition = var(intrachr) #是否在染色体内部
condition = abs(var(pos1) - var(pos2)) < 1Mb #var返回的是起始和终止的中部区域
show = no
</rule>
<rule>
condition = var(intrachr)
condition = abs(var(pos1) - var(pos2)) < 50Mb
bezier_radius = 1.1r #贝塞尔半径大于1 弧度偏向外
bezier_radius_purity = 0.25
color = dblue_a2
</rule>
<rule>
condition = var(intrachr)
condition = abs(var(pos1) - var(pos2)) < 100Mb
bezier_radius = 0.75r
crest = 2
bezier_radius_purity = 0.25
color = dgreen_a2
</rule>
<rule>
condition = var(intrachr)
show = no
</rule>
<rule>
condition = var(interchr)
condition = on(hs2,65Mb,75Mb)
z = 60
color = red_a1
thickness = 5
radius = 1r
</rule>
<rule>
condition = max(var(size1),var(size2)) > 40kb
z = 50
color = black
thickness = 5
</rule>
<rule>
condition = max(var(size1),var(size2)) > 10kb
z = 45
color = dgrey
thickness = 4
</rule>
<rule>
condition = max(var(size1),var(size2)) > 5kb
z = 40
color = grey
thickness = 3
</rule>
<rule>
condition = max(var(size1),var(size2)) > 1kb
z = 35
color = lgrey
thickness = 2
</rule>
</rules>
</link>
</links>
···