Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bowtie memory warning

    Hi

    I'm mapping solid reads using bowtie and received the following warnings for MANY reads. Mapping % is very low as well.
    Am I running out of memory?

    command line:
    bowtie -C -S -k 3 --best --un unmapped_list genome_index -f seq.csfasta -Q seq.qual bowtie.sam



    Warning: Exhausted best-first chunk memory for read 301_728_1008_F3 (patid 11858997); skipping read

    # reads processed: 19323816
    # reads with at least one reported alignment: 4408352 (22.81%)
    # reads that failed to align: 14915464 (77.19%)
    Reported 5326165 alignments to 1 output stream(s)

    Thanks

    charles

  • #2
    I get similar warnings when I use paired end data. But if I align read1 and read2 separately, I do not get any warnings. Solutions/pointers are appreciated.
    Thanks,
    Vishal

    Comment


    • #3
      Yes, I believe you are running out of memory. If you use --chunkmbs 200 you should stop getting the warning. This may help improve your alignment statistics.

      Comment


      • #4
        In fact, I have tried using --chunkmbs 1000 and still managed to get these errors ... only for paired-end data though.

        Comment


        • #5
          Try higher I guess. 200 worked for me when working with non-paired end data.

          Comment


          • #6
            bowtie memory warning

            All

            Thanks for the advice:

            Using --chunkmbs 200 eliminated the memory warnings.

            I'm still not mapping a significant fraction of the reads. I'm going to trim the unmapped colorspace reads from 3'end to see if I can get any of these to map - other suggestions?

            bowtie -C -S -k 3 --best


            # reads processed: 19323816
            # reads with at least one reported alignment: 4408352 (22.81%)
            # reads that failed to align: 14915464 (77.19%)
            Reported 5326165 alignments to 1 output stream(s)

            Comment


            • #7
              Bowtie: Memory error warning

              Same as the quoted question. When I run the Bowtie command as below to generate SAM with the 1 best hit from fasta input, it has such warning messages,

              $ bowtie -f -S -m 1 -k 1 --best h_sapiens_37_asm Input.fa > Output.sam
              ------------------
              Warning: Exhausted best-first chunk memory for read seq_5023_1 (patid 5023); skipping read
              Warning: Exhausted best-first chunk memory for read seq_5041_1 (patid 5041); skipping read
              Warning: Exhausted best-first chunk memory for read seq_20071_1 (patid 20071); skipping read
              ......
              ------------------

              Which flag options should I use?

              Many thanks.

              Ben


              Originally posted by crh View Post
              Hi

              I'm mapping solid reads using bowtie and received the following warnings for MANY reads. Mapping % is very low as well.
              Am I running out of memory?

              command line:
              bowtie -C -S -k 3 --best --un unmapped_list genome_index -f seq.csfasta -Q seq.qual bowtie.sam



              Warning: Exhausted best-first chunk memory for read 301_728_1008_F3 (patid 11858997); skipping read

              # reads processed: 19323816
              # reads with at least one reported alignment: 4408352 (22.81%)
              # reads that failed to align: 14915464 (77.19%)
              Reported 5326165 alignments to 1 output stream(s)

              Thanks

              charles

              Comment


              • #8
                I had this problem as well. I tried --chunkmbs 200 as well, and it didn't help, but --chunkmbs 2000 seems like it might have fixed the bug being reported, but I am still getting very low alignmment proportions (5% of reads).

                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, 03-27-2024, 06:37 PM
                0 responses
                12 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-27-2024, 06:07 PM
                0 responses
                11 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                52 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                68 views
                0 likes
                Last Post seqadmin  
                Working...
                X