У меня есть два набора данных из MedIP-seq, которые пытаются проанализировать и сравнить их. Мой рабочий процесс использовал bwa для их выравнивания, используя метилQA и вызовы macs peak. Прямо сейчас у меня есть вывод вызова macs peak в формате bed и bdg. Я не знаю, как мне их анализировать и призывать к дифференциально метилированным регионам. Мне также нужно знать, насколько эти два образца похожи в пиках. Не могли бы вы помочь мне разобраться в этом вопросе?