Question: How to identify bowtie2 unique mapping using sam flag value?
0
gravatar for jinxinhao1988
19 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 • 728 views
ADD COMMENTlink modified 9 months ago by khalid.belkhir20 • written 19 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 19 months ago by Devon Ryan70k
0
gravatar for d-cameron
9 months ago by
d-cameron770
Australia
d-cameron770 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 9 months ago by d-cameron770
0
gravatar for khalid.belkhir
9 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 9 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 9 months ago by Devon Ryan70k
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: 1447 users visited in the last hour