Bedtools

Opis programu:
Collectively, the bedtools utilities are a swiss-army knife of tools for a wide-range of genomics analysis tasks. The most widely-used tools enable genome arithmetic: that is, set theory on the genome. For example, bedtools allows one to intersect, merge, count, complement, and shuffle genomic intervals from multiple files in widely-used genomic file formats such as BAM, BED, GFF/GTF, VCF. While each individual tool is designed to do a relatively simple task (e.g., intersect two interval files), quite sophisticated analyses can be conducted by combining multiple bedtools operations on the UNIX command line.

Wersje:

Name:           bedtools
Version:        2.21.0
Platform:       x86_64
Category:       Applications/Engineering
URL     http://bedtools.readthedocs.org/en/latest/index.html
Provides:       bedtools

Poniżej znajduje się informacja jak załadować moduł do oprogramowania:
Odnośnik do ładowania modułów

Przykładowy skrypt uruchomieniowy SLURM:

#!/bin/bash
#SBATCH --nodes=1
#SBATCH --ntasks-per-node=4
#SBATCH --mem=4gb
#SBATCH --time=01:00:00

# Ustawiamy sciezki lub ladujemy odpowiednie moduly
module load bedtools
module load gcc

# Ustawiamy zmienna $TMPDIR
export TMPDIR=$HOME/$SLURM_JOB_ACCOUNT/scratch/$USER/$SLURM_JOB_ID

# Ustawiamy zmienne aplikacji
export SCR=${TMPDIR}

# Ustawiamy zmienne pomocnicze
INPUT_DIR="input"
OUTPUT_DIR="output"
OUTPUT_FILE="OUTPUT"

# Tworzymy katalog tymczasowy
mkdir -p ${TMPDIR}

# Kopiujemy dane wejsciowe do katalogu wskazywanego zmienna $TMPDIR
cp ${SLURM_SUBMIT_DIR}/${INPUT_DIR}/* ${TMPDIR}

# Przechodzimy do katalogu $TMPDIR
cd $TMPDIR

# Wykonujemy obliczenia
bedtools intersect -a cpg.bed -b exons.bed | head -5 > wynik.txt

# Stopka
cat << EOF
-------------------------------------------------------------------------------

End of calculations [$(date)].

-------------------------------------------------------------------------------
EOF

# Konczymy obliczenia, zawartosc katalogu $TMPDIR/output kopiujemy 
# do katalogu z ktorego zakolejkowano zadanie.
mkdir $SLURM_SUBMIT_DIR/${OUTPUT_DIR}
cp -r $TMPDIR/* $SLURM_SUBMIT_DIR/${OUTPUT_DIR}/

# Czyscimy katalog roboczy
rm -rf $TMPDIR

Poniżej znajdują się pliki które powinny znaleźć się w katalogu input.
Plik:Bedtools input.zip