PICARD ValidateSamFile terminating prematurely, without errors, processes only portion of total reads
0
0
Entering edit mode
6.4 years ago
freuv ▴ 20

Hi,

I am processing ~200 million WGS aligned reads with Picard's ValidateSamFile, with the following options:

MODE=VERBOSE MAX_OUTPUT=10000000 TMP_DIR=[/local2/tmp] IGNORE_WARNINGS=false VALIDATE_INDEX=true IS_BISULFITE_SEQUENCED=false MAX_OPEN_TEMP_FILES=8000 VERBOSITY=INFO QUIET=false VALIDATION_STRINGENCY=STRICT COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false

However, the job ends prematurely without any errors and says it is complete -- after processing only 300K or so reads.

I have used this tool before without any issues for alignments with multimillions of reads.

Please advise, any help would be much appreciated!

freuv

sequencing Assembly genome • 1.6k views
ADD COMMENT
0
Entering edit mode

version of picard ? full command line ?

ADD REPLY
0
Entering edit mode

Picard 1.67

Command Line:

java -Xmx12G -Xms12G -jar ValidateSamFile.jar I=input.bam MODE=VERBOSE O=results.txt MAX_OUTPUT=10000000 TMP_DIR=/local2/tmp

ADD REPLY
0
Entering edit mode

your version of picard is outdated. 1.67 there is more than 500 bug fixes/commit since 1.67. Update picard and try again please. https://github.com/broadinstitute/picard/releases/tag/2.16.0

ADD REPLY
0
Entering edit mode

I am unable to install the latest version, as the cluster I work on has java v-1.7.0_151. Please advise.

Thanks! F

ADD REPLY
0
Entering edit mode

java 1..7 has been deprecated: http://www.oracle.com/technetwork/java/javase/downloads/java-se-jre-7-download-432155.html

This release is no longer the most current release of Java SE

update java to 1.8...

ADD REPLY

Login before adding your answer.

Traffic: 1522 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6