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的更多相关文章

  1. simple python code when @ simplnano

    code: import serial,time,itertools try: ser=serial.Serial(2,115200,timeout=0) except: print 'Open CO ...

  2. 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 ...

  3. Peer Code Reviews Made Easy with Eclipse Plug-In

    欢迎关注我的社交账号: 博客园地址: http://www.cnblogs.com/jiangxinnju/p/4781259.html GitHub地址: https://github.com/ji ...

  4. GOOGLE CODE ANDROID 开源项目 集合

    转:http://blog.csdn.net/dellheng/article/details/7163333 1.        ZXing  http://code.google.com/p/zx ...

  5. 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 ...

  6. Coroutine in Java - Quasar Fiber实现--转载

    转自 https://segmentfault.com/a/1190000006079389?from=groupmessage&isappinstalled=0 简介 说到协程(Corout ...

  7. Unit Testing with NSubstitute

    These are the contents of my training session about unit testing, and also have some introductions a ...

  8. Atitit webservice的发现机制 discover机制

    Atitit webservice的发现机制 discover机制 1.1. Ws disconvert 的组播地址和端口就是37021 1.2. Ws disconvert的发现机制建立在udp组播 ...

  9. Atitti  onvif 设备发现与原理

    Atitti  onvif 设备发现与原理 1.1. ,有以下几个步骤:1 1.2. 设备搜索原理及编程技巧:2 1.3. Ws disconvert 的组播地址和端口就是37022 1)发现ipca ...

随机推荐

  1. Python3+Selenium获取session和token供Requests使用教程

    一.背景说明 之前写了一款简单的api模糊测试工具,之前系统可以使用http Base认证现在改成session形式并加上了token. 最简单的改造方法,是自己先在浏览器手动登录,然后提取出sess ...

  2. linux常用命令 运算命令

    linux的运算命令 expr命令,对整数进行运算 > expr的运算必须用空格间隔开 > \* 表示转义字符 > 保持先乘除后加减,如果需要优先运算则需要加命令替换符 > 也 ...

  3. Win10系列:C#应用控件进阶6

    路径 路径(Path)可以用来定义任意形状的曲线和几何图形,当然这种任意性也带来了复杂性.为了方便的绘制几何图形,微软在Visual Studio 2012安装包中为程序开发者提供了免费的Blend ...

  4. View事件体系

    View事件体系 文章目录 View事件体系 一.Android View基础知识 1.1 View简介 1.2 View分类 1.3 View的结构 1.4 View的坐标 1.4.1 Androi ...

  5. c语言五子棋

    #include <stdio.h>#include <stdlib.h>#include <windows.h>#include <conio.h> ...

  6. 常用的lamp环境以及一些依赖包的安装

  7. 安装SQl Server 报错 "需要 Microsoft.NET Framework 3.5 ServicePack 1" 解决方法

    前言 之前装Sql Server都没遇到过这样的问题, 昨天重装了系统之后, 然后安装SQl Server 报错,提示 "需要 Microsoft.NET Framework 3.5 Ser ...

  8. java完整并发知识结构图

    一张大的java并发知识结构图,梳理清楚知识的脉络,知识不再零散

  9. js事件、事件流以及target、currentTarget、this那些事

    你是如此简单我却将你给遗忘   前面面试被问到js的事件机制  target.currentTarget.碰巧今天有时间来拔一拔,顺便记下.

  10. css实现文本块在容器中垂直居中

    效果如图: css代码: .textContainer { display: -webkit-box; -webkit-box-orient: vertical; -webkit-box-pack: ...