bwa: fail to locate the index files
1
1
Entering edit mode
2.2 years ago

Hi all,

I'm trying to align a fastq file to a reference fasta file using bwa mem and for the life of me I keep getting the following error message. Any thoughts? Thanks!

bwa index -a bwtsw myref.fasta -p myref
ls
myref.amb  myref.bwt  myref.pac
myref.ann  myref.fasta    myref.sa
  
bwa mem -M 1 -R '@RG\tID:{foo}.{bar}\tPU:{foo}.{bar}.{999999}\tSM:{999999}\tPL:{foo}\tLB{bar}' myref myfastq.fastq.gz

[E::bwa_idx_load_from_disk] fail to locate the index files
alignment bwa • 6.2k views
ADD COMMENT
0
Entering edit mode

Sorry about the weird formatting; this is my first post to Biostars.

ADD REPLY
0
Entering edit mode

Did you try myref.fasta instead of just myref?

ADD REPLY
2
Entering edit mode
2.2 years ago
h.mon 32k

I am not sure, but I think the cause of the error is the -M parameter doesn't take any values:

         -M        Mark shorter split hits as secondary (for Picard compatibility).
  

Try with:

bwa mem -M -R '@RG\tID:{foo}.{bar}\tPU:{foo}.{bar}.{999999}\tSM:{999999}\tPL:{foo}\tLB{bar}' myref myfastq.fastq.gz
ADD COMMENT
0
Entering edit mode

Yup, that was it! Thank you!!!

ADD REPLY
0
Entering edit mode

Nice, you can upvote (thumb up) an answer that was useful, or even accept that answer.

ADD REPLY

Login before adding your answer.

Traffic: 1159 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