Ecosyste.ms: Awesome

An open API service indexing awesome lists of open source software.

Awesome Lists | Featured Topics | Projects

https://github.com/arq5x/bedtools2

bedtools - the swiss army knife for genome arithmetic
https://github.com/arq5x/bedtools2

Last synced: 9 days ago
JSON representation

bedtools - the swiss army knife for genome arithmetic

Awesome Lists containing this project

README 

        
![shamalamadingdong](https://github.com/arq5x/bedtools2/blob/master/docs/bedtools.swiss.png "foo")



bedtools - the swiss army knife for genome arithmetic         

=====================================================



[Download current version](https://github.com/arq5x/bedtools2/releases/latest)



[Documentation](http://bedtools.readthedocs.org)



[Cheat-sheet from Ilya Levantis](https://gist.github.com/ilevantis/6d6ecf8718a5803acff736c2dffc933e)



Summary

-------

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.



Performance

-----------

As of version 2.18, ``bedtools`` is substantially more scalable thanks to improvements we have made in the algorithm used to process datasets that are pre-sorted

by chromosome and start position. As you can see in the plots below, the speed and memory consumption scale nicely

with sorted data as compared to the poor scaling for unsorted data. The current version of bedtools intersect is as fast as (or slightly faster) than the ``bedops`` package's ``bedmap`` which uses a similar algorithm for sorted data.  The plots below represent counting the number of intersecting alignments from exome capture BAM files against CCDS exons.

The alignments have been converted to BED to facilitate comparisons to ``bedops``. We compare to the bedmap ``--ec`` option because similar error checking is enforced by ``bedtools``.



**Note:** bedtools could not complete when using 100 million alignments and the R-Tree algorithm used for unsorted data.



![Speed Comparison](http://bedtools.readthedocs.org/en/latest/_images/speed-comparo.png)

![Memory Comparison](http://bedtools.readthedocs.org/en/latest/_images/memory-comparo.png)



Details

-------

First created through urgency and adrenaline by Aaron Quinlan Spring 2009. 

Maintained by the Quinlan Laboratory at the University of Virginia.



1. **Lead developers**:           Aaron Quinlan, Hao Hou, Brent Pedersen, Neil Kindlon

2. **Significant contributions**: John Marshall, Assaf Gordon, Royden Clark, Ryan Dale

3. **Repository**:                https://github.com/arq5x/bedtools2

4. **Stable releases**:           https://github.com/arq5x/bedtools2/releases

5. **Documentation**:             http://bedtools.readthedocs.org

6. **License**:                   Released under MIT license



Citation

--------

*Please cite the following article if you use BEDTools in your research*:

  * Quinlan AR and Hall IM, 2010. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 26, 6, pp. 841–842. 



Also, if you use *pybedtools*, please cite the following.

  * Dale RK, Pedersen BS, and Quinlan AR. Pybedtools: a flexible Python library for manipulating genomic datasets and annotations. Bioinformatics (2011). doi:10.1093/bioinformatics/btr539