熊朝亮
如何将多个Bed文件进行合并?
2014-11-29 23:54
阅读:9415
假如有1.bed, 2.bed, 3.bed,如何对这三个bed文件进行合并,去掉冗余的条目,得到一个并集all.bed?

答案是bedtools的mergeBed
[noncode@gnode105 script]$ mergeBed -h
Tool: bedtools merge (aka mergeBed)
Version: v2.16.2
Summary: Merges overlapping BED/GFF/VCF entries into a single interval.

Usage: bedtools merge [OPTIONS] -i <bed/gff/vcf>
Options:
-s Force strandedness. That is, only merge features
that are the same strand.
- By default, merging is done without respect to strand.

-n Report the number of BED entries that were merged.
- Note: "1" is reported if no merging occurred.

-d Maximum distance between features allowed for features
to be merged.
- Def. 0. That is, overlapping & book-ended features are merged.
- (INTEGER)

-nms Report the names of the merged features separated by semicolons.

-scores Report the scores of the merged features. Specify one of
the following options for reporting scores:
sum, min, max,
mean, median, mode, antimode,
collapse (i.e., print a semicolon-separated list),
- (INTEGER)

Notes:
(1) All output, regardless of input type (e.g., GFF or VCF)
will in BED format with zero-based starts
(2) The input file (-i) file must be sorted by chrom, then start.[/mw_shl_code]


转载本文请联系原作者获取授权,同时请注明本文来自熊朝亮科学网博客。

链接地址:https://wap.sciencenet.cn/blog-1509670-847318.html?mobile=1

收藏

分享到:

当前推荐数:0
推荐到博客首页
网友评论0 条评论
确定删除指定的回复吗?
确定删除本博文吗?