Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • seq_GA
    Senior Member
    • Feb 2009
    • 124

    picard error

    While I am trying to remove duplicates for pair end data using picard,
    I am getting th following errors.

    I start picard as below:
    Code:
    java -jar /tools/SAM/picard-tools-1.13/MarkDuplicates.jar I=bin1_sorted.bam O=bin1_PCR_removed.bam M=metrics.txt ASSUME_SORTED=true REMOVE_DUPLICATES=true VALIDATION_STRINGENCY=SILENT
    Runtime.totalMemory()=1031405568
    Exception in thread "main" java.lang.IllegalStateException: Invalid values for pairs and unique pairs: 1574422, 2050
    at net.sf.picard.sam.DuplicationMetrics.estimateLibrarySize(DuplicationMetrics.java:90)
    at net.sf.picard.sam.DuplicationMetrics.calculateDerivedMetrics(DuplicationMetrics.java:64)
    at net.sf.picard.sam.MarkDuplicates.doWork(MarkDuplicates.java:176)
    at net.sf.picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:142)
    at net.sf.picard.sam.MarkDuplicates.main(MarkDuplicates.java:79)


    Thanks.
  • seq_GA
    Senior Member
    • Feb 2009
    • 124

    #2
    Can anyone tell me why this is happening?

    Comment

    • Lv Ray
      Member
      • Jun 2014
      • 42

      #3
      try :
      samtools view -H bin1_sorted.bam |less -S
      inspect the header of your .bam file. I got the same error.

      Comment

      Latest Articles

      Collapse

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 06-05-2026, 10:09 AM
      0 responses
      14 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-04-2026, 08:59 AM
      0 responses
      24 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-02-2026, 12:03 PM
      0 responses
      28 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-02-2026, 11:40 AM
      0 responses
      23 views
      0 reactions
      Last Post SEQadmin2  
      Working...