Question: More than one possible mapping place for a read
0
gravatar for biotech
3.0 years ago by
biotech520
United States
biotech520 wrote:

When there is more than one possible mapping place for a read, being both places of equal score,

Where 'bwa mem' mapper finally maps it?

What about the mapping quality scores for the read?

bwa • 757 views
ADD COMMENTlink modified 3.0 years ago by WouterDeCoster38k • written 3.0 years ago by biotech520
0
gravatar for WouterDeCoster
3.0 years ago by
Belgium
WouterDeCoster38k wrote:

It depends on your settings, but I think the default is multimapping the reads. So the same read will be present in multiple locations. Mapping quality will be 0 (meaning the read maps equally well in multiple places).

ADD COMMENTlink written 3.0 years ago by WouterDeCoster38k

Are you certain about multi-mapping? bwa mem manual entry does not say anything specific but the older sampe/samse explicitly say that

Repetitive read pairs will be placed randomly.

Repetitive hits will be randomly chosen.

ADD REPLYlink written 3.0 years ago by genomax65k
1

You are giving me doubts about my answer. ;-) Might very well be possible that random mapping happens.

Well, the easiest approach to settle this would be to test it. Have a single multimapping read, map it, count the reads in the samfile.

ADD REPLYlink written 3.0 years ago by WouterDeCoster38k

I really don't care about random mapping as long as having mapping quality of zero for them. In fact, multi-mapping reads should be discarded for some analysis such as RNA-seq.

ADD REPLYlink written 3.0 years ago by biotech520

Great proof test @decosterwouter, I will try to see what's going on!

ADD REPLYlink written 3.0 years ago by biotech520

I agree about that mapping quality is 0, but not sure if repetitive read pairs will be placed randomly.

ADD REPLYlink written 3.0 years ago by biotech520
Please log in to add an answer.

Help
Access

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