Piping samtools merge to samtools calmd

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

samtools pipe

• 3,748 views • link

written 9.4 years ago by rubic

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

• 0 views • link

written 9.4 years ago by Carlo Yague

0 answers

No answers yet.

More posts like this

will using samtools merge to multiple bam files also remove the redundancy ?
written 2.8 years ago by naveedhasan2000 •

So i have merged some bam files using the samtools merge -n command, i wanna know if the merged bam file which is made now …
Extract primary STAR mapped reads to specific chromosomes
written 4.5 years ago by rubic

Hi, I have RNA-seq reads which I mapped with STAR and I'd like to extract from the bam file the primary alignments that map to …
Merge bam files and get VCF file
written 6.9 years ago by khatami.mahshid •

I have 2 sorted Bam files and I want to merge them and get a VCF file from them using samtools pileup. first i merge …
Are merged bam files and merge fastq file -> bam files same?
written 8.4 years ago by woongjaej •

Hi, guys I'm processing NGS data and have a question. I need to make my data have over 100,000,000 reads, so when my first processing …
samtools sorting and indexing
written 9.1 years ago by ggman

Hi friends, I am attempting to sort my bam files that I obtained from my bowtie sam files. I am not indexing them appropriate according …
Piping a bam file to samtools calmd
written 9.6 years ago by rubic

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 …
Merge BAM files without retaining sample information
written 10.2 years ago by yampolsk •

I want to merge several groups of BAM files in such a way that when I then create an mpileup from the resulting merged files, …
Merge and fixmate of multiple BAM files generated from pair-end DNA sequences
written 10.9 years ago by lisadavic66 •

Dear All, I have multiple name-sorted BAM files generated from pair-end DNA sequences which look like: `R1-1.bam`, `R1-2.bam`, ..., `R2-1.bam`, `R2-2.bam`, ..., and I am …
Size of BAM file reduces after sorting with samtools
written 11.1 years ago by shakeelbiochemist •

I have 3 BAM files of the same specie, each of ~7GB, from three experimental runs. I merged the three BAM files to produce a …
How To Get Ref.Fasta
written 15.1 years ago by Zhshqzyc •

Hi, I want to use samtools command <pre><code>samtools faidx <ref.fasta> [region1 [...]] </code></pre> My question: where can I get ref.fasta or how to create ref.fasta …