Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • xied75
    replied
    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.

    Leave a comment:


  • rwarren
    replied
    i confirm that junk in the custom sequence database generated the error.
    issue resolved.

    Leave a comment:


  • rwarren
    replied
    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.

    Leave a comment:


  • xied75
    replied
    That block looks like XML? Is that really a FA file?

    Leave a comment:


  • rwarren
    started a topic bwa sampe error ... bad database format?

    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

Latest Articles

Collapse

  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin


    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM
  • 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...
    04-22-2024, 07:01 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
25 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
45 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-09-2024, 02:46 PM
0 responses
59 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-07-2024, 06:57 AM
0 responses
45 views
0 likes
Last Post seqadmin  
Working...
X