Question: Picard MarkDuplicates /SAM validation error
1
gravatar for whaiyu06
9 months ago by
whaiyu0620
whaiyu0620 wrote:

Hello,everyone

I used Picardtools/MarkDuplicates as following:

java -jar picard.jar \
    INPUT=A01-2m-D6.sort.bam \
    OUTPUT=A01-2m-D6.bam \
        METRICS_FILE=A01-2m-D6.metrics CREATE_INDEX=true 2>A01-2m-D6.picard_MarkDuplicates.stderr

then I met thus error:

picard.sam.markduplicates.MarkDuplicates done. Elapsed time: 2.38 minutes.
Runtime.totalMemory()=2058354688
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread "main" htsjdk.samtools.SAMFormatException: SAM validation error: ERROR: Record 78004, Read name E00477:238:hggvjccxy:8:2224:1519:63050, Zero-length read without FZ, CS or CQ tag
        at htsjdk.samtools.SAMUtils.processValidationErrors(SAMUtils.java:454)
        at htsjdk.samtools.BAMFileReader$BAMFileIterator.advance(BAMFileReader.java:812)
        at htsjdk.samtools.BAMFileReader$BAMFileIterator.next(BAMFileReader.java:797)
        at htsjdk.samtools.BAMFileReader$BAMFileIterator.next(BAMFileReader.java:765)
        at htsjdk.samtools.SamReader$AssertingIterator.next(SamReader.java:576)
        at htsjdk.samtools.SamReader$AssertingIterator.next(SamReader.java:548)
        at picard.sam.markduplicates.MarkDuplicates.buildSortedReadEndLists(MarkDuplicates.java:495)
        at picard.sam.markduplicates.MarkDuplicates.doWork(MarkDuplicates.java:232)
        at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:269)
        at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:98)
        at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:108)

How can I deal with this problem ?Thanks

sequencing alignment • 750 views
ADD COMMENTlink modified 9 months ago by Pierre Lindenbaum115k • written 9 months ago by whaiyu0620
6
gravatar for Pierre Lindenbaum
9 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum115k wrote:

use https://broadinstitute.github.io/picard/command-line-overview.html#StandardOptions

VALIDATION_STRINGENCY=LENIENT
ADD COMMENTlink written 9 months ago by Pierre Lindenbaum115k

Thanks a heap Pierre!

ADD REPLYlink written 9 months ago by andrew.phipps.000
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2147 users visited in the last hour