This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How Can I Split Bam Into Chromosome (In A Loop) Using Samtools?

i want to split my bam file into chromosome and then index it.

split bam chromosome
• 21,038 views • link
updated by Dan • Dan • written by Hmm • Hmm

3 answers

Depending on your chromosomes, something like:

for chrom in `seq 1 22` X Y
do

    samtools view -bh $BAM chr${chrom} | samtools sort - chr${chrom}
    samtools index chr${chrom}.bam
done
• 178 views • link
updated by Ram Ram • written by brentp brentp

do we need to sort bam everytime we split it...or are we doing it to be on the safe side? Thanks for your response.

• 0 views • link
written by Hmm • Hmm

if your input bam is sorted, then you don't need to resort.

• 0 views • link
written by brentp brentp

bamtools has a "split" command for exactly this purpose.

• 136 views • link
updated by Ram Ram • written by Aaronquinlan Aaronquinlan

Cool. I didn't know that. It's a bit sparse on docs

• 0 views • link
updated by Ram Ram • written by brentp brentp

You can use this tool to do exactly that.

http://code.google.com/p/chrom-bed/

• 0 views • link
written by Dan • Dan

Log in to answer this question.

More posts like this
  • Read File Format .Bam
    written by hudarembang • hudarembang
    <p>My research is related to the manufacture of snp database. constraints that occurs is how to read the file. <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a>. in this case, I …
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • How To Split Chr11 Bam File In To 2 Halves ?
    written by biorepine biorepine
    <p>Dear Biostars</p> <p>I am aware of <a href='http://samtools.sourceforge.net/'>samtools</a> splitting of <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file based on chromosome. But is there any way I can divide a …
  • Picard Mark Duplicate Reads In Galaxy
    written by Tonyzeng • Tonyzeng
    <p>HI, I want to use picard mark duplicate reads to just mark my read with potential duplication. from the manual of picard, we need to …
  • Samtools View Of A Sam File Entire Chromosome
    written by dfernan dfernan
    <p>Hi, </p> <p>I have (what I think) is a valid <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file.</p> <p>When I do:</p> <pre><code>samtools view -S F1i1F.C57.2.Aligned.out.sam chr19_C57 </code></pre> <p>I get the …
  • Error In Converting Sam To Bam By Samtools
    written by Raghav Raghav
    <p>hello every one,</p> <p>after mapping,i encountered an error during converting my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> out file into <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> out put file.</p> <p>[cdac@nbri <a href='http://samtools.sourceforge.net/'>samtools</a>-0.1.19]$ <a …
  • Soap Result Converted To Bam Not Working In Samtool
    written by Sarwar • Sarwar
    <p>I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …
  • Bam File: Change Chromosome Notation
    written by Lisanne • Lisanne
    <p>Dear All,</p> <p>I am using <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files for chip-seq analysis. The chromosome notation in a usual <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file is like: chr1. In my …
  • 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 …