Bowtie: Reporting one of multiple valid mapping locations
1
0
Entering edit mode
9.1 years ago
Sander ▴ 20

Dear Biostars,

I'm slightly confused by the way Bowtie reports reads that align to multiple locations. If I understand correctly from the manual, by default only 1 alignment is returned (-k 1). But this may not necessarily be the best alignment, because then --best must be specified as well. So is this a random or just a 'good' alignment? And what is the difference with -M 1? According to the manual this gives 1 random alignment from the best alignments, but is this not the same as -k 1?

Kinds regards,
Sander

bowtie mapping • 2.4k views
ADD COMMENT
0
Entering edit mode
9.1 years ago

-k is not related to 'best' alignment. -k is related to multimapped reads. A read may map at various places with v.high quality/best mapping, but only one position is reported by default.

ADD COMMENT
0
Entering edit mode

So with -k 1, in the case that a read maps to various places, which position is reported?

ADD REPLY

Login before adding your answer.

Traffic: 2036 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