ChIP-Seq,MeRIP-seq峰(peak),eccDNA等染色体分布可视化
人类基因组由1-22、X、Y等染色体构成,染色体经过细胞学处理后会呈现出深浅不同的染色带。 染色带的数目、部位、宽窄和着色深浅均相对稳定性,所以每一条染色体都有固定的分带模式,即称带型。 染色体带型是鉴别染色体的重要依据。
图1. 人类染色体
通常我们可以将ChIP-seq、MeRIP-seq、eccDNA等高通量测序的分析结果简化为染色体上的一些区域,即chr:start-end。 将获得的结果在染色体上进行可视化不仅能够看出感兴趣区域的染色体分布和密度,而且能够展示多个条件下差异情况,非常形象。 常见的展示方式有条形和圆形(circos)。 今天我们来看下条形展示方式。
1,打开染色体分布绘图页面
首先,使用浏览器(推荐chrome或者edge)打开染色体分布绘图页面。左侧为常见作图导航,中间为数据输入框和可选参数,右侧为描述和结果示例。也可以在主页搜索框中搜索chromosome,找到绘图页面。
https://www.bioinformatics.com.cn/plot_basic_up_down_gene_chromsome_distribution_plot_095
图2.染色体分布绘图页面
2,示例数据
点击右侧“示例数据”链接下载excel格式的示例数据。
图3. 输入数据示例
示例数据(仅供参考)包括两组,即最多绘制染色体上面和下面两组数据:
第1列是染色体号,去掉了chr
第2列是基因组坐标,对于小片段来说,起始坐标,终止坐标及中点坐标在染色体尺度上基本没有区别,因此,这里可以直接使用获得的peak的中点坐标。
3,粘贴示例数据
直接复制示例数据中的up两列数据,和down两列数据(可选),分别粘贴到输入框。
注意:不是拷贝excel文件,是拷贝excel文件里边的数据。另外粘贴到输入框后,格式乱了没关系,只要在excel中是整齐的就行。并且数据矩阵中不能有空的单元格,中文字符等。
图4.输入
4,修改参数,并提交
我们设置了图片尺寸,颜色,物种,字体等参数。
图5.可调参数
5,提交出图
图6.预览与下载
图中红色为数据1,例如上调peak;蓝色为数据2,例如下调peak。
没有预览就是没有出图,这时请参考示例数据,检查自己输入数据的格式。
遇到文字截断,需要修改字体、调整字体大小等,。
微生信助力发文章,谷歌引用630+,知网引用470+
ChIP-Seq,MeRIP-seq峰(peak),eccDNA等染色体分布可视化相关推荐
- seq to seq
NLP(3)--seq to seq 什么是Seq2Seq网络? 在Seq2Seq模型中采用了这种 Encoder-Decoder架构,其中 Encoder 是一个RNNCell(RNN ,GRU,L ...
- NLP(3)| seq to seq 模型
NLP(1) | 词向量one hot编码词向量编码思想 NLP(2) | 中文分词分词的概念分词方法分类CRFHMM分词 什么是Seq2Seq网络? 在Seq2Seq模型中采用了这种 Encoder ...
- matlab如何表示三峰正态分布,正态分布及常用分布的matlab编程实现
<正态分布及常用分布的matlab编程实现>由会员分享,可在线阅读,更多相关<正态分布及常用分布的matlab编程实现(3页珍藏版)>请在人人文库网上搜索. 1.functio ...
- shell实例第8讲:seq命令
shell脚本之seq命令 seq 用于生成从一个数到另一个数之间的所有整数. (1)seq 5或seq 1 5 (2)-f 指定输出格式,seq -f "%3g"1 3 数字位 ...
- R语言seqm_R语言seq()函数用法
1.seq() 用来生成一组数字的函数. Usage: ## Default S3 method: seq(from = 1, to = 1, by = ((to - from)/(length.ou ...
- 使用Seq搭建免费的日志服务
Seq简介 Seq是老外开发的一个针对.NET平台非常友好的日志服务.支持容器部署,提供一个单用户免费的开发版本. 官网:https://datalust.co/seq 使用文档:https://do ...
- EPICS记录参考--Sequence记录(seq)
Sequence记录用于触发最多16个其它记录运行并且发送这些值到那些记录.除了它将获取一个输入值并且写一个输出值而不是仅仅运行一个forward链接的集合外,它类似于fanout记录.它也能够指定若 ...
- TCP协议中的Ack和Seq号
TCP协议中的Ack和Seq号 一.基本概念 seq:表示本次发送数据的偏移量,也就是从哪里开始发送数据. len:表示本次tcp携带的数据长度. ack:首先意味着已经收到对方多少字节数据,其次告诉 ...
- R语言中seq函数的用法
seq(from,to,length), 该函数的意思是生成一组数字,从from开始,到to结束,每两个数间的间隔是length,如 seq(2,10,2),会生成一组数:2 4 6 8 10 seq ...
- scala的Seq与java的List相互转换
在使用sparksql的dataset时有些地方需要将seq和list进行互相转化使用,scala集合提供相关的转化操作: //Seq 转 List List<Column> list = ...
最新文章
- Canvas基础知识
- Oracle忘记用户名和密码的解决方案
- come type6 定义_COMe-B6101龙芯3A3000 COM Express Type6模块
- Spring之RestTemplate如何返回map结果集
- 入门 | 初学者必读:解读14个深度学习关键词
- 影视感悟专题---1、B站-魔兽世界代理及其它乱七八糟
- Jmeter 监控多台服务器CPU、内存、i/o等资源
- Git笔记(12) 分支使用
- RS 纠删码为什么可以提高分布式存储可靠性?| 原力计划
- 拓端tecdat|R语言圆填充( Circle packing)算法圆堆图圆形空间填充算法可视化
- 多线程(what,why,when)
- JAX-RS之Jersey入门
- 主机宝iis版_主机宝iis版下载|
- SPSS基础教程:统计分析前的准备
- 电力系统104规约帧报文解析
- C#配置文件操作】程序配置文件App.Config操作
- java学习笔记(10) 第10 章(上) - 面向对象编程(高级部分)
- Linux oracle 怎么安装,oracle在linux下怎么安装
- c语言scandf用法,scan的用法总结大全
- Html转义字符表详细记录