BEDOPS

BEDOPS is a toolkit that performs highly efficient and scalable Boolean and other set operations, statistical calculations, archiving, conversion and other management of genomic data of arbitrary scale.

BEDOPS is available as a module on Apocrita.

Usage

To run the default installed version of BEDOPS, simply load the bedops module:

module load bedops

For usage documentation, pass the --help switch to any BEDOPS program.

Example job

Serial job

Here is an example job running on 1 core and 1GB of memory:

#!/bin/bash
#SBATCH --ntasks=1        # (or -n 1) Request 1 core
#SBATCH --mem-per-cpu=1G  # Request 1GB RAM per core (1GB total)
#SBATCH --time=1:0:0      # (or -t 1:0:0) Request 1 hour runtime

module load bedops

# Convert a vcf file to bed format
vcf2bed < test.vcf > test.bed

# Sort the input data
sort-bed test.bed > test_sorted.bed

# Count number of single-base or greated overlaps
bedmap --ec \
       --delim "\t" \
       --bp-ovr 1 \
       --echo \
       --count \
       test_sorted.bed > test_matrix.bed

References

• BEDOPS GitHub
• BEDOPS website