#tail -n +2 AR_full_test1.txt | sort -k1,1 > motif.sorted
#sort -k4,4 ../MergedChipseq.bed > peaks.sorted
#tail -n +2 AR_full_test1.txt | sort -k1,1 | join -1 1 -2 4 -t $'\t' - <(sort -k4,4 ../MergedChipseq.bed)\
 | awk 'BEGIN{OFS="\t"}{offset=$2+0;seq=$3;strand=$5;score=$6;chr=$7;start=$8+0;end=$9+0;peak_center=start+int((end-start)/2);if(strand=="+")\
{motif_start=int(peak_center+offset);motif_end=motif_start+length(seq)}else{motif_end=int(peak_center+offset);\
motif_start=motif_end-length(seq)};print chr,motif_start,motif_end,$1,score,strand}' > Motif_absolute_coords.bed