Comparing bed files for shifts and enlargements
1
0
Entering edit mode
2.5 years ago
Deniz • 0

I was wondering how I can analyze the changes in two different bed files such as:

  1. Shifts (if start-end shifted from, for example, 10-20 to 15-25
  2. Enlargements (if start-end shifted from, for example, 10-20 to 5-25) 2.a and in which direction (did it enlarge both directions or only single)
  3. Shrinkage (if start-end shifted from, for example, 10-20 to 12-18)

Any help would be very appreciated. Thank you so much!

bed bedtools chip-seq • 969 views
ADD COMMENT
0
Entering edit mode

id assume you want to visualize them?

ADD REPLY
0
Entering edit mode

Do you want to compare the first entry of the first file with the first entry in the second file, the second with the second, and so on?

ADD REPLY
0
Entering edit mode

I want to compare all the entries to each other (more like a double for loop)

ADD REPLY
0
Entering edit mode
2.5 years ago

If you're aiming to visualize differences, you could throw datasets into a HiGlass viewer and lock two tracks (one "control", the other "treatment"). Then browse your regions of interest.

Otherwise, you could pipe bedmap to awk to measure offsets, percentage or fold differences between elements between tracks. Spit the values out to to a text file and then feed that into R to summarize and visualize statistics.

ADD COMMENT

Login before adding your answer.

Traffic: 1825 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6