Question: Pacbio Metagenome alignment
gravatar for 934963534
20 months ago by
93496353420 wrote:

Can I use bwa to align metagenome pacbio reads to several reference genome? Or are there any other tools to achieve metagenome alignment for pacbio reads?

alignment sequence • 642 views
ADD COMMENTlink modified 17 months ago by drkennetz440 • written 20 months ago by 93496353420

If you already know the genome(s) that could be present in your sample then you can use minimap2 to align reads to the genome(s).

Another approach would be to carry out denovo assembly using canu and then search the contigs against reference databases using either BLAST or DIAMOND.

ADD REPLYlink modified 17 months ago • written 17 months ago by Sej Modha4.5k
gravatar for drkennetz
17 months ago by
drkennetz440 wrote:

You don't use bwa to run against multiple genomes at once, but you could of course manually align your sequencing data to multiple genomes 1 by 1 using bwa mem or minimap2. Furthermore, for actual metagenome alignment tools, I would look at this paper:

As your question is more complicated than a simple yes or no. To do what you are asking properly there are several things to take into account (which are all discussed in this review).

You will want to look at the readQC metrics as well as the metagenome alignment section (and maybe even metagenome assembly).

ADD COMMENTlink written 17 months ago by drkennetz440
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1173 users visited in the last hour