Error message using ReorderSam

JuttaJutta Bonn, GermanyMember

Dear GATK team,

I am still following the Best practice guide for SNP calling of RNA-seq data.
For one of my samples, I obtained the following error message when applying the following program code:

java -jar /opt/picard/picard.jar ReorderSam I=sample_dedupped.bam O=sample_reordered.bam R=Zea_mays.AGPv3.31.dna.genome.edit.fa CREATE_INDEX=TRUE

For all my other samples this program code worked fine, and I treaded all the samples similar in any previous steps.

INFO 2016-12-19 16:03:33 ReorderSam Processing 7
INFO 2016-12-19 16:04:03 ReorderSam Wrote 2554799 reads
INFO 2016-12-19 16:04:03 ReorderSam Processing 5
INFO 2016-12-19 16:04:44 ReorderSam Wrote 3838774 reads
INFO 2016-12-19 16:04:44 ReorderSam Processing 3
[Mon Dec 19 16:04:47 CET 2016] picard.sam.ReorderSam done. Elapsed time: 3,19 minutes.
To get help, see
Exception in thread "main" htsjdk.samtools.SAMFormatException: SAM validation error: ERROR: Record 228651, Read name K00114:299:H5Y5TBBXX:6:1202:27407:1508, bin field of BAM record does not equal value computed based on alignment start and end, and length of sequence to which read is aligned
at htsjdk.samtools.SAMUtils.processValidationErrors(
at htsjdk.samtools.BAMFileReader$BAMFileIterator.advance(
at htsjdk.samtools.BAMFileReader$
at htsjdk.samtools.BAMFileReader$
at htsjdk.samtools.BAMFileReader$BAMQueryFilteringIterator.advance(
at htsjdk.samtools.BAMFileReader$
at htsjdk.samtools.BAMFileReader$
at htsjdk.samtools.SamReader$
at htsjdk.samtools.SamReader$
at picard.sam.ReorderSam.writeReads(
at picard.sam.ReorderSam.doWork(
at picard.cmdline.CommandLineProgram.instanceMain(
at picard.cmdline.PicardCommandLine.instanceMain(
at picard.cmdline.PicardCommandLine.main(

I had a look at several websites, but I couldn't find any solution for this problem.
So, I would be very happy, if there are any suggestions how to solve this error?

Greetings Jutta

Best Answers


Sign In or Register to comment.