Filtering BED file for interval length

Wondering about a quick and easy way to filter a 3 column bed file for interval length, lets say to toss all of the intervals above 1000bp in length. Thanks.

bedtools bed bedops

• 3,719 views • link

written 6.9 years ago by rbronste

1 answer

awk should work fine:

awk '{if($3-$2 <= 1000) print}' file.bed > file_filtered.bed

• 0 views • link

written 6.9 years ago by alex.zaccaron

More posts like this

awk for interval midpoints and length
written 6.5 years ago by rbronste

Kind of a basic awk question but just trying to take the following intervals: chr1 10001 10150 chr1 10008 10090 chr1 10025 10207 chr1 10058 …
Combining two BED files based on shared criteria
written 6.7 years ago by rbronste

I have a somewhat tricky BED file convergence that I need to produce, well tricky for me anyway. The first set of intervals looks as …
Limiting interval sizes with awk
written 6.9 years ago by rbronste

Kind of a basic awk question. I have a bed file where I have been doing the following to limit interval sizes to those under …
Pulling out interval adjoining regions
written 7.7 years ago by rbronste

Looking for a good way to take a set of intervals and print out an interval set (bed file) that represents regions just upstream and …
Read depth interval comparison
written 7.7 years ago by rbronste

I was wondering about a good method to compare read depth at specific positions (bed file intervals) between two bam files and get as an …
FIMO output that preserves padj values
written 7.9 years ago by rbronste

Hi I am generating a fasta file from a DESeq2 output in the following way: bedtools getfasta -fi /son/genomes/mm10.fa -bed treatment_pw_padj.001_sorted.bed -fo treatment_pw_padj.001_sorted.fa The above …
Excluding intervals with repeats
written 8.4 years ago by rbronste

Is there a quick and easy method for excluding intervals from a bed file that contain any kind of repetitive elements? Thanks.
Filter HT-Seq output for gene subclasses
written 8.5 years ago by rbronste

I was wondering about a straightforward way to filter RNA-seq datasets specifically those coming out of HT-Seq for highly expressed genes in specific subclasses, for …
Extracting interval from fasta header
written 8.5 years ago by rbronste

Hi, Im wondering about the most straightforward way to extract the interval information contained in a fasta header such as the one below, thanks! Also …
Bedops --difference errors
written 9.3 years ago by rbronste

Hi bedops --difference giving me some errors in the final bed file and Im wondering why, how to get a workaround? chr8 9 69893500 5.00E+07 …