cmd = MergedChipseq_peaks.fa hg38 output