This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MarkDuplicates error for RNA-Seq

Hi, I am trying to run MarkDuplicates but keep getting the error listed below, although I added the TMP_DIR=yourdirectory/tmp in command line.

Please advice Thanks a bunch!

**Command**
 *java -jar MarkDuplicates.jar INPUT= Read_sorted.bam OUTPUT=Read_sortremove.bam METRICS_FILE=metrics.txt REMOVE_DUPLICATES=true VALIDATION_STRINGENCY=LENIENT TMP_DIR=yourdirectory/tmp*


**Error**
----------
Exception in thread "main" htsjdk.samtools.SAMException: /home/as/SNPiR/yourdirectory/tmp/CSPI.1714605993024191960.tmp/13316.tmpnot found
    at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(FileAppendStreamLRUCache.java:63)
    at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(FileAppendStreamLRUCache.java:49)
    at htsjdk.samtools.util.ResourceLimitedMap.get(ResourceLimitedMap.java:76)
    at htsjdk.samtools.CoordinateSortedPairInfoMap.getOutputStreamForSequence(CoordinateSortedPairInfoMap.java:180)
    at htsjdk.samtools.CoordinateSortedPairInfoMap.put(CoordinateSortedPairInfoMap.java:164)
    at picard.sam.DiskReadEndsMap.put(DiskReadEndsMap.java:67)
    at picard.sam.MarkDuplicates.buildSortedReadEndLists(MarkDuplicates.java:449)
    at picard.sam.MarkDuplicates.doWork(MarkDuplicates.java:177)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:183)
    at picard.sam.MarkDuplicates.main(MarkDuplicates.java:161)
Caused by: java.io.FileNotFoundException: /home/as/SNPiR/yourdirectory/tmp/CSPI.1714605993024191960.tmp/13316.tmp (Too many open files)
    at java.io.FileOutputStream.open0(Native Method)
    at java.io.FileOutputStream.open(FileOutputStream.java:270)
    at java.io.FileOutputStream.<init>(FileOutputStream.java:213)
    at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(FileAppendStreamLRUCache.java:60)
rna-seq software error
• 1,406 views • link
updated by GenoMax GenoMax • written by KG • KG

Normally there is no need to mark/remove duplicates in RNAseq data.

• 0 views • link
written by GenoMax GenoMax

Picard is not finding the temporary directory:

Exception in thread "main" htsjdk.samtools.SAMException: /home/as/SNPiR/yourdirectory/tmp/CSPI.1714605993024191960.tmp/13316.tmpnot found

But listen to genomax suggestion, probably there is no need to mark/remove duplicates in your data.

• 0 views • link
written by h.mon h.mon

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Picard MarkDuplicates shows error (removing pcr duplicates)
    written by ashaneev07 • ashaneev07
    Hiii.... i got the following error while running picards markduplicates. Does anyone have any experience with using this command in picard? Need help.. > java …
  • Error in picard tool: Mark Duplicates
    written by jcanitz • jcanitz
    Hey, I have a problem with the picard tool: Mark Duplicates To get an overview: I have a *.bam input file which is sorted. And …
  • error when running MarkDuplicates
    written by kk.mahsa kk.mahsa
    when i run Picard to mark duplicates i get this error: To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp Exception in thread "main" htsjdk.samtools.SAMException: /tmp/CSPI.61127736197136018.tmp/7237.tmpnot found at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(FileAppendStreamLRUCache.java:63) …