Hi I believe the merge function may not be functioning properly.
I'm executing the following command:
bedtools merge -s -c 6 -o distinct -d 200 -i Athaliana_hmm_e00005.bed
My input data:
1 11896548 11896704 . . -
1 11896785 11897244 . . -
1 12678495 12679080 . . +
1 12723287 12723953 . . -
1 12833400 12834042 . . -
1 13187287 13188100 . . -
1 13269707 13270250 . . -
1 13611137 13611404 . . -
1 14049097 14049739 . . +
1 14226222 14226468 . . -
1 14226532 14226739 . . -
And I get the following output:
1 11896785 11897244 -
1 12678495 12679080 +
1 12723287 12723953 -
1 12833400 12834042 -
1 13187287 13188100 -
1 13269707 13270250 -
1 13611137 13611404 -
1 14049097 14049739 +
1 14226532 14226739 -
I have signaled in bold the problematic lines. As it can be seen in the output, instead of merging the lines, and keeping the lowest start position and the highest end of close ranges, the program only deletes a line and doesn't change the values.
I hope you can help me see if I'm making a mistake when using the program or if this is actually a bug and can be fixed. Thank you!
Navei.