Bowtie2 adds insertions instead of allowing extention
0
0
Entering edit mode
4.5 years ago
chrisgr ▴ 20

Say, I'm aligning, with Bowtie2 default settings, this short read:

>read
AGTGATCAGTGGTACACCATGTACGAATTCTCGCCTGCGAGACGA

to this short reference:

>oligo40
AGTGATCAGTGGTACACCATGTACGAATTCTCGCCTGCGA

The read is a bit longer than the reference. The problem: Bowtie aligns it resulting in this CIGAR string: 36M5I4M Which looks like:

reference:AGTGATCAGTGGTACACCATGTACGAATTCTCGCCT     GCGA
read:     AGTGATCAGTGGTACACCATGTACGAATTCTCGCCTGCGAGACGA

Bowtie2 'made up' an insertion only for the purpose to fit the whole read on the reference. But I want to allow for the read to extend the reference (perhaps assign a penalty to it, but not throw in a insertion just like that). How do I do that? Can't seem to find it anywhere.

Many thanks in advance!

bowtie2 short reads aligning insertions • 768 views
ADD COMMENT
1
Entering edit mode

Doesn't seem bowtie2 is the correct tool for this task. What exactly are you trying to accomplish? If you are trying to extend the reference with reads, you can try tadpole.sh from the BBTools / BBMap package suite.

ADD REPLY
1
Entering edit mode

Don't know if you are able to use an alternate aligner but bbmap.sh from BBMap suite aligns the data soft-clipping the not matching part like this

read    0       oligo40 1       30      40=5S   *       0       0       AGTGATCAGTGGTACACCATGTACGAATTCTCGCCTGCGAGACGA   *       NM:i:0  AM:i:30
ADD REPLY

Login before adding your answer.

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