Version 2.14.1-3 (2-Nov-2011) Bug Fixes ========= 1. Corrected the help for closestBed. It now correctly reads -io instead of -no. 2. Fixed regression in closestBed injected in version 2.13.4 whereby B features to the right of an A feature were missed. New tool. ============ 1. Added the multiIntersectBed tool for reporting common intervals among multiple **sorted** BED/GFF/VCF files. Version 2.13.4 (26-Oct-2011) Bug Fixes ========= 1. The -sorted option (chromsweep) in intersectBed now obeys -s and -S. I had neglected to implement that. Thanks to Paul Ryvkin for pointing this out. 2. The -split option was mistakenly splitting of D CIGAR ops. 3. The Makefile was not including zlib properly for newer versions of GCC. Thanks to Istvan Albert for pointing this out and providing the solution. Improvements ============ 1. Thanks to Jacob Biesinger for a new option (-D) in closestBed that will report _signed_ distances. Moreover, the new option allows fine control over whether the distances are reported based on the reference genome or based on the strand of the A or B feature. Many thanks to Jacob. 2. Thanks to some nice analysis from Paul Ryvkin, I realized that the -sorted option was using way too much memory in certain cases where there is a chromosome change in a sorted BED file. This has been corrected. Version 2.13.3 (30-Sept-2011) Bug Fixes ============ 1. intersectBed detected, but did not report overlaps when using BAM input and -bed. Other ===== 1. Warning that -sorted trusts, but does not enforce that data is actually sorted. Version 2.13.2 (23-Sept-2011) New algorithm ============= 1. Preliminary release of the chrom_sweep algorithm. New options =========== 1. genomeCoverageBed no longer requires a genome file when working with BAM input. It instead uses the BAM header. 2. tagBam now has a -score option for annotating alignments with the BED "scores" field in annotation files. This overrides the default behavior, which is to use the -labels associated with the annotation files passed in on the command line. Bug fixes ========= 1. Correct a bug that prevented proper BAM support in intersectBed. 2. Improved detection of GFF features with negative coordinates. Version 2.13.1 (6-Sept-2011) New options =========== 1. tagBam now has -s and -S options for only annotating alignments with features on the same and opposite strand, respectively. 2. tagBam now has a -names option for annotating alignments with the "name" field in annotation files. This overrides the default behavior, which is to use the -labels associated with the annotation files passed in on the command line. Currently, this works well with BED files, but given the limited metadata support for GFF files, annotating with -names and GFF files may not work as well as wished, depending on the type of GFF file used. Version 2.13.0 (1-Sept-2011) New tools ========= 1. tagBam. This tool annotates a BAM file with custom tag fields based on overlaps with BED/GFF/VCF files. For example: $ tagBam -i aln.bam -files exons.bed introns.bed cpg.bed utrs.bed \ -tags exonic intonic cpg utr \ > aln.tagged.bam For alignments that have overlaps, you should see new BAM tags like "YB:Z:exonic", "YB:Z:cpg;utr" 2. multiBamCov. The new tool counts sequence coverage for multiple bams at specific loci defined in a BED/GFF/VCF file. For example: $ multiBamCov -bams aln.1.bam aln.2.bam aln3.bam -bed exons.bed chr1 861306 861409 SAMD11 1 + 181 280 236 chr1 865533 865718 SAMD11 2 + 249 365 374 chr1 866393 866496 SAMD11 3 + 162 298 322 where the last 3 columns represent the number of alignments overlapping each interval from the three BAM file. The following options are available to control which types of alignments are are counted. -q Minimum mapping quality allowed. Default is 0. -D Include duplicate-marked reads. Default is to count non-duplicates only -F Include failed-QC reads. Default is to count pass-QC reads only -p Only count proper pairs. Default is to count all alignments with MAPQ greater than the -q argument, regardless of the BAM FLAG field. 3. nucBed. This new tool profiles the nucleotide content of intervals in a fasta file. The following information will be reported after each original BED/GFF/VCF entry: 1) %AT content 2) %GC content 3) Number of As observed 4) Number of Cs observed 5) Number of Gs observed 6) Number of Ts observed 7) Number of Ns observed 8) Number of other bases observed 9) The length of the explored sequence/interval. 10) The sequence extracted from the FASTA file. (optional, if -seq is used) 11) The number of times a user defined pattern was observed. (optional, if -pattern is used.) For example: $ nucBed -fi ~/data/genomes/hg18/hg18.fa -bed simrep.bed | head -3 #1_usercol 2_usercol 3_usercol 4_usercol 5_usercol 6_usercol 7_pct_at 8_pct_gc 9_num_A 10_num_C 11_num_G 12_num_T 13_num_N 14_num_oth 15_seq_len chr1 10000 10468 trf 789 + 0.540598 0.459402 155 96 119 98 0 0 468 chr1 10627 10800 trf 346 + 0.445087 0.554913 54 55 41 23 0 0 173 One can also report the sequence itself: $ nucBed -fi ~/data/genomes/hg18/hg18.fa -bed simrep.bed -seq | head -3 #1_usercol 2_usercol 3_usercol 4_usercol 5_usercol 6_usercol 7_pct_at 8_pct_gc 9_num_A 10_num_C 11_num_G 12_num_T 13_num_N 14_num_oth 15_seq_len 16_seq chr1 10000 10468 trf 789 + 0.540598 0.459402 155 96 119 98 0 0 468 ccagggg... chr1 10627 10800 trf 346 + 0.445087 0.554913 54 55 41 23 0 0 173 TCTTTCA... Or, one can count the number of times that a specific pattern occur in the intervals (reported as the last column): $ nucBed -fi ~/data/genomes/hg18/hg18.fa -bed simrep.bed -pattern CGTT | head #1_usercol 2_usercol 3_usercol 4_usercol 5_usercol 6_usercol 7_pct_at 8_pct_gc 9_num_A 10_num_C 11_num_G 12_num_T 13_num_N 14_num_oth 15_seq_len 16_user_patt_count chr1 10000 10468 trf 789 + 0.540598 0.459402 155 96 119 98 0 0 468 0 chr1 10627 10800 trf 346 + 0.445087 0.554913 54 55 41 23 0 0 173 0 chr1 10757 10997 trf 434 + 0.370833 0.629167 49 70 81 40 0 0 240 0 chr1 11225 11447 trf 273 + 0.463964 0.536036 44 86 33 59 0 0 222 0 chr1 11271 11448 trf 187 + 0.463277 0.536723 37 69 26 45 0 0 177 0 chr1 11283 11448 trf 199 + 0.466667 0.533333 37 64 24 40 0 0 165 0 chr1 19305 19443 trf 242 + 0.282609 0.717391 17 57 42 22 0 0 138 1 chr1 20828 20863 trf 70 + 0.428571 0.571429 10 7 13 5 0 0 35 0 chr1 30862 30959 trf 79 + 0.556701 0.443299 35 22 21 19 0 0 97 0 New options =========== 1. Support for named pipes and FIFOs. 2. "-" is now allowable to indicate that data is being sent via stdin. 3. Multiple tools. Added new -S option to annotateBed, closestBed, coverageBed, intersectBed, pairToBed, subtractBed, and windowBed (-Sm). This new option does the opposite of the -s option: that is, overlaps are only processed if they are on _opposite_ strands. Thanks to Sol Katzman for the great suggestion. Very useful for certain RNA-seq analyses. 4. coverageBed. Added a new -counts option to coverageBed that only reports the count of overlaps, instead of also computing fractions, etc. This is much faster and uses much less memory. 5. fastaFromBed. Added a new -full option that uses the full BED entry when naming each output sequence. Also removed the -fo option such that all output is now written to stdout. 6. genomeCoverageBed. - Added new -scale option that allows the coverage values to be scaled by a constant. Useful for normalizing coverage with RPM, RPKM, etc. Thanks to Ryan Dale for the useful suggestion. - Added new -5, -3, -trackline, -trackopts, and -dz options. Many thanks to Assaf Gordon for these improvements. -5: Calculate coverage of 5" positions (instead of entire interval) -3: Calculate coverage of 3" positions (instead of entire interval). -trackline: Adds a UCSC/Genome-Browser track line definition in the first line of the output. -trackopts: rites additional track line definition parameters in the first line. -dz: Report the depth at each genome position with zero-based coordinates, instead of zero-based. 7. closestBed. See below, thanks to Brent Pedersen, Assaf Gordon, Ryan Layer and Dan Webster for the helpful discussions. - closestBed now reports _all_ features in B that overlap A by default. This allows folks to decide which is the "best" overlapping feature on their own. 2. closestBed now has a "-io" option that ignores overlapping features. In other words, it will only report the closest, non-overlapping feature. An example: $ cat a.bed chr1 10 20 $ cat b.bed chr1 15 16 chr1 16 40 chr1 100 1000 chr1 200 1000 $ bin/closestBed -a a.bed -b b.bed chr1 10 20 chr1 15 16 chr1 10 20 chr1 16 40 $ bin/closestBed -a a.bed -b b.bed -io chr1 10 20 chr1 100 1000 Updates ======= 1. Updated to the latest version of BamTools. This allows greater functionality and will facilitate new options and tools in the future. - Bug Fixes ========= 1. GFF files cannot have zero-length features. 2. Corrected an erroneous check on the start coordinates in VCF files. Thanks to Jan Vogel for the correction. 3. mergeBed now always reports output in BED format. 3. Updated the text file Tokenizer function to yield 15% speed improvement. 4. Various tweaks and improvements. Version 2.12.0 (April-3-2011) New Tool ======== 1. Added new tool called "flankBed", which allows one to extract solely the flanking regions that are upstream and downstream of a given feature. Unlike slopBed, flankBed does not include the original feature itself. A new feature is created for each flabking region. For example, imagine the following feature: chr1 100 200 The following would create features for solely the 10 bp regions flanking this feature. $ bin/flankBed -i a.bed -b 10 -g genomes/human.hg18.genome chr1 90 100 chr1 200 210 In contrast, slopBed would return: bin/slopBed -i a.bed -b 10 -g genomes/human.hg18.genome chr1 90 210 FlankBed has all of the same features as slopBed. New Features ============ 1. Added new "-scores" feature to mergeBed. This allows one to take the sum, min, max, mean, median, mode, or antimode of merged feature scores. In addition, one can use the "collapse" operation to get a comma-separated list of the merged scores. 2. mergeBed now tolerates multiple features in a merged block to have the same feature name. 3. Thanks to Erik Garrison's "fastahack" library, fastaFromBed now reports its output in the order of the input file. 4. Added a "-n" option to bed12ToBed6, which forces the score field to be the 1-based block number from the original BED12 feature. This is useful for tracking exon numbers, for example. 5. Thanks to Can Alkan, added a new "-mc" option to maskFastaFromBed that allows one to define a custom mask character, such as "X" (-n X). Bug Fixes ========= 1. Thanks to Davide Cittaro, intersectBed and windowBed now properly capture unmapped BAM alignments when using the "-v" option. 2. ClosestBed now properly handles cases where b.end == a.start 3. Thanks to John Marshall, the default constructors are much safer and less buggy. 4. Fixed bug in shuffleBed that complained about a lack of -incl and -excl. 5. Fixed bug in shuffleBed for features that would go beyond the end of a chromosome. 6. Tweaked bedToIgv to make it more Windows friendly. Version 2.11.2 (January-31-2010) Fixed a coordinate reporting bug in coverageBed. Added "max distance (-d)" argument back to the new implementation of mergeBed. Version 2.11.0 (January-21-2010) Enhancements: ============= 1. Support for zero length features (i.e., start = end) - For example, this allows overlaps to be detected with insertions in the reference genome, as reported by dbSNP. 2. Both 8 and 9 column GFF files are now supported. 3. slopBed can now extend the size of features by a percentage of it's size (-pct) instead of just a fixed number of bases. 4. Two improvements to shuffleBed: 3a. A -f (overlapFraction) parameter that defines the maximum overlap that a randomized feature can have with an -excl feature. That is, if a chosen locus has more than -f overlap with an -excl feature, a new locus is sought. 3b. A new -incl option (thanks to Michael Hoffman and Davide Cittaro) that, defines intervals in which the randomized features should be placed. This is used instead of placing the features randomly in the genome. Note that a genome file is still required so that a randomized feature does not go beyond the end of a chromosome. 5. bamToBed can now optionally report the CIGAR string as an additional field. 6. pairToPair can now report the entire paired feature from the B file when overlaps are found. 7. complementBed now reports all chromosomes, not just those with features in the BED file. 8. Improved randomization seeding in shuffleBed. This prevents identical output for runs of shuffleBed that occur in the same second (often the case). Bug Fixes: ========== 1. Fixed the "BamAlignmentSupportData is private" compilation issue. 2. Fixed a bug in windowBed that caused positions to run off the end of a chromosome. Major Changes: ============== 1. The groupBy command is now part of the filo package (https://github.com/arq5x/filo) and will no longer be distributed with BEDTools. Version 2.10.0 (September-21-2010) ==New tools== 1. annotateBed. Annotates one BED/VCF/GFF file with the coverage and number of overlaps observed from multiple other BED/VCF/GFF files. In this way, it allows one to ask to what degree one feature coincides with multiple other feature types with a single command. For example, the following will annotate the fraction of the variants in variants.bed that are covered by genes, conservaed regions and know variation, respectively. $ annotateBed -i variants.bed -files genes.bed conserv.bed known_var.bed This tool was suggested by Can Alkan and was motivated by the example source code that he kindly provided. ==New features== 1. New frequency operations (freqasc and freqdesc) added to groupBy. These operations report a histogram of the frequency that each value is observed in a given column. 2. Support for writing uncompressed bam with the -ubam option. 3. Shorthand arguments for groupBy (-g eq. -grp, -c eq. -opCols, -o eq. -opCols). 4. In addition, all BEDTools that require only one main input file (the -i file) will assume that input is coming from standard input if the -i parameter is ignored. For example, the following are equivalent: $ cat snps.bed | sortBed –i stdin $ cat snps.bed | sortBed As are these: $ cat data.txt | groupBy -i stdin -g 1,2,3 -c 5 -o mean $ cat data.txt | groupBy -g 1,2,3 -c 5 -o mean ==Bug fixes== 1. Increased the precision of the output from groupBy. Version 2.9.0 (August-16-2010) ==New tools== 1. unionBedGraphs. This is a very powerful new tool contributed by Assaf Gordon from CSHL. It will combine/merge multiple BEDGRAPH files into a single file, thus allowing comparisons of coverage (or any text-value) across multiple samples. ==New features== 1. New "distance feature" (-d) added to closestBed by Erik Arner. In addition to finding the closest feature to each feature in A, the -d option will report the distance to the closest feature in B. Overlapping features have a distance of 0. 2. New "per base depth feature" (-d) added to coverageBed. This reports the per base coverage (1-based) of each feature in file B based on the coverage of features found in file A. For example, this could report the per-base depth of sequencing reads (-a) across each capture target (-b). ==Bug Fixes== 1. Fixed bug in closestBed preventing closest features from being found for A features with start coordinates < 2048000. Thanks to Erik Arner for pointing this out. 2. Fixed minor reporting annoyances in closestBed. Thanks to Erik Arner. 3. Fixed typo/bug in genomeCoverageBed that reported negative coverage owing to numeric overflow. Thanks to Alexander Dobin for the detailed bug report. 4. Fixed other minor parsing and reporting bugs/annoyances. Version 2.8.3 (July-25-2010) 1. Fixed bug that caused some GFF files to be misinterpreted as VCF. This prevented the detection of overlaps. 2. Added a new "-tag" option in bamToBed that allows one to choose the _numeric_ tag that will be used to populate the score field. For example, one could populate the score field with the alignment score with "-tag AS". 3. Updated the BamTools API. Version 2.8.2 (July-18-2010) 1. Fixed a bug in bedFile.h preventing GFF strands from being read properly. 2. Fixed a bug in intersectBed that occasionally caused spurious overlaps between BAM alignments and BED features. 3. Fixed bug in intersectBed causing -r to not report the same result when files are swapped. 4. Added checks to groupBy to prevent the selection of improper opCols and groups. 5. Fixed various compilation issues, esp. for groupBy, bedToBam, and bedToIgv. 6. Updated the usage statements to reflect bed/gff/vcf support. 7. Added new fileType functions for auto-detecting gzipped or regular files. Thanks to Assaf Gordon. Version 2.8.1 (July-05-2010) 1. Added bedToIgv. Version 2.8.0 (July-04-2010) 1. Proper support for "split" BAM alignments and "blocked" BED (aka BED12) features. By using the "-split" option, intersectBed, coverageBed, genomeCoverageBed, and bamToBed will now correctly compute overlaps/coverage solely for the "split" portions of BAM alignments or the "blocks" of BED12 features such as genes. 2. Added native support for the 1000 Genome Variant Calling Format (VCF) version 4.0. 3. New bed12ToBed6 tool. This tool will convert each block of a BED12 feature into discrete BED6 features. 4. Useful new groupBy tool. This is a very useful new tool that mimics the "groupBy" clause in SQL. Given a file or stream that is sorted by the appropriate "grouping columns", groupBy will compute summary statistics on another column in the file or stream. This will work with output from all BEDTools as well as any other tab-delimited file or stream. Example summary operations include: sum, mean, stdev, min, max, etc. Please see the help for the tools for examples. The functionality in groupBy was motivated by helpful discussions with Erik Arner at Riken. 5. Improvements to genomeCoverageBed. Applied several code improvements provided by Gordon Assaf at CSHL. Most notably, beyond the several efficiency and organizational changes he made, he include a "-strand" option which allows one to specify that coverage should only be computed on either the "+" or the "-" strand. 6. Fixed a bug in closestBed found by Erik Arner (Riken) which incorrectly reported "null" overlaps for features that did not have a closest feature in the B file. 7. Fixed a careless bug in slopBed also found by Erik Arner (Riken) that caused an infinite loop when the "-excl" option was used. 8. Reduced memory consumption by ca. 15% and run time by ca. 10% for most tools. 9. Several code-cleanliness updates such as templated functions and common tyedefs. 10. Tweaked the genome binning approach such that 16kb bins are the most granular. Version 2.7.1 (May-06-2010) Fixed a typo that caused some compilers to fail on closestBed. Version 2.7.0 (May-05-2010) General: 1. "Gzipped" BED and GFF files are now supported as input by all BEDTools. Such files must end in ".gz". 2. Tools that process BAM alignments now uniformly compute an ungapped alignment end position based on the BAM CIGAR string. Specifically, "M", "D" and "N" operations are observed when computing the end position. 3. bamToBed requires the BAM file to be sorted/grouped by read id when creating BEDPE output. This allows the alignments end coordinate for each end of the pair to be properly computed based on its CIGAR string. The same requirement applies to pairToBed. 4. Updated manual. 5. Many silent modifications to the code that improve clarity and sanity-checking and facilitate future additions/modifications. New Tools: 1. bedToBam. This utility will convert BED files to BAM format. Both "blocked" (aka BED12) and "unblocked" (e.g. BED6) formats are acceptable. This allows one to, for example, compress large BED files such as dbSNP into BAM format for efficient visualization. Changes to existing tools: intersectBed 1. Added -wao option to report 0 overlap for features in A that do not intersect any features in B. This is an extension of the -wo option. bamToBed 1. Requires that BAM input be sorted/grouped by read name. pairToBed 1. Requires that BAM input be sorted/grouped by read name. 2. Allows use of minimum mapping quality or total edit distance for score field. windowBed 1. Now supports BAM input. genomeCoverageBed 1. -bga option. Thanks to Gordon Assaf for the suggestion. 2. Eliminated potential seg fault. Acknowledgements: 1. Gordon Assaf: for suggesting the -bga option in genomeCoverageBed and for testing the new bedToBam utility. 2. Ivan Gregoretti: for helping to expedite the inclusion of gzip support. 3. Can Alkan: for suggesting the addition of the -wao option to intersectBed. 4. James Ward: for pointing out that bedToBam did not need to create "dummy" seq and qual entries. Version 2.6.1 (Mar-29-2010) 1. Fixed a careless command line parsing bug in coverageBed. Version 2.6.0 (Mar-23-2010) ***Specific improvements / additions to tools*** 1. intersectBed * Added an option (-wo) that reports the number of overlapping bases for each intersection b/w A and B files. -- Not sure why this wasn't added sooner; it's obvious. 2. coverageBed * native BAM support * can now report a histogram (-hist) of coverage for each feature in B. Useful for exome sequencing projects, for example. -- thanks for the excellent suggestion from Jose Bras * faster 3. genomeCoverageBed * native BAM support * can now report coverage in BEDGRAPH format (-bg) -- thanks for the code and great suggestion from Gordon Assaf, CSHL. 4. bamToBed * support for "blocked" BED (aka BED12) format. This facilitates the creation of BED entries for "split" alignments (e.g. RNAseq or SV) -- thanks to Ann Loraine, UNCC for test data to support this addition. 5. fastaFromBed * added the ability to extract sequences from a FASTA file according to the strand in the BED file. That is, when "-" the extracted sequence is reverse complemented. -- thanks to Thomas Doktor, U. of Southern Denmark for the code and suggestion. 6. ***NEW*** overlap * newly added tool for computing the overlap/distance between features on the same line. -- For example: $ cat test.out chr1 10 20 A chr1 15 25 B chr1 10 20 C chr1 25 35 D $ cat test.out | overlaps -i stdin -cols 2,3,6,7 chr1 10 20 A chr1 15 25 B 5 chr1 10 20 C chr1 25 35 D -5 ***Bug fixes*** 1. Fixed a bug in pairToBed when comparing paired-end BAM alignments to BED annotations and using the "notboth" option. 2. Fixed an idiotic bug in intersectBed that occasionally caused segfaults when blank lines existed in BED files. 3. Fixed a minor bug in mergeBed when using the -nms option. ***General changes*** 1. Added a proper class for genomeFiles. The code is much cleaner and the tools are less sensitive to minor problems with the formatting of genome files. Per Gordon Assaf's wise suggestion, the tools now support "chromInfo" files directly downloaded from UCSC. Thanks Gordon---I disagreed at first, but you were right. 2. Cleaned up some of the code and made the API a bit more streamlined. Will facilitate future tool development, etc. Version 2.5.4 (Mar-3-2010) 1. Fixed an insidious bug that caused malform BAM output from intersectBed and pairToBed. The previous BAM files worked fine with samtools as BAM input, but when piped in as SAM, there was an extra tab that thwarted conversion from SAM back to BAM. Many thanks to Ivan Gregoretti for reporting this bug. I had never used the BAM output in this way and thus never caught the bug! Version 2.5.3 (Feb-19-2010) 1. Fixed bug to "re-allow" track and "browser" lines. 2. Fixed bug in reporting BEDPE overlaps. 3. Fixed bug when using type "notboth" with BAM files in pairToBed. 4. When comparing BAM files to BED/GFF annotations with intersectBed or pairToBed, the __aligned__ sequence is used, rather than the __original__ sequence. 5. Greatly increased the speed of pairToBed when using BAM alignments. 6. Fixed a bug in bamToBed when reporting edit distance from certain aligners. Version 2.5.2 (Feb-2-2010) 1. The start and end coordinates for BED and BEDPE entries created by bamToBed are now based on the __aligned__ sequence, rather than the original sequence. It's obvious, but I missed it originally...sorry. 2. Added an error message to mergeBed preventing one from using "-n" and "-nms" together. 3. Fixed a bug in pairToBed that caused neither -type "notispan" nor "notospan" to behave as described. Version 2.5.1 (Jan-28-2010) 1. Fixed a bug in the new GFF/BED determinator that caused a segfault when start = 0. Version 2.5.0 (Jan-27-2010) 1. Added support for custom BED fields after the 6th column. 2. Fixed a command line parsing bug in pairToBed. 3. Improved sanity checking. Version 2.4.2 (Jan-23-2010) 1. Fixed a minor bug in mergeBed when -nms and -s were used together. 2. Improved the command line parsing to prevent the occasional segfault. Version 2.4.1 (Jan-12-2010) 1. Updated BamTools libraries to remove some compilation issues on some systems/compilers. Version 2.4.0 (Jan-11-2010) 1. Added BAM support to intersectBed and pairToBed 2. New bamToBed feature. 3. Added support for GFF features 4. Added support for "blocked" BED format (BED12) 6. Wrote complete manual and included it in distribution. 7. Fixed several minor bugs. 8. Cleaned up code and improved documentation. Version 2.3.3 (12/17/2009) Rewrote complementBed to use a slower but much simpler approach. This resolves several bugs with the previous logic. Version 2.3.2 (11/25/2009) Fixed a bug in subtractBed that prevent a file from subtracting itself when the following is used: $ subtractBed -a test.bed -b test.bed Version 2.3.1 (11/19/2009) Fixed a typo in closestBed that caused all nearby features to be returned instead of just the closest one. Version 2.3.0 (11/18/2009) 1. Added four new tools: - shuffleBed. Randomly permutes the locations of a BED file among a genome. Useful for testing for significant overlap enrichments. - slopBed. Adds a requested number of base pairs to each end of a BED feature. Constrained by the size of each chromosome. - maskFastaFromBed. Masks a FASTA file based on BED coordinates. Useful making custom genome files from targeted capture experiment, etc. - pairToPair. Returns overlaps between two paired-end BED files. This is great for finding structural variants that are private or shared among samples. 2. Increased the speed of intersectBed by nearly 50%. 3. Improved corrected some of the help messages. 4. Improved sanity checking for BED entries. Version 2.2.4 (10/27/2009) 1. Updated the mergeBed documentation to describe the -names option which allows one to report the names of the features that were merged (separated by semicolons). Version 2.2.3 (10/23/2009) 1. Changed windowBed to optionally define "left" and "right" windows based on strand. For example by default, -l 100 and -r 500 will add 100 bases to the left (lower coordinates) of a feature in A when scanning for hits in B and 500 bases to the right (higher coordinates). However if one chooses the -sw option (windows bases on strandedness), the behavior changes. Assume the above example except that a feature in A is on the negative strand ("-"). In this case, -l 100, -r 500 and -sw will add 100 bases to the right (higher coordinates) and 500 bases to the left (lower coordinates). In addition, there is a separate option (-sm) that can optionally force hits in B to only be tracked if they are on the same strand as A. ***NOTE: This replaces the previous -s option and may affect existing pipelines***. Version 2.2.2 (10/20/2009) 1. Improved the speed of genomeCoverageBed by roughly 100 fold. The memory usage is now less than 2.0 Gb. Version 2.2.1 1. Fixed a very obvious bug in subtractBed that caused improper behavior when a feature in A was overlapped by more than one feature in B. Many thanks to folks in the Hannon lab at CSHL for pointing this out. Version 2.2.0 === Notable changes in this release === 1. coverageBed will optionally only count features in BED file A (e.g. sequencing reads) that overlap with the intervals/windows in BED file B on the same strand. This has been requested several times recently and facilitates CHiP-Seq and RNA-Seq experiments. 2. intersectBed can now require a minimum __reciprocal__ overlap between intervals in BED A and BED B. For example, previously, if one used -f 0.90, it required that a feature in B overlap 90% of the feature in A for the "hit" to be reported. If one adds the -r (reciprocal) option, the hit must also cover 90% of the feature in B. This helps to exclude overlaps between say small features in A and large features in B: A ========== B ********************************************************** -f 0.50 (Reported), whereas -f 0.50 -r (Not reported) 3. The score field has been changed to be a string. While this deviates from the UCSC definition, it allows one to track much more meaningful information about a feature/interval. For example, score could now be: 7.31E-05 (a p-value) 0.334577 (mean enrichment) 2:2.2:40:2 (several values encoded in a string) 4. closestBed now, by default, reports __all__ intervals in B that overlap equally with an interval in A. Previously, it merely reported the first such feature that appeared in B. Here's a cartoon explaining the difference. **Prior behavior** A ============== B.1 ++++++++++++++ B.2 ++++++++++++++ B.3 +++++++++ ----------------------------------------- Result = B.1 ++++++++++++++ **Current behavior** A ============== B.1 ++++++++++++++ B.2 ++++++++++++++ B.3 +++++++++ ----------------------------------------- Result = B.1 ++++++++++++++ B.2 ++++++++++++++ Using the -t option, one can also choose to report either the first or the last entry in B in the event of a tie. 5. Several other minor changes to the algorithms have been made to increase speed a bit. VERSION 2.1.2 1. Fixed yet another bug in the parsing of "track" or "browser" lines. Sigh... 2. Change the "score" column (i.e. column 5) to b stored as a string. While this deviates from the UCSC convention, it allows significantly more information to be packed into the column. VERSION 2.1.1 1. Added limits.h to bedFile.h to fix compilation issues on some systems. 2. Fixed bug in testing for "track" or "browser" lines. VERSION 2.1.0 1. Fixed a bug in peIntersectBed that prevented -a from being correctly handled when passed via stdin. 2. Added new functionality to coverageBed that calculates the density of coverage. 3. Fixed bug in geneomCoverageBed. VERSION 2.0.1 1. Added the ability to retain UCSC browser track/browser headers in BED files. VERSION 2.0 1. Sped up the file parsing. ~10-20% increase in speed. 2. Created reportBed() as a common method in the bedFile class. Cleans up the code quite nicely. 3. Added the ability to compare BED files accounting for strandedness. 4. Paired-end intersect. 5. Fixed bug that prevented overlaps from being reported when the overlap fraction requested is 1.0 VERSION 1.2, 04/27/2009. (1eb06115bdf3c49e75793f764a70c3501bb53f33) 1. Added subtractBed. A. Fixed bug that prevented "split" overlaps from being reported. B. Prevented A from being reported if >=1 feature in B completely spans it. 2. Added linksBed. 3. Added the ability to define separate windows for upstream and downstream to windowBed. VERSION 1.1, 04/23/2009. (b74eb1afddca9b70bfa90ba763d4f2981a56f432) Initial release.