Skip to content

Commit e94df2d

Browse files
committed
fix chrM issue - finally
1 parent eda677d commit e94df2d

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

modules/postalign_bam.bds

+2-2
Original file line numberDiff line numberDiff line change
@@ -191,7 +191,7 @@ string[] dedup_bam( string bam, string o_dir, string log_o_dir, string group, in
191191

192192
//# Remove duplicates
193193
//# Index final position sorted BAM
194-
sys samtools view -F 1804 -h -u -b $dupmark_bam | awk '($3 != "chrM")' | samtools view -@ 8 -Sb - > $nodup_bam
194+
sys samtools view -F 1804 -h $dupmark_bam | awk '($3 != "chrM")' | samtools view -@ 8 -Sb - > $nodup_bam
195195

196196
//# Index Final BAM file
197197
sys sambamba index -t $nth_dedup $nodup_bam
@@ -330,7 +330,7 @@ string[] dedup_bam_PE( string bam, string o_dir, string log_o_dir, string group,
330330
//# Index final position sorted BAM
331331
//# Create final name sorted BAM
332332
//# ============================
333-
sys samtools view -h -F 1804 -f 2 -u $dupmark_bam | awk '($3 != "chrM") '| samtools view -@ 8 -Sb - > $nodup_bam
333+
sys samtools view -h -F 1804 -f 2 $dupmark_bam | awk '($3 != "chrM") '| samtools view -@ 8 -Sb - > $nodup_bam
334334

335335
sys sambamba index -t $nth_dedup $nodup_bam
336336

0 commit comments

Comments
 (0)