Manage Deletions - C++
1
0
Entering edit mode
9.1 years ago
MarioRio ▴ 40

Hi, my program analyze mutations in a BAM file after a comparison with the FASTA. The problem comes when an alignment has a base deletion, because the next bases move back in position clinging at the rest of the read. Is there a a function or something else in C++ bamtools library that can automatically extract the deletion information and manage it? If not, how can i resolve it?

thanks in advance

c bam • 1.6k views
ADD COMMENT
1
Entering edit mode
9.1 years ago

You could try to align the sequence of the read with the reference using a Smith-Waterman algorithm.

The second result of https://www.google.com/search?q=Smith+Waterman+filetype%3Acpp is a link to BandedSmithWaterman.cpp

"CBandedSmithWaterman - aligns reads using a banded Smith-Waterman algorithm"

https://code.google.com/p/mosaik-aligner/source/browse/src/CommonSource/PairwiseAlignment/BandedSmithWaterman.cpp

ADD COMMENT

Login before adding your answer.

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