Simple circos code
According to the tutorials (http://circos.ca/documentation/tutorials/), to draw a graph using circos, at least three files needed.
Download the software and install: http://circos.ca/software/download/
running code: /PATH/circos/bin/circos -conf z1.conf
Files preparation:
# 1. genome file#
# can be in one file or one species one file. Following is the format.
chr - sb7 sb7 0 354 chr7
chr - si7 si7 0 322 chr7
chr - sb8 sb8 0 215 chr8
chr - si8 si8 0 251 chr8
chr - sb15 sb15 0 130 chr15
chr - si15 si15 0 147 chr15
# 2. hist.txt #
# this is not necessary to have.
"V1" "temp1" "temp"
"1" 1 0.5 555
"2" 1 1.5 157
"3" 1 2.5 0
"4" 1 3.5 3
"5" 1 4.5 3
"6" 1 5.5 0
"7" 1 6.5 0
"8" 1 7.5 0
"9" 1 8.5 0
"10" 1 9.5 0
"11" 1 10.5 0
"12" 1 11.5 0
"13" 1 12.5 0
"14" 1 13.5 0
"15" 1 14.5 0
"16" 1 15.5 0
"17" 1 16.5 0
"18" 1 17.5 0
"19" 1 18.5 0
"20" 1 19.5 0
"21" 1 20.5 0
"22" 1 21.5 8
"23" 1 22.5 209
"24" 2 0.5 0
"25" 2 1.5 1
"26" 2 2.5 0
"27" 2 3.5 0
"28" 2 4.5 0
"29" 2 5.5 0
"30" 2 6.5 0
# 3. link file#
linkzm-si-1-1 zm1 0 61 z=0,color=chr1
linkzm-si-1-1 si1 0 78 z=0,color=chr1
linkzm-si-1-9 zm1 71 125 z=0,color=chr9
linkzm-si-1-9 si9 0 73 z=0,color=chr9
linkzm-si-1-14 zm1 125 222 z=0,color=chr14
linkzm-si-1-14 si14 0 73 z=0,color=chr14
# 4.genome.conf#
#karyotype = /home/xlai/Desktop/UntitledFolder/motif/Newcluster/circos/zm.txt,/home/xlai/Desktop/UntitledFolder/motif/Newcluster/circos/sb.txt,/home/xlai/Desktop/UntitledFolder/motif/Newcluster/circos/si.txt
karyotype = /home/xlai/Desktop/UntitledFolder/motif/Newcluster/circos/sbsi.scaled.txt
#chromosomes_units = 100
#this histgram part is not necessary
########################################histogram
<plots>
<plot>
type = histogram
file = /home/xlai/Desktop/UntitledFolder/motif/Newcluster/circos/sbsi.hist.txt
fill_color = vdgrey
color = black
stroke_thickness = 3
r1 = 0.99r
r0 = 0.7r
orientation = out
extend_bin = yes
<rule>
condition = on()
show = yes
</rule>
</plot>
<plot>
type = histogram
file = /home/xlai/Desktop/UntitledFolder/motif/Newcluster/circos/clock.txt
fill_color = vdgrey
color = black
stroke_thickness = 1
r1 = 0.99r
r0 = 0.95r
orientation = in
extend_bin = yes
<rule>
condition = on()
show = yes
</rule>
</plot>
</plots>
############################z1.conf
<ideogram>
<spacing>
default = 0.001r
</spacing>
# Ideogram position, fill and outline
radius = 0.9r
thickness = 20p
fill = yes
stroke_color = dgrey
stroke_thickness = 2p
# Minimum definition for ideogram labels.
show_label = yes
# see etc/fonts.conf for list of font names
label_font = default
label_radius = 1r + 75p
label_size = 40
label_parallel = yes
</ideogram>
###############
show_links = yes
### links
<links>
show = conf(show_links)
ribbon = yes
flat = yes
radius = 0.7r-30p
bezier_radius = 0.3r
color = black_a5
<link>
file = /home/xlai/Desktop/UntitledFolder/motif/Newcluster/circos/lk.sbsi.scaled.txt
<rules>
<rule>
condition = var(size1) < 0
show = no
</rule>
<rule>
condition = 1
color = eval(sprintf("chr%s_a3",substr(var(chr2),2)))
z = eval(remap_int(var(size2),0,1,0,-100))
</rule>
</rules>
</link>
</links>
################################################################
<image>
# Included from Circos distribution.
<<include etc/image.conf>>
</image>
# RGB/HSV color definitions, color lists, location of fonts, fill patterns.
# Included from Circos distribution.
<<include etc/colors_fonts_patterns.conf>>
# Debugging, I/O an dother system parameters
# Included from Circos distribution.
<<include etc/housekeeping.conf>>
Simple circos code的更多相关文章
- simple python code when @ simplnano
code: import serial,time,itertools try: ser=serial.Serial(2,115200,timeout=0) except: print 'Open CO ...
- Transport Block Size, Throughput and Code rate-----http://www.simpletechpost.com/2012/12/transport-block-size-code-rate-protocol.html
Transport Block Size, Throughput and Code rate Since the size of transport block is not fixed, oft ...
- Peer Code Reviews Made Easy with Eclipse Plug-In
欢迎关注我的社交账号: 博客园地址: http://www.cnblogs.com/jiangxinnju/p/4781259.html GitHub地址: https://github.com/ji ...
- GOOGLE CODE ANDROID 开源项目 集合
转:http://blog.csdn.net/dellheng/article/details/7163333 1. ZXing http://code.google.com/p/zx ...
- C# How To Read .xlsx Excel File With 3 Lines Of Code
Download Excel.zip - 9.7 KB Download ExcelDLL.zip - 3.7 KB Introduction We produce professional busi ...
- Coroutine in Java - Quasar Fiber实现--转载
转自 https://segmentfault.com/a/1190000006079389?from=groupmessage&isappinstalled=0 简介 说到协程(Corout ...
- Unit Testing with NSubstitute
These are the contents of my training session about unit testing, and also have some introductions a ...
- Atitit webservice的发现机制 discover机制
Atitit webservice的发现机制 discover机制 1.1. Ws disconvert 的组播地址和端口就是37021 1.2. Ws disconvert的发现机制建立在udp组播 ...
- Atitti onvif 设备发现与原理
Atitti onvif 设备发现与原理 1.1. ,有以下几个步骤:1 1.2. 设备搜索原理及编程技巧:2 1.3. Ws disconvert 的组播地址和端口就是37022 1)发现ipca ...
随机推荐
- Python OS模块常用
python 读写.创建 文件 第二个:目录操作-增删改查 第三个:判断 第四个:PATH 第四个:os.mknod 创建文件(不是目录) import os os.chdir("/&quo ...
- AS3.0 给addEventListener里的方法加上参数传递
方法一:for(var i:int=1;i<=4;i++){this["btn"+i].addEventListener(MouseEvent.CLICK,EventUp(b ...
- 苹果手机的SB系列(2)为什么不能重命名?
为什么没有重命名? 在手机端不能重命名,在WINOWS端文件是只读的,连他TM的只读属性都无法改,不能重命名,你让我怎么备份? 我怎么知道哪些照片上次备份过了?又重头来过?还是要用苹果的MAC?这种态 ...
- 二叉树之AVL树
高度为 h 的 AVL 树,节点数 N 最多2^h − 1: 最少N(h)=N(h− 1) +N(h− 2) + 1. 最少节点数n 如以斐波那契数列可以用数学归纳法证明: 即: N(0) = 0 ( ...
- Webpack学习-Loader
什么是Loader? 继上两篇文章webpack工作原理介绍(上篇.下篇),我们了解到Loader:模块转换器,也就是将模块的内容按照需求装换成新内容,而且每个Loader的职责都是单一,只会完成一种 ...
- error connecting to master 'repl@192.168.1.107:3306' - retry-time: 60 retries: 86400
很多时候,就算thread 正在进行,也不代表没有错误,一定要看看具体表示错误的变量 mysql> show slave status\G *************************** ...
- 常见类 Object
1.在Java类继承结构中Java.lang.Object类位于顶端 2.如果定义一个Object类没有使用extends关键字声明其父类,则其父类为Java.lang.Object类 3.Objec ...
- CTF资料
题目大概有这么几个 web,密码学,pwn(综合渗透),misc(杂项),reverse(逆向),ppc(编程类) 入门: http://bobao.360.cn/ctf/ 360的CTF训练营ww ...
- Collections集合工具类
一.Collection与Collections Collection 是所有单列集合的根接口 Collection 是操作集合的工具类 二.Collections中常见的方法:(大都是static方 ...
- mp4文件解析(纯c解析代码)
参考链接:1. mp4文件格式解析 https://www.cnblogs.com/ranson7zop/p/7889272.html 2. MP4文件格式分析及分割实现(附源码) https: ...