Question: Understanding Bowtie output
0
gravatar for am3
8 weeks ago by
am370
am370 wrote:

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!

bowtie alignment • 147 views
ADD COMMENTlink modified 8 weeks ago by manuel.belmadani730 • written 8 weeks ago by am370
1
gravatar for manuel.belmadani
8 weeks ago by
Canada
manuel.belmadani730 wrote:

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.

ADD COMMENTlink modified 8 weeks ago • written 8 weeks ago by manuel.belmadani730
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: 1850 users visited in the last hour