Question: How to identify bowtie2 unique mapping using sam flag value?
gravatar for jinxinhao1988
21 months ago by
China/Beijing/Beijing Institute of Genomics,CAS
jinxinhao19880 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 • 921 views
ADD COMMENTlink modified 7 weeks ago by Biostar ♦♦ 20 • written 21 months ago by jinxinhao19880

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

ADD REPLYlink written 21 months ago by Devon Ryan71k
gravatar for d-cameron
11 months ago by
d-cameron850 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 11 months ago by d-cameron850
gravatar for khalid.belkhir
11 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 11 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 11 months ago by Devon Ryan71k
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: 771 users visited in the last hour