This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BAM to CRAM and BAM recover with smaller size

Hi, I was trying to decide which algorithm from samtools to use for CRAM compression and I noticed the BAM files recovered from CRAM are smaller (-10 GiB) than the original BAM file. Do you know what information I am losing with this transformation?

$ docker run -v $PWD:$PWD quay.io/biocontainers/samtools:1.19.2--h50ea8bc_1 samtools view -@ 14 -T $PWD/hg38.fa -C --output-fmt-option archive -o $PWD/SAMPLE-P_archive.cram $PWD/SAMPLE-P.bam
$ docker run -v $PWD:$PWD quay.io/biocontainers/samtools:1.19.2--h50ea8bc_1 samtools view -@ 14 -T $PWD/hg38.fa --input-fmt-option archive -o $PWD/SAMPLE-P_unarchive.bam $PWD/SAMPLE-P_archive.cram

-rw-rw-r-- 1 where where  61G Mar 14 21:25 SAMPLE-P.bam
-rw-rw-r-- 1 where where 9.0M Mar 14 21:25 SAMPLE-P.bam.bai
-rw-r--r-- 1 root     root      18G Mar 15 00:16 SAMPLE-P_archive.cram
-rw-r--r-- 1 root     root      51G Mar 15 01:02 SAMPLE-P_unarchive.bam
-rw-r--r-- 1 root     root      18G Mar 15 01:23 SAMPLE-P_small.cram
-rw-r--r-- 1 root     root      51G Mar 15 01:47 SAMPLE-P_unsmall.bam
-rw-r--r-- 1 root     root      19G Mar 15 01:58 SAMPLE-P_normal.cram
-rw-r--r-- 1 root     root      51G Mar 15 02:14 SAMPLE-P_unormal.bam
-rw-r--r-- 1 root     root      21G Mar 15 03:13 SAMPLE-P_fast.cram
-rw-r--r-- 1 root     root      51G Mar 15 03:35 SAMPLE-P_unfast.bam
bam samtools cram
• 1,285 views • link
updated by jkbonfield jkbonfield • written by geocarvalho geocarvalho

Don't depend on file sizes for any decisions. Look inside/compare the reads.

• 0 views • link
written by GenoMax GenoMax

Just as gzip -1 to gzip -9 can give different file sizes, so can two identical BAMs be very different in size. That may or may not be the cause. You'd have to uncompress to test. (Note there's no point in --input-fmt-option archive as the input format is self-describing, but I wonder if it somehow enabled archive mode for BAM output, which would indeed be something like gzip -9.)

Best thing though is to convert both to SAM and compare them. Htslib comes with a compare_sam.pl tool in the test directory to aid such things. It's slow as it's not designed for anything other than testing, but it'd maybe help give some confidence.

Also, if you know you'll be sticking to samtools/htslib/noodles derived tools for decoding CRAM then you could also try -O cram,archive,version=3.1 to get maximum compression. ALthough frankly "archive" is typically too extreme IMO. It's good compression, but "small" is often a better tradeoff. Try both and see.

• 0 views • link
written by jkbonfield jkbonfield

0 answers

No answers yet.

Log in to answer this question.

More posts like this