Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • ArtVandelay
    Junior Member
    • Sep 2017
    • 2

    #1

    IndexOutOfBoundsException in BBTool's kmercountexact.sh

    I'm trying to use BBTool's kmercountexact.sh (version 37.50) to produce a sketch from a test fasta file:
    kmercountexact.sh in=test1.fasta sketch=test1.sketch mincount=3 -overwrite=true

    and I'm getting an IndexOutOfBoundsException. Does anyone have any idea of what I can do differently or how I may otherwise get help?

    Here is the full output:

    > kmercountexact.sh in=test1.fasta sketch=test1.sketch mincount=3 -overwrite=true
    Max memory cannot be determined. Attempting to use 3200 MB.
    If this fails, please add the -Xmx flag (e.g. -Xmx24g) to your command,
    or run this program qsubbed or from a qlogin session on Genepool, or set ulimit to an appropriate value.
    java -ea -Xmx3200m -Xms3200m -cp /Users/droe/bin/bbmap/current/ jgi.KmerCountExact in=test1.fasta sketch=test1.sketch mincount=3 -overwrite=true
    Executing jgi.KmerCountExact [in=test1.fasta, sketch=test1.sketch, mincount=3, -overwrite=true]

    Initial:
    Memory: max=3215m, free=3182m, used=33m

    Executing kmer.KmerTableSet [in=test1.fasta, sketch=test1.sketch, mincount=3, -overwrite=true]

    Initial:
    Ways=31, initialSize=128000, prefilter=f, prealloc=f
    Memory: max=3215m, free=3182m, used=33m

    Estimated kmer capacity: 121523321
    After table allocation:
    Memory: max=3215m, free=3111m, used=104m

    After loading:
    Memory: max=3215m, free=3060m, used=155m

    Input: 1 reads 1000 bases.

    For K=31
    Unique Kmers: 970
    Load Time: 0.059 seconds.
    Generating sketch.
    Exception in thread "main" java.lang.IndexOutOfBoundsException: Index: 0, Size: 0
    at java.util.ArrayList.rangeCheck(ArrayList.java:653)
    at java.util.ArrayList.get(ArrayList.java:429)
    at sketch.SketchTool.toSketch(SketchTool.java:327)
    at sketch.SketchTool.toSketch_MT(SketchTool.java:225)
    at sketch.SketchTool.toSketch(SketchTool.java:167)
    at jgi.KmerCountExact.makeSketch(KmerCountExact.java:360)
    at jgi.KmerCountExact.process(KmerCountExact.java:264)
    at jgi.KmerCountExact.main(KmerCountExact.java:58)
    Attached Files
  • Brian Bushnell
    Super Moderator
    • Jan 2014
    • 2709

    #2
    I just tested and was able to replicate this... it's due to the "mincount=3" flag. That file has no kmers that occur at least 3 times, so it ends up trying to make a size-0 sketch, which isn't allowed. If you take away "mincount=3" it works fine. I'll clarify that error message.

    Comment

    • ArtVandelay
      Junior Member
      • Sep 2017
      • 2

      #3
      Makes sense. Thank you!

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        New Genomics Technologies Take Aim at Long-Standing Limits
        by SEQadmin2


        Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

        We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
        ...
        Yesterday, 10:25 AM
      • SEQadmin2
        How Immunogenomics Decodes Immunity’s Genetic Blueprint
        by SEQadmin2




        The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.

        This convergence of genetics, immunology, and computation...
        09-01-2026, 05:41 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, Today, 09:51 AM
      0 responses
      7 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 09-25-2026, 09:06 AM
      0 responses
      31 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 09-23-2026, 11:05 AM
      0 responses
      27 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 09-18-2026, 11:37 AM
      1 response
      47 views
      0 reactions
      Last Post pekgio
      by pekgio
       
      Working...