Question: How to identify bowtie2 unique mapping using sam flag value?
gravatar for jinxinhao1988
2.5 years ago by
jinxinhao198850 wrote:

I used BWA before,after mapping to genome or cDNA,I can filter unique mapping reads by sam flag value,but BWA mapping efficiency is not as good as bowtie2,so now I use bowtie2 instead.But I don't know the bowtie2 flag value that can indicate unique mapping.Could I get some answers?Thank you.

rna-seq sam bowtie2 flag • 1.5k views
ADD COMMENTlink modified 10 months ago by Biostar ♦♦ 20 • written 2.5 years ago by jinxinhao198850

Which definition of "unique" do you want to go by? I recommend that you just filter by an appropriate MAPQ.

ADD REPLYlink written 2.5 years ago by Devon Ryan81k
gravatar for d-cameron
20 months ago by
d-cameron1.8k wrote:

If you are running bowtie2 with the -k or -a options, then any read that has only a single output alignment can be considered uniquely mapped. If you're using the default reporting only the best alignment, the you should set a threshold on the MAPQ score reported based on how uniquely mapping you required your alignment.

ADD COMMENTlink written 20 months ago by d-cameron1.8k
gravatar for khalid.belkhir
20 months ago by
khalid.belkhir20 wrote:

Among the optional fields reported by bowtie2 you have :

XS:i:<n> aln score for next best aln (only reported if >1 aln found)

So the absence of this field for a read (in addition to the mapping flag and MAPQ) can be a sign of uniqueness ?

ADD COMMENTlink written 20 months ago by khalid.belkhir20

If your definition of "unique" is that there is only one alignment passing the --score-min threshold then yes.

ADD REPLYlink written 20 months ago by Devon Ryan81k
Please log in to add an answer.


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