bbmap, mapq and uniquelly mapped reads filtering
1
0
Entering edit mode
3 months ago
boczniak767 ▴ 790

Hi,

I'd like to know, how do postalignment filtering to get only uniquelly mapped reads. I know that it is possible at the bbmap level, i.e. during mapping. I've read, that for aligners like BWA there are correspondence between mapq value and uniqueness of alignment.

  1. Does for bbmap similar correspondence exists?
  2. If yes, what mapq threshold should be used to get only unique mappings?
  3. Is there any general method (no matter what aligner was used) to get only uniquelly mapped reads from sam file, for example using flags and samtools?
alignment mapq • 223 views
ADD COMMENT
0
Entering edit mode
3 months ago
GenoMax 104k

Many threads exist that show how to extract uniquely mapped reads. Here are a couple of examples:
Counting and extracting in a new file the unique mapped reads from BWA 0.7.15-r1140 mapping
How To Filter Mapped Reads With Samtools

Implementation of mapq values is not consistent. What is consistent is multi-mappers will have 0 or low values.

If you are using ambig= option with bbmap.sh to choose a random location

random  (select one top-scoring site randomly)

then you are selecting a site randomly during mapping step itself.

ADD COMMENT

Login before adding your answer.

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