File: fix_sort_samtools13

package info (click to toggle)
reapr 1.0.18%2Bdfsg-5
  • links: PTS, VCS
  • area: main
  • in suites: bookworm, bullseye
  • size: 1,088 kB
  • sloc: cpp: 4,816; perl: 1,467; sh: 160; makefile: 122
file content (16 lines) | stat: -rw-r--r-- 671 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
Description: adapt samtools sort usage to 1.3 syntax 
 Samtools 1.3 changed the syntax of the 'samtools sort' command, breaking
 REAPR's smaltmap command in the process.
 Discussed and forwarded to upstream in person, hence no URL.
Author: Sascha Steinbiss <satta@debian.org>
--- a/src/task_smaltmap.pl
+++ b/src/task_smaltmap.pl
@@ -155,7 +155,7 @@
                . " | $samtools view -S -T $assembly -b - > $raw_bam";
 
 # sort the bam by coordinate
-push @commands, "$samtools sort $raw_bam $raw_bam.sort";
+push @commands, "$samtools sort $raw_bam -O bam -o $raw_bam.sort.bam";
 
 # remove duplicates
 push @commands, "$samtools rmdup $raw_bam.sort.bam $rmdup_bam";