This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Newbie question: how to create BWA index from multiple .fasta files

I would like to use BWA MEM to align short reads against the entire hg19 human genome. To do that, I assume that I must create a BWA index. I understand how to create an index from a single file, but how does one create an index from multiple files? Is it best practice to just concatenate all files together? Or is there a better way?

Thanks in advance.

sequencing bwa index
• 6,734 views • link
updated by swbarnes2 swbarnes2 • written by phonybone • phonybone

Easy way may be to get a sequence/annotation/index bundle from the iGenomes site.

• 0 views • link
written by GenoMax GenoMax

2 answers

I'm assuming you currently have each chromosome as a separate fasta file. If you don't have a scientific or logistical reason to avoid concatenating the files, then the most straightforward approach would be to concatenate those fasta files into a single one, and then compute your index. Otherwise you'll need to merge alignment data downstream, and you could miss out on things like discordant alignments and large structural abnormalities such as chromosomal translocation.

• 0 views • link
written by Dan D Dan D

Instead of making a catted file like Dan suggests, you might be able to get away with catting them and piping that straight into bwa. Something like

cat *.fa | bwa index -p catted -
• 0 views • link
written by swbarnes2 swbarnes2

Piping is awesome. But wouldn't that preclude running subsequent alignments using that index, unless the complete catted fasta file was saved somewhere?

This would work, I think:

cat *.fa | tee bwa index -p catted - > catted.fasta

• 0 views • link
written by Dan D Dan D

Yes, this way does not make a catted fasta, but I think you can still use the index without it.

• 0 views • link
written by swbarnes2 swbarnes2

Log in to answer this question.

More posts like this