bacteria genome mapping using bwa
0
0
Entering edit mode
20 months ago
hirakuda • 0

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?

sequence genome • 712 views
ADD COMMENT
1
Entering edit mode

What is the error?

First you need to create an index for your bacterial genome reference.

Then you can align your sample(s) to the reference genome.

bwa index ref.fa

bwa mem ref.fa reads.fq > aln-se.sam

samtools view -S -b aln-se.sam > aln-se.bam

Note that the output format is not .out

http://bio-bwa.sourceforge.net/bwa.shtml

http://quinlanlab.org/tutorials/samtools/samtools.html

ADD REPLY

Login before adding your answer.

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