cmd = step3.4_peakcall_OA_Hip_R1_peaks.narrowPeak hg38 OA_Hip_open_peaks_motif