bowtie1 parameter -k
1
0
Entering edit mode
8.7 years ago
biolab ★ 1.4k

Hi, everyone, I have a simple question regarding to the -k parameter of bowtie1. The manual states "report up to <int> good alignments per read (default: 1)". If I set -k 1, for a multiple read, is it randomly assigned to a position?
Thank you very much!

bowtie1 • 1.8k views
ADD COMMENT
1
Entering edit mode
8.7 years ago
mastal511 ★ 2.1k

Yes. It just reports the first alignment.

Quote from the bowtie1 manual:

http://bowtie-bio.sourceforge.net/manual.shtml#example-2--k-3

"Leaving the reporting options at their defaults causes bowtie to report the first valid alignment it encounters. Because --best was not specified, we are not guaranteed that bowtie will report the best alignment, and in this case it does not (the 1-mismatch alignment from the previous example would have been better). The default reporting mode is equivalent to -k 1."

ADD COMMENT
0
Entering edit mode

Thank you, mastal511. It's very helpful!

ADD REPLY

Login before adding your answer.

Traffic: 1638 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6