Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Error encountered when running kmercountexact.sh from bbmap

    Hi all,

    I am encountering a java.lang.AssertionError when trying to run kmercountexact.sh.

    The two input files are 820 Mb each. I also tried running this script with subsetted data of the original datasets (using just the first 250 fastq records) but I got the same error message.

    The help documentation says "Last modified August 2, 2017" so it is a fairly recent version.

    Any advice is appreciated!
    ~Lina

    Code:
    $ ~/bbmap/kmercountexact.sh in1=${R1} in2=${R1} khist=${HIST} peaks=${PEAKS}
    java -ea -Xmx89748m -Xms89748m -cp /home/lina/bbmap/current/ jgi.KmerCountExact in1=7064890_R1.fastq in2=7064890_R1.fastq khist=khist.txt peaks=peaks.txt
    Executing jgi.KmerCountExact [in1=7064890_R1.fastq, in2=7064890_R1.fastq, khist=khist.txt, peaks=peaks.txt]
    
    Initial:
    Memory: max=90186m, free=88774m, used=1412m
    
    Executing kmer.KmerTableSet [in1=7064890_R1.fastq, in2=7064890_R1.fastq, khist=khist.txt, peaks=peaks.txt]
    
    Initial:
    Ways=101, initialSize=128000, prefilter=f, prealloc=f
    Memory: max=90186m, free=87363m, used=2823m
    
    Estimated kmer capacity:        3509022984
    After table allocation:
    Memory: max=90186m, free=86421m, used=3765m
    
    Exception in thread "main" java.lang.AssertionError: 7064890_R1.fastq, 7064890_R1.fastq
            at stream.ConcurrentReadInputStream.getReadInputStream(ConcurrentReadInputStream.java:111)
            at stream.ConcurrentReadInputStream.getReadInputStream(ConcurrentReadInputStream.java:45)
            at kmer.KmerTableSet.loadKmers(KmerTableSet.java:381)
            at kmer.AbstractKmerTableSet.loadKmers(AbstractKmerTableSet.java:288)
            at kmer.AbstractKmerTableSet.processInput(AbstractKmerTableSet.java:135)
            at kmer.AbstractKmerTableSet.process(AbstractKmerTableSet.java:80)
            at jgi.KmerCountExact.process2(KmerCountExact.java:285)
            at jgi.KmerCountExact.process(KmerCountExact.java:259)
            at jgi.KmerCountExact.main(KmerCountExact.java:58)

  • #2
    Hi Lina,

    I think you might be trying to give your forward reads twice to both the in1 and in2 arguments. Maybe your in2 argument should be:

    in2=${R2}

    I guess bbmap checks that the F and R reads have different file names and this is why you get an 'assertion' error.

    Cheers,

    Matt.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    67 views
    0 likes
    Last Post seqadmin  
    Working...
    X