Update human_peaks_analysis.sh

human_peaks_analysis.sh

@@ -14,7 +14,7 @@ bedtools closest -D ref -id -a Human.q0.01.80%.multi.2.mergedpeaks.bed -b Human_
 
 cat Human_new_nearest_allgenes_upstream.txt Human_new_nearest_allgenes_downstream.txt | cut -f7 | sort -u > human.closest.genes.txt
 
-at Human_new_nearest_allgenes_upstream.txt Human_new_nearest_allgenes_downstream.txt | cut -f7 | sort | uniq -c | awk '{OFS="\t"; print $1,$2}' | sort -k1,1nr > Human_number_of_peaks_close_to_genes.txt
+cat Human_new_nearest_allgenes_upstream.txt Human_new_nearest_allgenes_downstream.txt | cut -f7 | sort | uniq -c | awk '{OFS="\t"; print $1,$2}' | sort -k1,1nr > Human_number_of_peaks_close_to_genes.txt
 
 cut -f 2 Human_number_of_peaks_close_to_genes.txt | grep -f - <(grep -v "^#" ../../gencode.v33.chr_patch_hapl_scaff.annotation.gtf | awk '{OFS="\t"; if ($3=="gene") print}') | awk '{print $14}' | sed 's/"//g' | sed 's/;//g' | sort -u > list_of_human_closest_genes.txt