This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Merge Runlevel Bam Files?

Illumina HiSeq data: My data was run on 2 flowcells on illumina HiSeq data. I have processed the data and at a point where i have flow cell A and B runlevel data merged. Can anyone tell me the command/tool which i can use to merge flowcell A bam to flow cell B bam. My bams are indexed and sorted.

Thanks

illumina hiseq bam merge
• 3,453 views • link
updated by Pablo Pablo • written by Hmm • Hmm

1 answer

You should use "samtools merge"

Usage: samtools merge out.bam in1.bam in2.bam

For details, here is the man page: http://samtools.sourceforge.net/samtools.shtml

• 0 views • link
written by Pablo Pablo

do people use MergeSamFiles.jar from the picard tools merge files?

• 0 views • link
written by Hmm • Hmm

Log in to answer this question.

More posts like this
  • snakemake to merge bams from multiple samples
    written by samuel samuel
    Hi, I have multiple sample bams which are split over two lanes. e.g. for 2 samples: flowcells: ["CCGB1ANXX-", "CCHNTANXX-"] samples: ["3666-02-45-1_S2", "3666-03-45-1_S3"] lanes: ["L005", "L006"] …
  • Does bamtools split preserve sort and rmdup?
    written by ccnn • ccnn
    Taking over someone else's code, I've discovered a bottleneck. This is a portion of a script that accepts a BAM and outputs smaller, sorted, rmdup-ed, …
  • Merging two fastq files
    written by addilynn.beach • addilynn.beach
    I used two lanes on a flow cell for the same library. I have the fastq files now, and want to merge the data together …
  • Check BAMs are complete and clean
    written by John John
    Hello all :) Is there anyway to check if a BAM file is complete and obeys all the BAM standards? I have created 100 processed …
  • Understanding Illumina Hisq Output
    written by progression1990 • progression1990
    I have output from an illumina HiSeq run, a series of html, zip, and fastq.gz files which are confusing to me. 'X' takes the value …
  • Illumina CASAVA runParameters.xml
    written by rasool.aysha • rasool.aysha
    Hello, In the `runParameters.xml` script of CASAVA 1.8, line 5 of the setup `<PairEndFC>false</PairEndFC>` is confusing to me since I am entering the identity of …
  • Fastq Files From Different Flowcells
    written by hellbio hellbio
    <p>Hi,</p> <p>For a single sample, i have several paired-end fastq files from four different flowcells. i.e. fastq files from different lanes from each flowcell. Instead …
  • Somatic Snp Mutation Igv Codon
    written by Hmm • Hmm
    <p>I am looking at a tumor/normal pair illumina HiSeq data paired end. When i open my normal and tumor <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> files in IGV, i …
  • Svdetect-Structure And Cnv Detection
    written by Hmm • Hmm
    <p>I have patient tumor/normal HiSeq paired end data. I have 2 <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> files which were cleaned/sorted/duplicates removed.One <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> for tumor and other for …
  • Is My Bam File Sorted ?
    written by Pierre Lindenbaum Pierre Lindenbaum
    <p>Hi all,</p> <p>I wonder if there is any tool (in <a href='http://samtools.sourceforge.net/'>samtools</a> ?) to quickly know whether a <strong>BAM</strong> file has been <strong>sorted</strong> or not.</p> …