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. Java读取txt文件和写入txt文件

    package com.nickwong.code; import java.io.*; /** * Created by Nickwong on 31/07/2018. * 根据1-8楼的建议,优化 ...

  2. 【技巧】EasyUI分页组件pagination显示项控制

    我们使用easyui分页时,分页组件显示项有很多,默认如下是这样的: 有时候我们并不想显示这么多控制项,我们可以通过属性来控制:如下来自EasyUI官网: 如下写法,在datagrid中使用,如下控制 ...

  3. react球形文字旋转标签

    /* * 球形文字旋转标签模块 * */ import React, {Component, PropTypes} from "react"; import ReactDOM fr ...

  4. python输入整数

    #!/usr/bin/env python#ecoding=utf-8'''Created on 2017年11月2日 @author: James zhan''' def fun(n): if n= ...

  5. Android 音视频深入 十六 FFmpeg 推流手机摄像头,实现直播 (附源码下载)

    源码地址https://github.com/979451341/RtmpCamera/tree/master 配置RMTP服务器,虽然之前说了,这里就直接粘贴过来吧 1.配置RTMP服务器 这个我不 ...

  6. iSlide——图标库、图示库的用法

    iSlide中,有一个“图示库”功能,主要功能是同时排列多块文字或多张图片.单击插图库,会弹出一个新的对话框.从中,可以选择权限.分类.数量数据和样式,也可以直接搜索.   下面就举一个例子:我要开一 ...

  7. C# foreach 值类型及引用类型迭代变量改变的方式

    C#中foreach不能改变迭代变量的值 然而此种说法只适用与值类型,更改值类型时会改变在栈上的内存分布 引用类型由于是引用地址的变更,不影响内存分布,所以能够在foreach中更改 至于引用类型中的 ...

  8. tinyproxy代理配置

    tinyproxy代理配置 应用场景: 生产机处于内网,无法直接访问外网,程序安装和漏洞修复等操作需要进行联网操作:通过在办公网(可访问外网)上设置代理服务器,生产机通过代理由办公网访问外网 代理服务 ...

  9. github下载项目代码到本地,不能运行 本地改完代码 再上传

    首先用git bash here,在指定目录下执行, git clone 将项目拉取下来, 试运行: 发现需要配置idea的SDK/jdk, 还要选择language level, 建立输出目录tar ...

  10. vue搭建

    ~感觉再不学点新技能,马上就要被淘汰了== 1.先安装node,直接点击下一步就OK  node官网   https://nodejs.org/en/ 2.安装完成之后,打开命令行工具,快捷键是cmd ...