More than one possible mapping place for a read
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?
• 2,196 views
•
link
updated
by
WouterDeCoster
•
written
by
biotech
•
written
by
biotech
1 answer
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).
• 0 views
•
link
written
by
WouterDeCoster
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.
• 0 views
•
link
written
by
WouterDeCoster
Log in to answer this question.
More posts like this
-
bwa mem -T (alignment score) not doing anything
written by chris.bird •
-
bwa mem and multiply mapped reads
written by lfreeborn •
-
BWA-MEM output for two species
written by maheetha.b •
-
Meaning of BWA-MEM MAPQ values
written by Isaac C. Joseph •
-
bwa.kit: run-bwamem or bwa mem
written by iraun
-
Proper use of BWA-backtrack and BWA-MEM
written by cfarmeri
-
How can I best use -c in bwa mem for alignment and to find polymorphisms?
written by kirannbishwa01
-
RNA-seq reads mapping to paralogs with high MAPQ. How to ensure this?
written by biotech
-
When and why is bwa aln better then bwa mem?
written by dariober
-
Does The Mapping Of A Read'S Pair Affect Its Own Mapping Score In Bwa?
written by jhrf •
