Question: Bwa Indexing
0
gravatar for Kssr
8.0 years ago by
Kssr110
Kssr110 wrote:

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 • 6.8k views
ADD COMMENTlink modified 23 months ago by Leila.javanparast0 • written 8.0 years ago by Kssr110
2

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

ADD REPLYlink written 8.0 years ago by lh332k

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 REPLYlink written 8.0 years ago by Houkto210

check http://seqanswers.com/forums/showthread.php?t=9808 this link for an answer.

ADD REPLYlink written 8.0 years ago by Houkto210

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 REPLYlink written 23 months ago by Leila.javanparast0

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 REPLYlink written 23 months ago by ATpoint44k

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

ADD REPLYlink written 22 months ago by Leila.javanparast0
0
gravatar for Matt Shirley
8.0 years ago by
Matt Shirley9.5k
Cambridge, MA
Matt Shirley9.5k wrote:

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 COMMENTlink modified 8.0 years ago • written 8.0 years ago by Matt Shirley9.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1553 users visited in the last hour
_