This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BWA error [bns_restore_core] Parse error reading

Hi,

I have downloaded some bacteria genomes on ncbi and then I cat them.

I have built an index with these genome after trying to map my reads I got this error:

[bns_restore_core] Parse error reading db/contaminants.amb

Someone know what is going on ? Thanks

bwa index
• 10,659 views • link
updated by said3427 said3427 • written by Picasa Picasa

I have the same problem with bwa-0.7.12 (I tried 0.7.4 and it works). Did you find a solution? my genome file does not contain spaces.

• 0 views • link
written by lironyoffe • lironyoffe

3 answers

Based on previous encounters with such issues, please check if your genome file contains any spaces. Remove using sed or tr (if you find any).

• 0 views • link
written by lakhujanivijay lakhujanivijay

sed -i 's/\s*$//g' db/contaminants

bwa index db/contaminants

• 0 views • link
written by said3427 said3427

If you have spaces in your header, use this instead:

sed -i 's/^[^>]\s*$//g' db/contaminants
• 1 views • link
written by bird77 • bird77

See the answer from Vijay Lakhujani above. More specifically, remove whitespace from the fasta reference sequences (whitespace in fasta header is okay). See also: Burrows-Wheeler Aligner / Mailing Lists

https://sourceforge.net/p/bio-bwa/mailman/message/33230220/

https://sourceforge.net/p/bio-bwa/mailman/message/34430833/

• 0 views • link
written by tshtatland tshtatland

Log in to answer this question.

More posts like this