Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • lfaller
    Junior Member
    • Jul 2017
    • 4

    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)
  • neavemj
    Member
    • Feb 2014
    • 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
      New Genomics Tools and Methods Shared at AGBT 2025
      by seqadmin


      This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

      The Headliner
      The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
      03-03-2025, 01:39 PM
    • seqadmin
      Investigating the Gut Microbiome Through Diet and Spatial Biology
      by seqadmin




      The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
      02-24-2025, 06:31 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 05:03 AM
    0 responses
    16 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-19-2025, 07:27 AM
    0 responses
    15 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-18-2025, 12:50 PM
    0 responses
    16 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-03-2025, 01:15 PM
    0 responses
    185 views
    0 reactions
    Last Post seqadmin  
    Working...