Question: bedtools closest default?
0
gravatar for vctrm67
7 days ago by
vctrm6710
vctrm6710 wrote:

I'm confused about the default behavior of bedtools closest. I ran the command bedtools closest -a a.bed -b b.bed -k 1 > closest.bed and got this as a part of the output (the columns correspond to chrom, start in a.bed and chrom, start in b.bed):

9   5819608 9   5819608
9   5819608 9   6080573
9   5819608 9   6164857
9   5819608 9   6278938
9   5819608 9   6506183
9   5819608 9   6507818
9   5819608 9   6832603
9   5819608 9   6957619
9   5819608 9   7010222

However, the only region in a.bed this corresponds to is (the columns correspond to chrom, start, end):

9       5819608 8043573

and the only segments in b.bed that are relevant are (the columns correspond to chrom, start, end):

9   5086998 5326535
9   5326536 5819607
9   5819608 6080572
9   6080573 6164856
9   6164857 6278937
9   6278938 6506182
9   6506183 6507817
9   6507818 6832602
9   6832603 6957618
9   6957619 7010221
9   7010222 8043572

But it was my understanding that, for default behavior, bedtools closest only reports multiple regions in b.bed if there is a tie. However, I don't see a tie here; the only closest regions should be the first and last regions of the segments posted in b.bed, no?

bedtools • 51 views
ADD COMMENTlink modified 7 days ago • written 7 days ago by vctrm6710
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1576 users visited in the last hour