Dear all,
I would like to use BWA mem to output the best alignment for each read instead of all alignments, just like the "-k" argument below in Bowtie2, however, I haven't found the same function in BWA mem, which my version is 0.7.17-r1188, did I miss something or how can I tackle it?
-k <int> report up to <int> alns per read; MAPQ not meaningful
OR
-a/--all report all alignments; very slow, MAPQ not meaningful
•
written
by
zhangdengwei
1 answer
did I miss something
Reporting only the best alignment is the default in bwa mem (and bowtie2). By default, bwa reports some alternate alignments in the XA tag but you can ignore that. Note that for some reads the best alignment is chimeric so multiple alignment record will be reported for that read, linked by the SA tag. Note that these record collectively form a single alignment and are not considered to be alternate alignments.
Yes.
XA indicates an alternative alignment location. Read with XA indicate there could be other equally/almost as good alignment locations for this read.
SA indicates that this alignment is part of a chimeric alignment (aka split read alignment). Reads with SA imply the read spans a structural variant with a partial alignment on either side.
Log in to answer this question.
-
How does featureCounts handle multimapped reads from Bowtie2 -k 100 in default mode?
written by rachana •
-
problem to aligning FASTQ through BWA aligner
written by modarzi
-
Bowtie alignment reporting
written by VineethVenumadhavan •
-
Getting Bowtie2 or BWA to output more than one top alignment (not just a random one)
written by anm5579 •
-
Parallel for bwa mem - problem with -R argument for ID and SM
written by Korsocius
-
MAPQ values are 60 in bwa mem
written by John
-
bwa mem has more MAPQ values as 0 and 60
written by John
-
BWA-mem MAPQ Calculation
written by rfriedman22 •
-
Bwa mem is giving me alignements that are causing noise in subsequent analysis. How I can filter t…
written by shinken123
-
bowtie1 parameter -k
written by biolab
