Difference between revisions of "BEDTools"
Moskalenko (talk | contribs) m (Text replace - "|{{#vardefine:intel|}} <!-- E.g. "11.1" --> |{{#vardefine:mpi|}} <!-- E.g. "openmpi/1.3.4" -->" to "") |
Moskalenko (talk | contribs) m (Text replace - "]] {| <!--Main settings - REQUIRED-->" to "]] {|<!--Main settings - REQUIRED-->") |
||
Line 2: | Line 2: | ||
__NOEDITSECTION__ | __NOEDITSECTION__ | ||
[[Category:Software]][[Category:Bioinformatics]][[Category:Genomics]] | [[Category:Software]][[Category:Bioinformatics]][[Category:Genomics]] | ||
− | + | {|<!--Main settings - REQUIRED--> | |
− | {| | ||
− | <!--Main settings - REQUIRED--> | ||
|{{#vardefine:app|bedtools}} | |{{#vardefine:app|bedtools}} | ||
|{{#vardefine:url|http://code.google.com/p/bedtools/}} | |{{#vardefine:url|http://code.google.com/p/bedtools/}} |
Revision as of 18:10, 10 August 2012
|
Description
The BEDTools utilities allow one to address common genomics tasks such as finding feature overlaps and computing coverage. The utilities are largely based on four widely-used file formats: BED, GFF/GTF, VCF, and SAM/BAM. Using BEDTools, one can develop sophisticated pipelines that answer complicated research questions by "streaming" several BEDTools together. The following are examples of common questions that one can address with BEDTools.
- Intersecting two BED files in search of overlapping features.
- Culling/refining/computing coverage for BAM alignments based on genome features.
- Merging overlapping features.
- Screening for paired-end (PE) overlaps between PE sequences and existing genomic features.
- Calculating the depth and breadth of sequence coverage across defined "windows" in a genome.
- Screening for overlaps between "split" alignments and genomic features.
The fact that all of the BEDTools accept input from “standard input (stdin)” allows one to “stream / pipe” several commands together to facilitate more complicated analyses. Also, the tools allow fine control over how output is reported. Most recently, I have added support for sequence alignments in BAM ([1]) format, as well as for features in VCF and GFF, as well as “blocked” BED format. The tools are quite fast and typically finish in a matter of a few seconds, even for large datasets. Upstream documentation for bedtools.
Required Modules
Serial
- bedtools
System Variables
- HPC_{{#uppercase:bedtools}}_DIR - installation directory
- HPC_BEDTOOLS_BIN - the bin directory
- HPC_BEDTOOLS_DATA - BEDTools data directory
- HPC_BEDTOOLS_GENOMES - BEDTools genomes directory
- HPC_BEDTOOLS_SCRIPTS - BEDTools scripts directory