Does GATK SetNmMdAndUqTags reduces the size of a CRAM?
Entering edit mode
9 months ago
Sd • 0

I performed GATK SetNmMdAndUqTags on a CRAM file for Whole Genome Sequencing after completing the MarkDuplicates step. The initial size of the CRAM file was 19GB, and after performing the SetNmMdAndUqTags operation, its size reduced to 8GB. The following is the log file I obtained after completing the SetNmMdAndUqTags step. I am unsure whether it executed successfully or not. Can someone help me with this?

> [Tue Oct 17 22:13:26 EDT 2023] picard.sam.SetNmMdAndUqTags done. Elapsed time: 164.26 minutes.
To get help, see
java.lang.ArrayIndexOutOfBoundsException: Index 90338345 out of bounds for length 90338345
        at htsjdk.samtools.util.SequenceUtil.sumQualitiesOfMismatches(
        at picard.sam.AbstractAlignmentMerger.fixUq(
        at picard.sam.AbstractAlignmentMerger.fixNmMdAndUq(
        at picard.sam.SetNmMdAndUqTags.fixRecord(
        at picard.sam.SetNmMdAndUqTags.lambda$doWork$0(
        at java.base/$15$1.accept(
        at java.base/java.util.Iterator.forEachRemaining(
        at java.base/java.util.Spliterators$IteratorSpliterator.forEachRemaining(
        at java.base/
        at java.base/
        at java.base/$ForEachOp.evaluateSequential(
        at java.base/$ForEachOp$OfRef.evaluateSequential(
        at java.base/
        at java.base/
        at picard.sam.SetNmMdAndUqTags.doWork(
        at picard.cmdline.CommandLineProgram.instanceMain(
        at org.broadinstitute.hellbender.cmdline.PicardCommandLineProgramExecutor.instanceMain(
        at org.broadinstitute.hellbender.Main.runCommandLineProgram(
        at org.broadinstitute.hellbender.Main.mainEntry(
        at org.broadinstitute.hellbender.Main.main(
CRAM WGS sam GATK SetNmMdAndUqTags • 350 views
Entering edit mode

No theres's an error.

this error happens when the END position of the read is greater than the chromosome.

how did you map the reads ? check the BAM with


Login before adding your answer.

Traffic: 3569 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6