Question: How to identify bowtie2 unique mapping using sam flag value?
0
gravatar for jinxinhao1988
23 months ago by
China/Beijing/Celloud
jinxinhao198820 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.0k views
ADD COMMENTlink modified 3 months ago by Biostar ♦♦ 20 • written 23 months ago by jinxinhao198820

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

ADD REPLYlink written 23 months ago by Devon Ryan74k
0
gravatar for d-cameron
13 months ago by
d-cameron1.0k
Australia
d-cameron1.0k 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 13 months ago by d-cameron1.0k
0
gravatar for khalid.belkhir
13 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 13 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 13 months ago by Devon Ryan74k
Please log in to add an answer.

Help
Access

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