This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Alignment with deletion

Hello,

I have a DNA sequence with some deletions. I have a reference sequence, with deletions as well.

I want to align these reads pairwise. However, when I'm using Align packages (or others), it considers the deletions as mismatches: enter image description here

Is there a method to count deletions as matches if performed in both sequences?

Thank you!

match deletion alignment

I want to align these reads pairwise.

How will the information be carried over if you are doing pairwise alignments? e.g. how would the aligner know to put a gap in the sequences if it is simply aligning the two sequences you show in the example above.

1 answer

You can use EMBOSS Needle for DNA:

EMBOSS Needle

Thank you. Do you familiar with python package that I can integrate into a pipeline?

Log in to answer this question.