Size of BAM file reduces after sorting with samtools
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 single 22 GB bam file, using samtools merge -r option. Then I sorted this merged bam file with samtools sort, and i got 11 GB merged bam. Is is possible to reduce the size of merged bam file by 50%??
• 12,207 views
•
link
updated
by
Sean Davis
•
written
by
shakeelbiochemist •
•
written
by
shakeelbiochemist •
1 answer
When you sort by coordinate, you bring reads with similar sequences next to each other, allowing the compression algorithm to see more compressible content. It is worth checking, though, that the number of sequences is what you expect using samtools flagstat or simply samtools view and a wc -l.
• 0 views
•
link
updated
by
Ram
•
written
by
Sean Davis
•
written
by
Sean Davis
Log in to answer this question.
More posts like this
-
Merging bam files from different runs
written by shpak.max •
-
samtools index error - "Unsorted positions on sequence" even though file is sorted
written by Geoffrey •
-
Merge bam files and get VCF file
written by khatami.mahshid •
-
Sorting an already sorted file with samtools reduced its size
written by shhsvri •
-
samtools sort doesn't finish sorting bam file, just hangs for hours
written by dec986
-
How group merged VCFs in one GT column
written by miqrom •
-
Merge and fixmate of multiple BAM files generated from pair-end DNA sequences
written by lisadavic66 •
-
Problem With Sorting Bam File: Samtools
written by lyfsa •
-
Samtools Sort
written by GPR
-
Is My Bam File Sorted ?
written by Pierre Lindenbaum
Yes
You can use
samtools flagstat .bamto check read counts etc. for the different files.