Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: More than one possible mapping place for a read
0
gravatar for biotech
3.0 years ago by
biotech • 520
United States
biotech • 520 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 COMMENT • link •
modified 3.0 years ago by WouterDeCoster ♦ 38k • written 3.0 years ago by biotech • 520
0
gravatar for WouterDeCoster
3.0 years ago by
WouterDeCoster ♦ 38k
Belgium
WouterDeCoster ♦ 38k 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 COMMENT • link written 3.0 years ago by WouterDeCoster ♦ 38k

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 REPLY • link written 3.0 years ago by genomax ♦ 65k
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 REPLY • link written 3.0 years ago by WouterDeCoster ♦ 38k

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 REPLY • link written 3.0 years ago by biotech • 520

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

ADD REPLY • link written 3.0 years ago by biotech • 520

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

ADD REPLY • link written 3.0 years ago by biotech • 520
Please log in to add an answer.

Similar posts • Search »

Content
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