Error running Picard Tools v2.10.9

Hello- I'm trying to run MarkDuplicates in Picard Tools v2.10.9. Java version is 1.8.0_144. The data is a .bam file that has been sorted and indexed in samtools. (exome capture sequences from HiSeq 4000).

I used the following command line call to try the run:
java -Xmx5g -jar /Users/paa9/Desktop/picard.jar MarkDuplicates REMOVE_DUPLICATES=true

At first it looks like it is running successfully, but then I get a series of error messages:
INFO 2017-08-14 16:03:59 MarkDuplicates Start of doWork freeMemory: 1015051672; totalMemory: 1029177344; maxMemory: 4772593664
INFO 2017-08-14 16:03:59 MarkDuplicates Reading input file and constructing read end information.
INFO 2017-08-14 16:03:59 MarkDuplicates Will retain up to 17292006 data points before spilling to disk.
INFO 2017-08-14 16:04:08 MarkDuplicates Read 1,000,000 records. Elapsed time: 00:00:07s. Time for last 1,000,000: 7s. Last read position: JXUM01S000636:264,164
INFO 2017-08-14 16:04:08 MarkDuplicates Tracking 30491 as yet unmatched pairs. 4430 records in RAM.
INFO 2017-08-14 16:04:21 MarkDuplicates Read 2,000,000 records. Elapsed time: 00:00:20s. Time for last 1,000,000: 12s. Last read position: JXUM01S001819:138,692
INFO 2017-08-14 16:04:21 MarkDuplicates Tracking 41344 as yet unmatched pairs. 5 records in RAM.
INFO 2017-08-14 16:04:36 MarkDuplicates Read 3,000,000 records. Elapsed time: 00:00:36s. Time for last 1,000,000: 15s. Last read position: JXUM01S004018:43,724
INFO 2017-08-14 16:04:36 MarkDuplicates Tracking 48378 as yet unmatched pairs. 9 records in RAM.

Them jere is the error message I get…

To get help, see
Exception in thread "main" htsjdk.samtools.SAMException: /var/folders/t4/1tm2l_xd5r9c9lwhrjsfzwb80000gn/T/CSPI.4085908425341095212.tmp/4388.tmpnot found
at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(
at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(
at htsjdk.samtools.util.ResourceLimitedMap.get(
at htsjdk.samtools.CoordinateSortedPairInfoMap.getOutputStreamForSequence(
at htsjdk.samtools.CoordinateSortedPairInfoMap.ensureSequenceLoaded(
at htsjdk.samtools.CoordinateSortedPairInfoMap.remove(
at picard.sam.markduplicates.util.DiskBasedReadEndsForMarkDuplicatesMap.remove(
at picard.sam.markduplicates.MarkDuplicates.buildSortedReadEndLists(
at picard.sam.markduplicates.MarkDuplicates.doWork(
at picard.cmdline.CommandLineProgram.instanceMain(
at picard.cmdline.PicardCommandLine.instanceMain(
at picard.cmdline.PicardCommandLine.main(
Caused by: /var/folders/t4/1tm2l_xd5r9c9lwhrjsfzwb80000gn/T/CSPI.4085908425341095212.tmp/4388.tmp (Too many open files in system)
at Method)
at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(
... 11 more

To check and see if it is something wrong with the input file, I ran the “ValidateSamFile” in Picard Tools using the following call:

java -jar /Users/paa9/Desktop/picard.jar ValidateSamFile MODE=SUMMARY

And here is the terminal output I get:

HISTOGRAM java.lang.String

Error Type Count

Thank you for any suggestions!!!

Best Answers


Sign In or Register to comment.