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简介及应用的更多相关文章

  1. ASP.NET Core 1.1 简介

    ASP.NET Core 1.1 于2016年11月16日发布.这个版本包括许多伟大的新功能以及许多错误修复和一般的增强.这个版本包含了多个新的中间件组件.针对Windows的WebListener服 ...

  2. MVVM模式和在WPF中的实现（一）MVVM模式简介

    MVVM模式解析和在WPF中的实现(一) MVVM模式简介 系列目录: MVVM模式解析和在WPF中的实现(一)MVVM模式简介 MVVM模式解析和在WPF中的实现(二)数据绑定 MVVM模式解析和在 ...

  3. Cassandra简介

    在前面的一篇文章<图形数据库Neo4J简介>中,我们介绍了一种非常流行的图形数据库Neo4J的使用方法.而在本文中,我们将对另外一种类型的NoSQL数据库——Cassandra进行简单地介 ...

  4. REST简介

    一说到REST,我想大家的第一反应就是“啊,就是那种前后台通信方式.”但是在要求详细讲述它所提出的各个约束,以及如何开始搭建REST服务时,却很少有人能够清晰地说出它到底是什么,需要遵守什么样的准则. ...

  5. Microservice架构模式简介

    在2014年,Sam Newman,Martin Fowler在ThoughtWorks的一位同事,出版了一本新书<Building Microservices>.该书描述了如何按照Mic ...

  6. const,static,extern 简介

    const,static,extern 简介 一.const与宏的区别: const简介:之前常用的字符串常量,一般是抽成宏,但是苹果不推荐我们抽成宏,推荐我们使用const常量. 执行时刻:宏是预编 ...

  7. HTTPS简介

    一.简单总结 1.HTTPS概念总结 HTTPS 就是对HTTP进行了TLS或SSL加密. 应用层的HTTP协议通过传输层的TCP协议来传输,HTTPS 在 HTTP和 TCP中间加了一层TLS/SS ...

  8. 【Machine Learning】机器学习及其基础概念简介

    机器学习及其基础概念简介 作者:白宁超 2016年12月23日21:24:51 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现的深入理解.本系列文章是作者结 ...

  9. Cesium简介以及离线部署运行

    Cesium简介 cesium是国外一个基于JavaScript编写的使用WebGL的地图引擎,一款开源3DGIS的js库.cesium支持3D,2D,2.5D形式的地图展示,可以自行绘制图形,高亮区 ...

随机推荐

  1. 关于IOS给我的启发

    用了将近一年半的iOS,从4到4S,iOS5到iOS6.这里谈谈自己对iOS的一些看法,以及这款移动操作系统给我的启发.我知道这个帖子发出来可能有点“危险”.我从不发水贴,这些积分都是大家给的,不是灌 ...

  2. JDBC进行简单的增删改查

    一.准备工作(一):MySQL安装配置和基础学习 二.准备工作(二):下载数据库对应的jar包并导入 三.JDBC基本操作 (1)定义记录的类(可选) (2)连接的获取 (3)insert (4)up ...

  3. flutter初探

    这两天看了下flutter,感觉这两年可能会爆发,所以尝试在mac和win10上面跑了下hello world... 移动技术简介 原生开发 跨平台技术简介 H5+原生(Cordova.Ionic.微 ...

  4. 发送邮件(单独文字)的方法(网易邮箱 OR QQ邮箱)

    # coding:utf-8import smtplibfrom email.mime.text import MIMEText # 发邮件相关的参数# 网易邮箱用这个# smtpserver=&qu ...

  5. R语言学习——R读取txt、csv、xls和xlsx格式文件

    最近项目中运用到了R读取文件数据,所以把相关好用的.经过验证的方法总结了一下,有效避免下次入坑. 1. R读取txt文件 使用R读取txt文件直接使用read.table()方法进行读取即可,不需要加 ...

  6. 微信小程序 GMT+0800 (中国标准时间) WXSS 文件编译错误

    请尝试在控制台输入openVendor() ,清除里面的wcsc wcsc.exe 然后重启工具

  7. CA单向认证和双向认证的区别？

    1:单向认证,内容会被串改吗?

  8. ossim 4.1 安装

    在这说明一下,本人也是第一次安装这东东很多我自己也不懂,只能把自己的安装步骤通过贴图呈现出来了!!!! 我用root用户登录的,顺便查看ip 在浏览器打开网址https://192.168.86.15 ...

  9. solr 7+tomcat 8 + mysql实现solr 7基本使用(安装、集成中文分词器、定时同步数据库数据以及项目集成)

    基本说明 Solr是一个开源项目,基于Lucene的搜索服务器,一般用于高级的搜索功能: solr还支持各种插件(如中文分词器等),便于做多样化功能的集成: 提供页面操作,查看日志和配置信息,功能全面 ...

  10. JavaScript中的闭包与匿名函数

    知识内容: 1.预备知识 - 函数表达式 2.匿名函数 3.闭包 一.函数表达式 1.定义函数的两种方式 函数声明: 1 function func(arg0, arg1, arg2){ 2 // 函 ...