Question: Piping a bam file to samtools calmd
0
gravatar for rubic
21 months ago by
rubic170
United States
rubic170 wrote:

Is there way to combine these samtools commands:

  1. samtools view -h -F 0x10 <bam>

  2. samtools calmd <bam_from_step_1> <ref.fasta>

I tried this:

samtools view -h -F 0x10 <bam> | samtools calmd <ref.fasta> | samtools view -bS - > <out.bam>

But it doesn't work.

pipe calmd samtools • 977 views
ADD COMMENTlink modified 21 months ago by Dan D6.6k • written 21 months ago by rubic170
1

For one thing, samtools view outputs data in .sam format, not .bam format. So if calmd requires a .bam format, add -b to samtools view top give it a bam.

ADD REPLYlink written 21 months ago by swbarnes24.0k
3
gravatar for Dan D
21 months ago by
Dan D6.6k
Tennessee
Dan D6.6k wrote:

In bash you need to use - to specify that the input is coming from a pipe. You should also add the -S parameter to specify that the input is SAM (or follow the suggestion of @swbarnes2 and output BAM from view instead):

samtools view -h -F 0x10 <bam> | samtools calmd -S - <ref.fasta> | samtools view -bS - > <out.bam>
ADD COMMENTlink modified 21 months ago • written 21 months ago by Dan D6.6k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1198 users visited in the last hour