BWA not creating rBWT file

[ramiro2k]

Hi everyone,

I have BWA version 0.6.1 and for some reason it is not creating the .rbwt file in the human genome:

bash-3.2$ ~/MPS/programs/bwa-0.6.1/bwa index -a bwtsw genome.fa

creates:

-rwxr-xr-x 1 rbarrant usr 3142044949 Mar 10 14:45 genome.fa
-rw-r--r-- 1 rbarrant usr 6144 Mar 10 16:26 genome.fa.amb
-rw-r--r-- 1 rbarrant usr 941 Mar 10 16:26 genome.fa.ann
-rw-r--r-- 1 rbarrant usr 3080436148 Mar 10 16:25 genome.fa.bwt
-rw-r--r-- 1 rbarrant usr 770109014 Mar 10 16:26 genome.fa.pac
-rw-r--r-- 1 rbarrant usr 1540218080 Mar 10 17:01 genome.fa.sa


any suggestions?

Thanks,
Ramiro

[nilshomer]

See the "NEWS" file that comes with BWA:

...the index data structure has been changed to support genomes longer than 4GB...

[aituka]

I ramiro2k,
I do have the same problem. In fact I don't care if there is or not rbwt file. the only problem is that without it seems I can't really continue with the gatk pipeline. Did you find something.

did I miss something?