BWA index error

[aldrinyim]

Dear all,

I never met this problem before not until I used a multiple fasta file (around 30) with each fasta sequence around 2000 to 3000 bp, is it not applicable to use short reads to map to fasta sequence of such kind?

the problem is

[bwa_index] Pack FASTA... 0.00 sec
[bwa_index] Reverse the packed sequence... 0.00 sec
[bwa_index] Construct BWT for the packed sequence...
[BWTIncConstructFromPacked] 10 iterations done. 18660 characters processed.
[BWTIncConstructFromPacked] 20 iterations done. 28804 characters processed.
[BWTIncConstructFromPacked] 30 iterations done. 32948 characters processed.
BWTIncSetBuildSizeAndTextAddr(): Not enough space allocated to continue construction!

Would anyone help?

[aldrinyim]

Oh hi everyone I just test the multiple fasta file and found that if we want an effective indexing of the fasta file, at least each fasta sequence should be >110000bp.

The program now works, thanks everyone.

[nilshomer]

Oh hi everyone I just test the multiple fasta file and found that if we want an effective indexing of the fasta file, at least each fasta sequence should be >110000bp.

The program now works, thanks everyone.

There are two algorithms to create indexes, one for short genomes/references, and one for long genomes/references. The default is the algorithm for long references so you have to specify "-a is" to run the algorithm for the short reference.

[aldrinyim]

Thanks nilshomer for reminding me of that, because i wrote a few scripts for BWA, I mistakenly used the script using dbwtsw for the short genome.

The problem should be fixed now!