Please someone suggest how to use bwa to map bacteria genomes to a sample? I want to use bwa mem. I did something bwa mem bacteria ref original sample > align.mem.out. and got some error. I dont know whether I am using it correctly. Can someone please suggest the correct command to use bwa to map severel bacteria together?