multi-intersect-bed
Quick Start
- Command:
multiIntersectBed -i file1.bed file2.bed [file3.bed ...] [options] - Local executable:
/home/vimalinx/miniforge3/envs/bio/bin/multiIntersectBed - Full reference: See
references/help.md
When To Use This Tool
- Partition the genome into sub-intervals labeled by which of several files overlap them.
- Compare shared and unique regions across multiple BED / GFF / VCF inputs in one pass.
- Generate presence/absence matrices over segmented interval space.
- Include empty regions relative to a genome file with
-g -empty.
Common Patterns
# 1) Multi-file overlap segmentation
multiIntersectBed \
-i sample1.bed sample2.bed sample3.bed
# 2) Add a header with readable file names
multiIntersectBed \
-header \
-names tumor normal blacklist \
-i tumor.bed normal.bed blacklist.bed
# 3) Include empty regions across the genome
multiIntersectBed \
-i a.bed b.bed \
-g genome.txt \
-empty
Recommended Workflow
- Sort every input file by chromosome and start before running the tool.
- Decide whether you need just overlapping segmentation or also empty regions from a genome definition.
- Use
-namesand-headerwhen the output will be read by humans or imported into tables. - Post-process the membership columns to define how many files must support a region for your biological question.
Guardrails
- Each interval file must be sorted by chromosome and start.
- This tool segments coordinate space; the output intervals are often smaller than the original input intervals.
-emptyrequires-g.-incl/-exclstyle constraints do not exist here; this is an overlap partitioning tool, not a shuffler.- Prefer
-hfor help; the captured--helppath in references is misleading because the wrapper expects other arguments first.