Ajeet1699

taoliu / merge_then_call_consensus.sh

Last active March 3, 2025 07:35

Merge multiple peak files (BED format) then call the consensus regions

	#!/bin/bash
	#
	# This script will find the consensus peak regions from peak files (in
	# BED format) of multiple samples by:
	#
	# 1. Converting the peak file of each sample into non-overlapping 3
	# cols BED file and concatenating them;
	#
	# 2. Sorting the concatenated file and Building a genome coverage
	# track in BedGraph, of which the value (the 3rd col) indicates the