Bwa Indexing
1
0
Entering edit mode
11.9 years ago
Kssr ▴ 110

I am trying to index a reference genome using bwa index.

bwa index -a bwtsw xxx.fasta

It runs without errors and generates the following files with extensions .amb,.ann,.bwt,.pac,.sa.

I have never really done indexing myself before using bwa.However, I remember indexing generates other files as well such as .rpac,.rbwt,.rpac,.rsa.What does r indicate here? What are the files generated usually when we index a reference genome ? Any help is appreciated.

alignment • 9.3k views
ADD COMMENT
2
Entering edit mode

0.5.x generates .r files, but 0.6.x does not.

ADD REPLY
0
Entering edit mode

Hi there,

I cannot answer your question regarding what are these files are used for ? But I can tell that you have index your reference genome successfully.

ADD REPLY
0
Entering edit mode
ADD REPLY
0
Entering edit mode

hi to all.I have a problem during index reference file. it stops and doesn't continue after 2 hours. I should mention that I'm doing the whole process on VMware on my windows 10. does it matter?

ADD REPLY
0
Entering edit mode

I removed the previous discussion. Bottom line was that you do not have enough memory (RAM). Consider getting more RAM or try to switch to a more powerful computer. In case the issue persists please open a new question here on Biostars instead of warming up old threads. Please make sure to check on RAM consumption before opening a new question.

ADD REPLY
0
Entering edit mode

ok. I can not use another computer. do u think using original Linux instead of VM, would fix my problem?

ADD REPLY
0
Entering edit mode
11.9 years ago

The .r* files are reverse sequence index, BW-transformed sequence, suffix-array index. I think it uses these for reads on the minus strand, but I'm not sure. See class pac_rev.

ADD COMMENT

Login before adding your answer.

Traffic: 2989 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6