cut -f1,2,3,4 MergedChipseq_summit1.overlap_rmsk > MergedChipseq1.bed
bedtools getfasta -fi /home/Resource/Genome/hg38/hg38.fa -bed MergedChipseq1.bed -split -name -fo MergedChipseq_peaks1.fa
######## Get the meme file for the AR motif
cp /home/zhangbo/Tools/motif_databases/JASPAR_2023/MA0007.3.meme .
#######Run fimo on hpc1
fimo --oc OUTPUT_fimo2 MA0007.3.meme MergedChipseq_peaks1.fa
