Hello,
I'm using bowtie2 to align some whole-exome sequencing reads and have a basic question about the output I get.
The summary text includes the line X% aligned concordantly >1 times. In my files, X is generally pretty high (40% or so).
However, I did not specify -k or -a when I ran the program, which as far as I understand means the program should be reporting 1 alignment per read, NOT reporting "all valid alignments." How am I then getting so many reads that align more than once?
Am I misunderstanding what is meant by "aligned concordantly >1 times" or by what is meant by "report 1 valid alignment"?
Thank you!
•
written
by
am3
1 answer
From the bowtie2 manual:
Default mode: search for multiple alignments, report the best one
and
When it finds a valid alignment, it generally will continue to look for alignments that are nearly as good or better. It will eventually stop looking, either because it exceeded a limit placed on search effort (see -D and -R) or because it already knows all it needs to know to report an alignment.
So the summary text is probably just be reporting the runtime statistics as described above, not what's (reported) in your output files. You could check your SAM file to see if that's true.
Log in to answer this question.
-
Bowtie2 and Hisat2 read mapping comparisons
written by Allan •
-
unusual alignment rates for IgG and H3K27me3 in cuttag
written by david •
-
Cannot find summary message after Bowtie2 mapping
written by BeeWork •
-
bowtie2 results with NGS data
written by adnan.sabar7865 •
-
All sequenced reads in piared end sequencing don't have mate pairs?
written by piyushjo
-
bowtie2 summary looks very weird. Neither concordant nor discordant.
written by zhangzhenyusz •
-
Hisat2 versus Bowtie2
written by nalpas.nicolas •
-
Bowtie2 Alignment summary
written by bassanio •
-
unexpected bowtie2 unpaired alignment behavior
written by mark.rose •
-
Unaligned Reads From Bowtie2
written by deepthithomaskannan
