Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BWA crashes in Fedora

    Lately I changed OS from Ubuntu to Fedora (KDE) (both 64bit).
    While in Ubuntu BWA worked perfectly it now crashes while trying to align a 7.4GB Illumina fastQ file which I could align perfectly in Ubuntu.

    The output is

    Code:
    [bwa_aln] 17bp reads: max_diff = 2
    [bwa_aln] 38bp reads: max_diff = 3
    [bwa_aln] 64bp reads: max_diff = 4
    [bwa_aln] 93bp reads: max_diff = 5
    [bwa_aln] 124bp reads: max_diff = 6
    [bwa_aln] 157bp reads: max_diff = 7
    [bwa_aln] 190bp reads: max_diff = 8
    [bwa_aln] 225bp reads: max_diff = 9
    [bwa_aln_core] calculate SA coordinate... Segmentation fault (core dumped)
    I use 16GB quadcore 64bit Fedora 14 and I tried every bwa version from 0.5.6 to 0.5.9.

    The command was:
    Code:
    bwa aln ~/RefSeq/hg19_070510/Chromosomes/hg19.fa test.fq > test.aln
    Any help appreciated (looking in your direction lh3...)

  • #2
    I had the same problem (BWA on Fedora 13) and gave up on it for a long time. Then someone passed me a BLOG link:

    Here's what I use for bwa alignment (without removing PCR dups). You can replace the paths with your own and put into a bash script for au...


    My problem seemed to be i was using a solid2fastq script from a different program (BFAST i think. Why was i doing this, because it is C-compiled therefore much faster!). Once i used solid2fast from the BWA package it started working.

    I hope this helps!

    Comment


    • #3
      Thank you for your help.
      I was just too stupid for reading the manual. I did new indexing of the genome with the -a is argument which is not recommended for large genomes. After applying the -a bwt option in indexing, alignment worked.

      Comment


      • #4
        Dear all,

        I got this error message, too, when running bwa samse (bwa-0.5.9rc1):

        Code:
        [bwa_aln_core] refine gapped alignments... Segmentation fault (core dumped)
        however, the difference to the previous posts is that
        1. my original SOLiD reads were reformatted using bwa's solid2fastq.pl alright, and
        2. I ran a batch of 12 samples that were processed all the same way. Only 7 gave this error, 5 were processed without error.

        Any ideas why this happens will be appreciated.

        thanks,
        Sophia

        Comment


        • #5
          Let me reply to my own question:

          What happened is that the solid2fastq script obviously did not process the reads completely (as a part of the batch process).
          Funnily enough, the last read of the fastq file being incomplete did not prevent bwa aln to run and generate a .sai file.
          Only in the bwa samse step, this caused trouble and led to said error message.

          Hope this helps.

          cheers,
          Sophia

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Non-Coding RNA Research and Technologies
            by seqadmin




            Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

            Nobel Prize for MicroRNA Discovery
            This week,...
            10-07-2024, 08:07 AM
          • seqadmin
            Recent Developments in Metagenomics
            by seqadmin





            Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
            09-23-2024, 06:35 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Today, 02:44 PM
          0 responses
          7 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-11-2024, 06:55 AM
          0 responses
          14 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-02-2024, 04:51 AM
          0 responses
          110 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-01-2024, 07:10 AM
          0 responses
          116 views
          0 likes
          Last Post seqadmin  
          Working...
          X