diff --git a/chipathlon/conf.py b/chipathlon/conf.py
index 38c2691de5a1f183ae0990ab368588a9f704c0ef..b76e1408f658c30f17ad5b717686055f5ab160cd 100644
--- a/chipathlon/conf.py
+++ b/chipathlon/conf.py
@@ -17,7 +17,7 @@ file_extensions = {
     "sai": ["sai"],
     "sam": ["sam"],
     "bam": ["bam"],
-    "bed": ["bed", "narrowPeak", "broadPeak"],
+    "bed": ["bed", "narrowPeak", "broadPeak", "tagAlign"],
     "bwa_genome": ["amb", "ann", "bwt", "pac", "sa"],
     "bowtie2_genome": ["1.bt2", "2.bt2", "3.bt2", "4.bt2", "rev.1.bt2", "rev.2.bt2"],
     "quality": ["quality"],
diff --git a/chipathlon/jobs/modules/peak_call.yaml b/chipathlon/jobs/modules/peak_call.yaml
index 3f68d0ab321cf4bd418abb0fec18637ff016285e..ea10a30691354a4b9381dc5b653513096b9cff14 100644
--- a/chipathlon/jobs/modules/peak_call.yaml
+++ b/chipathlon/jobs/modules/peak_call.yaml
@@ -1,12 +1,28 @@
 peak_call:
   - spp[tool]:
-    - r_spp_nodups:
+    - mv_bed_tagalign:
         inputs:
           - exp.bed:
               type: file
+        additional_inputs: null
+        outputs:
+          - exp.tagAlign:
+              type: file
+    - mv_bed_tagalign:
+        inputs:
           - control.bed:
               type: file
         additional_inputs: null
+        outputs:
+          - control.tagAlign:
+              type: file
+    - r_spp_nodups:
+        inputs:
+          - exp.tagAlign:
+              type: file
+          - control.tagAlign:
+              type: file
+        additional_inputs: null
         outputs:
           - results.narrowPeak:
               type: file
@@ -14,7 +30,7 @@ peak_call:
               type: file
           - results.ccscore:
               type: file
-    - zcat_awk_sort_peaks:
+    - cat_awk_sort_peaks:
         inputs:
           - results.narrowPeak:
               type: file
diff --git a/chipathlon/jobs/params/cat_awk_sort_peaks.yaml b/chipathlon/jobs/params/cat_awk_sort_peaks.yaml
new file mode 100644
index 0000000000000000000000000000000000000000..d4cafc0c49312f38c6372d8d9d6f1bea8fb8ef88
--- /dev/null
+++ b/chipathlon/jobs/params/cat_awk_sort_peaks.yaml
@@ -0,0 +1,21 @@
+cat_awk_sort_peaks:
+  inputs:
+    - bed:
+        type: file
+  additional_inputs: null
+  outputs:
+    - bed:
+        type: file
+  command: cat_spp
+  arguments:
+    - "$inputs.0":
+        changeable: false
+        required: true
+        has_value: false
+    - "$outputs.0":
+        changeable: false
+        required: true
+        has_value: false
+  walltime: 2000
+  memory: 2000
+  cores: 1
diff --git a/chipathlon/jobs/params/mv_bed_tagalign.yaml b/chipathlon/jobs/params/mv_bed_tagalign.yaml
new file mode 100644
index 0000000000000000000000000000000000000000..5b54086ee08470e84fec3718e049739b06e4c78b
--- /dev/null
+++ b/chipathlon/jobs/params/mv_bed_tagalign.yaml
@@ -0,0 +1,21 @@
+zcat_awk_sort_peaks:
+  inputs:
+    - bed:
+        type: file
+  additional_inputs: null
+  outputs:
+    - bed:
+        type: file
+  command: mv_bed_tagalign
+  arguments:
+    - "$inputs.0":
+        changeable: false
+        required: true
+        has_value: false
+    - "$outputs.0":
+        changeable: false
+        required: true
+        has_value: false
+  walltime: 2000
+  memory: 2000
+  cores: 1
diff --git a/chipathlon/jobs/params/sort_awk_sort_peaks.yaml b/chipathlon/jobs/params/sort_awk_sort_peaks.yaml
index 1271f400432f943c366ed65adbb897c3cd159950..3979e687bfd1d7e9caa343664a2d1676664337a1 100644
--- a/chipathlon/jobs/params/sort_awk_sort_peaks.yaml
+++ b/chipathlon/jobs/params/sort_awk_sort_peaks.yaml
@@ -6,7 +6,7 @@ sort_awk_sort_peaks:
   outputs:
     - bed:
         type: file
-  command: sort
+  command: sort_macs2
   arguments:
     - "$inputs.0":
         changeable: false
diff --git a/chipathlon/jobs/params/zcat_awk_sort_peaks.yaml b/chipathlon/jobs/params/zcat_awk_sort_peaks.yaml
index f226897669e410bc443fb593a8ed30a16d8df409..cebf20eaca6eb8cf2c762ad5fa2bdda66508a1d7 100644
--- a/chipathlon/jobs/params/zcat_awk_sort_peaks.yaml
+++ b/chipathlon/jobs/params/zcat_awk_sort_peaks.yaml
@@ -6,7 +6,7 @@ zcat_awk_sort_peaks:
   outputs:
     - bed:
         type: file
-  command: zcat
+  command: zcat_spp
   arguments:
     - "$inputs.0":
         changeable: false
diff --git a/chipathlon/jobs/scripts/cat_spp.sh b/chipathlon/jobs/scripts/cat_spp.sh
new file mode 100644
index 0000000000000000000000000000000000000000..54abba40b3501a37886b393035ccb46b4ce49855
--- /dev/null
+++ b/chipathlon/jobs/scripts/cat_spp.sh
@@ -0,0 +1,3 @@
+#!/bin/bash
+
+/bin/cat "$1" | awk 'BEGIN{OFS="\t"}{print $1,int($2),int($3),$4,$5,$6,$7,$8,$9,$10;}' > "$2"
diff --git a/chipathlon/jobs/scripts/mv_bed_tagalign.sh b/chipathlon/jobs/scripts/mv_bed_tagalign.sh
new file mode 100644
index 0000000000000000000000000000000000000000..f5def69cd0b6e830edbc76cc89d251811434bb4b
--- /dev/null
+++ b/chipathlon/jobs/scripts/mv_bed_tagalign.sh
@@ -0,0 +1,3 @@
+#!/bin/bash
+
+/bin/mv "$1" "$2"
diff --git a/chipathlon/jobs/scripts/sort_wrapper.sh b/chipathlon/jobs/scripts/sort_macs2.sh
similarity index 100%
rename from chipathlon/jobs/scripts/sort_wrapper.sh
rename to chipathlon/jobs/scripts/sort_macs2.sh
diff --git a/chipathlon/jobs/scripts/zcat_wrapper.sh b/chipathlon/jobs/scripts/zcat_spp.sh
similarity index 100%
rename from chipathlon/jobs/scripts/zcat_wrapper.sh
rename to chipathlon/jobs/scripts/zcat_spp.sh