Question: Bowtie2 adds insertions instead of allowing extention
0
gravatar for christiangriffioen
4 months ago by
christiangriffioen20 wrote:

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!

ADD COMMENTlink written 4 months ago by christiangriffioen20
1

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 REPLYlink modified 4 months ago • written 4 months ago by h.mon31k
1

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 REPLYlink modified 4 months ago • written 4 months ago by genomax91k
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: 1377 users visited in the last hour