美文网首页绘图
circos画图中links文件生成

circos画图中links文件生成

作者: 徒唤奈何_c5f0 | 来源:发表于2023-03-30 14:20 被阅读0次

文件格式要求


chr-a    start-a    end-a    chr-b    start-b    end-b options
chr9    4210121 4211225 chr10   19666298    19668816    color=140,105,20
chr9    4221308 4222192 chr10   19647682    19648089    color=140,105,20
chr9    4225247 4230980 chr10   19439605    19445179    color=140,105,20
chr9    20249400    20252312    chr11   798933  803498  color=140,115,20
chr9    20307676    20311487    chr11   792810  795351  color=140,115,20
##文件共有7列,使用tab分隔;每一行表示一对有联系的区域,前3列和后3列分别定义一个区域,最后一列为颜色

如何生成该格式?

  • 生成共线性文件

    NTX_1g01972 1972    NTX_1g01973 1973    1
    NTX_1g01980 1980    NTX_1g01979 1979    1
    NTX_1g01982 1982    NTX_1g01983 1983    -1
    NTX_1g01988 1988    NTX_1g01987 1987    1
    NTX_1g02003 2003    NTX_1g02002 2002    1
    NTX_1g02004 2004    NTX_1g02003 2003    1
    ##这是共线性文件(在wgdi的共线性分析中生成)
    
    1   NTX_1g00001 12237   16451   -   1   chr1.g00002.1
    1   NTX_1g00002 58474   59181   +   2   chr1.g00003.1
    1   NTX_1g00003 59836   67849   +   3   chr1.g00004.1
    1   NTX_1g00004 67865   69721   +   4   chr1.g00005.1
    ##这是wgdi生成的特殊比对gff文件
    
    awk -F "\t" '{print $1 "\t" $3}' 002.con | awk '/^NTX/{print $1}' | xargs -I {} grep {} ./NTX.gff | awk '{a=substr($2,1,8)}{print a "\t" $3 "\t" $4}' > 1.tmp
    awk -F "\t" '{print $1 "\t" $3}' 002.con | awk '/^NTX/{print $2}' | xargs -I {} grep {} ./NTX.gff | awk '{a=substr($2,1,8)}{print a "\t" $3 "\t" $4}' > 2.tmp
    paste 1.tmp 2.tmp >link.tmp
    ##把文件按照共线性文件生成新格式
    ##excel打开文件,去除掉染色体内的比对(if语句即可)
    
  • 给共线性文件添加颜色

##文件添加rgb颜色
##把文件输入到excel
#使color第二列为=MID(A1,4,2)*10+50,第四列为=MID(D3,4,2)*10等,可自行设置喜欢的颜色,之后到txt中去掉所有tab
#即可得到与两条染色体都相关的颜色代码,复制粘贴到links文件中,如下
chr1    4733065 4735993 chr18   13798474    13801492    color=60,185,20
chr1    4875770 4881299 chr18   14176762    14182099    color=60,185,20
chr1    24033076    24034330    chr2    25586015    25587734    color=60,25,20
chr1    24037473    24041511    chr2    25592131    25596678    color=60,25,20
#运行
circos -conf circos.conf

相关文章

  • circos 学习手册(十六)

    link 和 relationships(一) 1. 基本 links 的绘制 circos 核心用途之一就是展示...

  • circos作图四(links)

    通常circos的中间部分不是空白区域,会用一条条线进行连接,表示两个染色体部分区域有关系。对于link,circ...

  • 技巧 | BEDTools - genome file

    主题:如何简便地生成 bedtools, circos 等所需要的基因组染色体大小文件? 参考 Biostar |...

  • circos图

    一,安装circos 二,数据准备 1,计算染色体长度 2,生成染色体文件 7列 karyotype.txt 3,...

  • circos绘制彩色的links

    对于links而言,默认情况下所有的links 都是同一种颜色,通过color属性指定。当我们想要构建出彩色的li...

  • circos图

    绘制circos图一般需要上图中的三个文件。 一:染色体配置文件 主要是染色体的位置,名称,颜色的设置。其中第三列...

  • 彩票管理系统(C语言Linux)<3.文件的编译makefile

    首先来看下我的头文件和源文件: 我编写的makefile文件: 上图中: 绿色为生成的目标 黑色为生成目标所需要的...

  • 生成用于circos的基因密度文件的python脚本

    最近在做circos图,由于我的重复序列的gff文件不是标准格式,无法用软件生成,只好自己写一个python脚本,...

  • For circos模块---circos文件的准备

    circos软件在展示基因组信息方面具有十分重要的作用。关于该软件的使用可以在网上找到很多教程,但关于如何准备运行...

  • 如何使用Circos-4: Links and Relation

    说明,该系列原文写于2016年3月。 这一部分简单介绍links的基础配置,我们使用的文件依旧是 如何使用Circ...

网友评论

    本文标题:circos画图中links文件生成

    本文链接:https://www.haomeiwen.com/subject/gmfcddtx.html