bedtools简介及应用
1)背景
处理基因组数据中,比较基因组不同区域,例如寻找overlap等,是一种基本的且常见的问题。虽然UCSC 中‘Table Browser’或者Galaxy可以用来处理,但是当这些工具面对大的数据的时候就会显得力不从心。因此,需要一款快速、灵活的软件来批量处理数据集。
bedtools是一款用C++编写的小巧且灵活的软件来处理这些复杂的问题,可以用来比较、操作、注释bed和gff文件中的genomic features。它设计主要是在linux环境下,可以和awk、grep、sort 等实现无缝对接,且支持bam格式文件,可以将bam文件转换为bed文件,因此对于大数据处理来说非常方便。
由于在版本在不断的更新,我们将依据最初版本的参考文献给出每个示例图,让大家知道理解每个命令的含义,知道如何去学习一个软件。明白了原理,具体要到的时候,依据不同版本的帮助文档就很容易上手。
2)功能
下图是bedtool的可用操作(注意版本不同,命令可能会变,但大致功能不会改变,这里主要是展示其原始文献):虽然有些命令已经该换名字,例如fastaFromBed在后续的版本中的别名是getfasta,但是这些命令已经可以被执行,尽管后续版本不断有新的参数被加入。(在一再次注意,这里主要通过图解来说每个命令到底在做什么。不讲实际应用)
2.1)intersectBed
‘intersecting’ or ‘overlapping’指feature之间至少有1bp 的重叠:
基本用法:intersectBed [OPTIONS] [-a <BED> || -abam <BAM>] -b <BED>
2.1.1)找出overlap并输出
intersectBed –a A.bed –b B.bed
2.1.2)当有overlap输出原始的 “A” feature
intersectBed –a A.bed –b B.bed -wa
2.2)pairToBed
主要用于将BedPE feature或者paired-end BAM alignment 与另一个bed文件比较寻找overlap
基本用法:pairToBed [OPTIONS] [-a <BEDPE> || -abam <BAM>] -b <BED>
2.2.1)默认情况下,如果任何一端与B有交集则输出A(Report A if either end overlaps B.)
2.2.2)-type 参数可以用来控制输出
-type both: 当用该参数的时候,只有在双端都和B有交集的情况下才输出A(Report A only if both ends overlap B)
-type neither: 如果双端和B都没有交集则输出A(Report A only if neither end overlaps B.)
-type ispan: 输出A,如果A双端间区域和B有重叠(Report A if it’s “inner span” overlaps B.)
2.2.3) -f参数通过调节最少比对长度来控制输出
输出A,如果一端至少有50%区域与B重叠(eport A only at least 50% of one of the two ends is overlapped by B)
pairToBed -a A.bedpe -b B.bed -f 0.5
2.3) pairToPair
寻找两个BEDPE文件之间的重叠
基本用法:pairToPair [OPTIONS] -a <BEDPE> -b <BEDPE>
2.3.1)默认情况下输出A,如果A和B都有链信息,且A的双端和B有交集
2.4)bamToBed
用来将bam格式转化为Bed格式
基本用法:bamToBed [OPTIONS] -i <BAM> ,默认输出的bed 文件有6列
bamToBed -i RNA-seq.bam |head -6
2.5)windowBed
在用户指定的窗口内返回两个BED文件之间的overlap
基本用法:windowBed [OPTIONS] -a <BED> -b <BED> ,默认为上下游1M窗口
2.5.1)默认为上下游1M,例如:
2.5.2)可以通过(-l 和-r)参数来控制上下游,加入不同的window
2.6)closestBed
如果A和B之间没有overlap,则输出与A最近的B的信息
2.7)subtractBed
用处:主要输用于移除重叠部分及间隔区
基本用法:subtractBed [OPTIONS] -a <BED> -b <BED>
默认情况下
2.8)mergeBed
主要用于将重叠的feature合并
基本用法:mergeBed [OPTIONS] -i <BED>
例如:
、
2.9)coverageBed
主要是统计一个BED文件中要素相对于另一个的覆盖深度和广度
基本用法:coverageBed [OPTIONS] -a <BED> -b <BED>
例如:统计结果解释可以参考下图。
2.10) genomeCoverageBed
对特征文件的覆盖情况进行统计
基本用法:genomeCoverageBed [OPTIONS] -i <BED> -g <GENOME>
2.11)fastaFromBed
用于提取fasta文件中,与bed文件中匹配的序列,形成一个新的fasta文件
基本用法:fastaFromBed [OPTIONS] -fi <input FASTA> -bed <BED> -fo <output FASTA>
2.12)maskFastaFromBed
用于将bed 中指定的位置在fasta中进行mask
基本用法:maskFastaFromBed [OPTIONS] -fi <input FASTA> -bed <BED> -fo <output FASTA>
默认情况下是硬soft,而不是软soft
2.13) shuffleBed
用于置换feture在染色体中的location
shuffleBed [OPTIONS] -i <BED> -g <GENOME>
2.14)sortBed
对bed文件进行排序
基本用法:sortBed [OPTIONS] -i <BED>
默认情况下,对bed文件按照染色体,然后按照起始位置。
2.15)linksBed
为bed 文件建立html
基本用法:linksBed [OPTIONS] -i <BED> > <HTML file>
例如:
2.16)complementBed
输出在bed 文件中不包含的区域
基本用法:complementBed [OPTIONS] -i <BED> -g <GENOME>
例如
bedtools简介及应用的更多相关文章
- ASP.NET Core 1.1 简介
ASP.NET Core 1.1 于2016年11月16日发布.这个版本包括许多伟大的新功能以及许多错误修复和一般的增强.这个版本包含了多个新的中间件组件.针对Windows的WebListener服 ...
- MVVM模式和在WPF中的实现(一)MVVM模式简介
MVVM模式解析和在WPF中的实现(一) MVVM模式简介 系列目录: MVVM模式解析和在WPF中的实现(一)MVVM模式简介 MVVM模式解析和在WPF中的实现(二)数据绑定 MVVM模式解析和在 ...
- Cassandra简介
在前面的一篇文章<图形数据库Neo4J简介>中,我们介绍了一种非常流行的图形数据库Neo4J的使用方法.而在本文中,我们将对另外一种类型的NoSQL数据库——Cassandra进行简单地介 ...
- REST简介
一说到REST,我想大家的第一反应就是“啊,就是那种前后台通信方式.”但是在要求详细讲述它所提出的各个约束,以及如何开始搭建REST服务时,却很少有人能够清晰地说出它到底是什么,需要遵守什么样的准则. ...
- Microservice架构模式简介
在2014年,Sam Newman,Martin Fowler在ThoughtWorks的一位同事,出版了一本新书<Building Microservices>.该书描述了如何按照Mic ...
- const,static,extern 简介
const,static,extern 简介 一.const与宏的区别: const简介:之前常用的字符串常量,一般是抽成宏,但是苹果不推荐我们抽成宏,推荐我们使用const常量. 执行时刻:宏是预编 ...
- HTTPS简介
一.简单总结 1.HTTPS概念总结 HTTPS 就是对HTTP进行了TLS或SSL加密. 应用层的HTTP协议通过传输层的TCP协议来传输,HTTPS 在 HTTP和 TCP中间加了一层TLS/SS ...
- 【Machine Learning】机器学习及其基础概念简介
机器学习及其基础概念简介 作者:白宁超 2016年12月23日21:24:51 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现的深入理解.本系列文章是作者结 ...
- Cesium简介以及离线部署运行
Cesium简介 cesium是国外一个基于JavaScript编写的使用WebGL的地图引擎,一款开源3DGIS的js库.cesium支持3D,2D,2.5D形式的地图展示,可以自行绘制图形,高亮区 ...
随机推荐
- ng-model的用法
参考: http://www.cnblogs.com/guanglin/p/5200097.html http://www.runoob.com/angularjs/ng-ng-cloak.html ...
- You Arent Gonna Need It 理解
软件开发是一项无法驾驭全局的工作,其工作效率,工作精神都高度依赖需求反馈. 开发不存在需求的功能,实际上是除了要开发软件,还要构思需求,这样一正一反,严重消耗了脑力资源. 如果是自己开发软件,需求来自 ...
- jQuery 的noConflict()的使用.
我们项目现在需要用到两个js库.一个是jQuery库,还有一个是我们自己开发的轻量级的gys.js库. 而gys库对外提供的接口也是$符号.和jQuery库是一样的,这个时候,两个库就会发生冲突了,我 ...
- 探究CSS中border-top属性的使用
上一节我们介绍了CSS top属性的用法,那么这节关于CSS border-top属性用法学习起来就会轻松一些,border-top 简写属性把上边框的所有属性设置到一个声明中. 本文向大家描述一下C ...
- 学习MongoDB 八: MongoDB索引(索引限制条件)(二)
一.简介 我们上一篇介绍了索引基本操作,通过db.collection.createIndex(keys, options)语法创建索引,我们继续介绍地理空间索引.索引的限制,使我们在MongoDB时 ...
- python入门第0篇 Windows下python的安装及pip安装和使用
知识内容: 1. python2和python3安装 2. pip安装及pip命令使用 注:安装python3就可以了,python2除非项目开发需要否则不用安装,目前学习python就使用pytho ...
- vim编程设置
在终端下使用vim进行编辑时,默认情况下,编辑的界面上是没有显示行号.语法高亮度显示.智能缩进 等功能的.为了更好的在vim下进行工作,需要手动设置一个配置文件:.vimrc.在启动vim时,当前用户 ...
- c#面向对象基础4
一.namespace 命名空间 作用:解决不同类重名的问题 我们可以认为类是属于命名空间的 当我们需要再一个类中与另一个类建立关系时,通过命名空间来区别不同的类.所以需要我们这样做:导入命名空间 ...
- leetcode405
public class Solution { string ConverToHex(string bit) { var s = ""; switch (bit) { " ...
- js实现点击按钮弹出上传文件的窗口
转自:https://www.jb51.net/article/100916.htm 1.详细描述 在页面上设置一个“选择文件”按钮,点击该按钮,会弹出本地磁盘信息用于选择文件. 2.代码 ? 1 2 ...