I am getting BAM files with Bismark mapping with bowtie2. This gives you an unordered BAM file Thus, I am sorting the BAM file with the command
samtools sort -O bam -T tmp_ -o Name_sorted.bam Unordered.bam
The unordered bam file can have a 2,6Gb size, but the ordered one drops down to 1,2Gb. The same happens with all the bam files I am ordering. They are reduced to a half.
Is that normal? Is that caused by a different compression level ?
In any case, makes me feel sort of uncomfortable..
2 answers
Sorting helps to give a better compression ratio. There's an old discussion about it here, or a tip of the day here
edit: I should also mention that you'll see a much better compression ratio with RNA seq than with DNA, as RNA seq should have a lot of repeated sequences.
You are right. They have the same number of reads
Log in to answer this question.
-
samtools sort
written by Javier •
-
Sorting an already sorted file with samtools reduced its size
written by shhsvri •
-
samtools sort: legacy usage
written by c.e.chong •
-
Need some clarifications about the methylation mapper Bismark
written by Antonio R. Franco
-
bwa samtool giving empty bam files while piping
written by Goku •
-
RRBS: import Bismark output to RnBeads
written by cvu
-
samtools not sorting bam file
written by sanariyaz21 •
-
samtools sorting and indexing
written by ggman
-
Script to sort bam file not working
written by emth4657 •
-
Can A Bam File Be Sorted By The Read Names
written by ivivek_ngs
Use
samtools flagstatto check the number of reads. If they are the same, you're fine and it is due to compression. How did you create the unsorted bam file?I mentioned that. I am mapping BS-Seq sequences with Bismark
Yes, it is normal. Sorting can save a lot of space.