BWA aln - failed to locate the index
2
0
Entering edit mode
20 months ago
m.wekking ▴ 10

So I have alot of paired-end sequence samples in the formate sample_x_fw.fastq and sample_x_rv.fastq and want to align these with the BWA sampe command

so I downloaded the hg19.fasta and ran:

bwa index -p bwa_hg19 -a bwtsw hg19.fa

this resulted in:

  • bwa_hg19.amb
  • bwa_hg19.ann
  • bwa_hg19.bwt
  • bwa_hg19.pac
  • bwa_hg19.sa

after this I tried to run:

bwa aln bwa_hg19.fa sample_1_fw.fastq > aln_sa1.sai

and I get the error:

[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_aln] fail to locate the index

Does anyone have any idea what goes wrong?

bwa bioinformatics alignment • 2.0k views
ADD COMMENT
1
Entering edit mode

quick question : why bwa aln and not the more recent bwa mem ? Are your reads small ( < 50bp ? )

ADD REPLY
3
Entering edit mode
20 months ago
ATpoint 49k

It is bwa aln bwa_hg19 sample_1_fw.fastq > aln_sa1.sai as the prefix you gave for indexing is bwa_hg19, not bwa_hg19.fa

ADD COMMENT
0
Entering edit mode

That's it! thank you very much!

ADD REPLY
2
Entering edit mode
20 months ago

try

bwa aln bwa_hg19 sample_1_fw.fastq > aln_sa1.sai

or

bwa index -a bwtsw hg19.fa
bwa aln hg19.fa sample_1_fw.fastq > aln_sa1.sai
ADD COMMENT
0
Entering edit mode

That's it! thank you very much!

ADD REPLY

Login before adding your answer.

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