Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bwa sampe error ... bad database format?

    Hi,
    I am running bwa 0.5.9 on a custom, indexed, sequence database comprised of various genbank sequences.

    When I run:

    bwa aln gb.fa x_1.fq > a.sai
    bwa aln gb.fa x_2.fq > b.sai
    bwa sampe gb.fa a.sai b.sai x_1.fq x_2.fq > gb.sam

    bwa crashes on bwa sampe with:

    [bns_restore_core] fail to open file ' VirtualSequence=""></div><div id="viewercontent13" class="seq gbff" val="381286149" SequenceSize="29462" VirtualSequence=""></div><div id="viewercontent14" class="seq gbff" val="381283770" SequenceSize="35021" VirtualSequence=""></div><div id="viewercontent15" class="seq gbff" val="381286161" SequenceSize="29364" VirtualSequence=""></div><div id="viewercontent16" class="seq gbff" val="381282985" SequenceSize="9499" VirtualSequence=""></div><div id="viewercontent17" class="seq gbff" val="381216220" SequenceSize="16093" VirtualSequence=""></div><div id="viewercontent18" class="seq gbff" val="381147283" SequenceSize="5929" VirtualSequence=""></div><div id="viewercontent19" class="seq gbff" val="380307325" SequenceSize="10631" VirtualSequence=""></div><div id="viewercontent20" class="seq gbff" val="380862749" SequenceSize="6588" VirtualSequence=""></div>'. Abort!
    ./test.sh: line 3: 25171 Aborted bwa sampe gb.fa a.sai b.sai x_1.fq x_2.fq > gb.sam

    The 9GB gb.fa database file was indexed with default, as per the instructions provided and no errors were returned from the indexing.

    When I run the same sets of bwa commands using a different, human genome database, there are no errors.

    I was wondering whether characters others than ACGTN in the seq. db would make bwa sampe crash (there are such characters in the gb.fa)? Has anyone seen an error like this?

    Thanks for your help
    Rene

  • #2
    That block looks like XML? Is that really a FA file?

    Comment


    • #3
      it is, now reformatting to remove any offending characters.
      The interesting thing about this error is that, usually, when bwa returns a fail to open file rror, the file is in quotation, not some XML block:

      eg.
      [bns_restore_core] fail to open file '/projects/02/rwarren/nt_db/2012/humanNCRNA.fa.fa.ann'. Abort!
      ./testdb.sh: line 3: 30215 Aborted /home/pubseq/BioSw/bwa/bwa-0.5.9/bwa sampe /projects/02/rwarren/nt_db/2012/humanNCRNA.fa.fa 1.sai 2.sai unmappedSorted.bam unmappedSorted.bam > test.sam

      bizarre.

      Comment


      • #4
        i confirm that junk in the custom sequence database generated the error.
        issue resolved.

        Comment


        • #5
          I just did a test, I modified a .fasta with two illegal char 'R', 'S', and then run INDEX, ALN, and SAMPE, everything still works, although output SAM files are different. So BWA should be able to cope with illegal characters.

          And if you open the .amb file, you should see the illegal chars listed there, in my case:

          249240584 1 R
          249240585 1 S

          So if your fasta is big, and really messed with many illegal chars, then your .amb file should be very very big.

          Besides, your error msg is BWA complaining a file can't be opened, at this stage, .bwt, .sa, and .pac are not opened yet, so having illegal chars shouldn't be the reason for this error.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Essential Discoveries and Tools in Epitranscriptomics
            by seqadmin


            The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
            Today, 07:01 AM
          • seqadmin
            Current Approaches to Protein Sequencing
            by seqadmin


            Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
            04-04-2024, 04:25 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 04-11-2024, 12:08 PM
          0 responses
          37 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 10:19 PM
          0 responses
          39 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 09:21 AM
          0 responses
          35 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-04-2024, 09:00 AM
          0 responses
          54 views
          0 likes
          Last Post seqadmin  
          Working...
          X