Piping samtools merge to samtools calmd

Question

[rubic]

Hi,

Suppose I have 2 bam files: left.bam and right.bam that I want to merge and then pipe the outputed merged file to samtools calmd with a ref.fasta file. Will the appropriate command be:

samtools -n - left.bam right.bam | samtools calmd -b - ref.fasta

[Carlo Yague]

You forgot the "merge" but other than that, I don't see anything wrong in that cmd. Does it work ?

samtools merge -n - left.bam right.bam | samtools calmd -b - ref.fasta