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

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


to this short reference:


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


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

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 from the BBTools / BBMap package suite.

ADD REPLYlink modified 4 months ago • written 4 months ago by h.mon31k

Don't know if you are able to use an alternate aligner but 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.


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