Scroll to navigation

BEDMAP(1) User Commands BEDMAP(1)

NAME

bedmap - manual page for bedmap 2.4.41+dfsg

DESCRIPTION

bedmap

citation: http://bioinformatics.oxfordjournals.org/content/28/14/1919.abstract
https://doi.org/10.1093/bioinformatics/bts277
2.4.41 (typical)
Shane Neph & Scott Kuehn
USAGE: bedmap [process-flags] [overlap-option] <operation(s)...> <ref-file> [map-file]
Any input file must be sorted per the sort-bed utility. The program accepts BED and Starch file formats. You may use '-' for a BED file to indicate the input comes from stdin.
Traverse <ref-file>, while applying <operation(s)> on qualified, overlapping elements from
<map-file>.
Output is one line for each line in <ref-file>, sent to standard output. There
is no limit on the number of operations you can specify to compute in one bedmap call.
If <map-file> is omitted, the given file is treated as both the <ref-file> and <map-file>.
This usage is more efficient than specifying the same file twice.
Arguments may be given in any order before the input file(s).
Process Flags:

--------

Jump to and process data for given <chromosome> only.
Change output delimiter from '|' to <delim> between columns (e.g. '\t').
Error check all input files (slower).
(advanced) Strong input assumptions are made. Compatible with: --bp-ovr, --range, --fraction-both, and --exact overlap options only.
Accept headers (VCF, GFF, SAM, BED, WIG) in any input file.
Print this message and exit successfully.
Minimize memory usage (slower).
Change delimiter of multi-value output columns from ';' to <delim>.
Change the post-decimal precision of scores to <int>. 0 <= <int>.
Use scientific notation for score outputs.
Print no output for a row with no mapped elements.
Ensure <map-file> is read completely (helps to prevent broken pipes).
Print <val> on unmapped --echo-map* and --min/max-element* operations. The default is to print nothing.
Print program information.
By default, --bp-ovr 1 is used):

--------

Require <int> bp overlap between elements of input files.
First 3 fields from <map-file> must be identical to <ref-file>'s.
Both --fraction-ref <val> and --fraction-map <val> must be true to qualify as overlapping. Expect 0 < val <= 1.
Either --fraction-ref <val> or --fraction-map <val> must be true to qualify as overlapping. Expect 0 < val <= 1.
The fraction of the element's size from <map-file> that must overlap the element in <ref-file>. Expect 0 < val <= 1.
The fraction of the element's size from <ref-file> that must overlap an element in <map-file>. Expect 0 < val <= 1.
Grab <map-file> elements within <int> bp of <ref-file>'s element, where 0 <= int. --range 0 is an alias for --bp-ovr 1.
(Any number of operations may be used any number of times.)

----------

SCORE:
<ref-file> must have at least 3 columns and <map-file> 5 columns.
The result of --stdev divided by the result of --mean.
Generalized median. Report the value, x, such that the fraction <val> of overlapping elements' scores from <map-file> is less than x, and the fraction 1-<val> of scores is greater than x. 0 < val <= 1.
The median absolute deviation of overlapping elements in <map-file>. Multiply mad score by <mult>. 0 < mult, and mult is 1 by default.
The highest score from overlapping elements in <map-file>.
A (non-random) highest-scoring and overlapping element in <map-file>.
A random highest-scoring and overlapping element in <map-file>.
The average score from overlapping elements in <map-file>.
The median score from overlapping elements in <map-file>.
The lowest score from overlapping elements in <map-file>.
A (non-random) lowest-scoring and overlapping element in <map-file>.
A random lowest-scoring and overlapping element in <map-file>.
The square root of the result of --variance.
Accumulated scores from overlapping elements in <map-file>.
The mean score from overlapping elements in <map-file>, after ignoring the bottom <low> and top <hi> fractions of those scores. 0 <= low <= 1. 0 <= hi <= 1. low+hi <= 1.
The variance of scores from overlapping elements in <map-file>.
Weighted mean, scaled in proportion to the coverage of the <ref-file> element by each overlapping <map-file> element.

----------

NON-SCORE:
<ref-file> must have at least 3 columns. For --echo-map-id/echo-map-id-uniq, <map-file> must have at least 4 columns. For --echo-map-score, <map-file> must have at least 5 columns. For all others, <map-file> requires at least 3 columns.
The total number of overlapping bases from <map-file>.
The number of distinct bases from <ref-file>'s element covered by overlapping elements in <map-file>.
The fraction of distinct bases from <ref-file>'s element covered by overlapping elements in <map-file>.
The number of overlapping elements in <map-file>.
Print each line from <ref-file>.
List all overlapping elements from <map-file>.
List IDs from all overlapping <map-file> elements.
List unique IDs from overlapping <map-file> elements.
Print genomic range of overlapping elements from <map-file>.
List scores from overlapping <map-file> elements.
List the full length of every overlapping element.

--echo-overlap-size List lengths of overlaps.

Print the first 3 fields of <ref-file> using chrom:start-end format.
Print 'id-' followed by the line number of <ref-file>.
Print the length of each line from <ref-file>.
Print 1 if there exists an overlapping element in <map-file>, 0 otherwise.

SEE ALSO

The full documentation for bedmap is maintained as a Texinfo manual. If the info and bedmap programs are properly installed at your site, the command

info bedmap

should give you access to the complete manual.

March 2024 bedmap 2.4.41+dfsg