BWA to map reads to distant genome
0
0
Entering edit mode
18 months ago
BioGeek ▴ 150

Hi Stars, I wanted to map my paired-end reads against genome of phylogenetically related but distant species. What would be the best parameters to map it with bwa mem ?

genome Genome Mappings bwa-mem parameters • 505 views
ADD COMMENT
0
Entering edit mode

With the limited amount of information, this is difficult to answer (I admit, I don't know what I would need to give a good answer, anyway)

I would approach this in an iterative way, you can play with minimum seed length and bandwith, z-dropoff, then dive into match and mismatch scores and gap penalties. After all, much depends on how similar the genomes are. You will always align some xy percentage, but does that help you?

ADD REPLY
0
Entering edit mode

complementing Carambakaracho comment :

it's always advisable to start with the default settings and take it from there.

ADD REPLY
0
Entering edit mode

In principle, it's a bad idea to map to a different genome, but as suggested, start with the defaults and see how much can be used.

ADD REPLY

Login before adding your answer.

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