Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
bowtie2 local alignment frymor Bioinformatics 6 07-06-2015 10:39 AM
Techniques for estimating local alignment score distribution parameters cb1025 Bioinformatics 2 03-19-2015 06:55 AM
Local drop in Phred score around SNPs in mpileup -- are these false positives? r1cgw3 Bioinformatics 1 09-15-2013 02:53 PM
GATK recalibration score and local realignment: what is the first step? m_elena_bioinfo Bioinformatics 0 01-24-2011 03:59 AM
gap alignment and local alignment? mingkunli Illumina/Solexa 3 02-19-2009 12:13 PM

Thread Tools
Old 09-22-2015, 01:38 AM   #1
Junior Member
Location: United States

Join Date: Sep 2015
Posts: 1
Default seqan local alignment returns incorrect score

Hi, I am using SeqAn 1.3 installed using the debian package seqan-apps on Ubuntu 14.04. I am getting what appears to be an incorrect score when performing a Waterman-Eggert alignment.
Here is my code:
#include <seqan/align.h>

using namespace std;
using namespace seqan;
int main(void) {
Align< String<Dna> > ali2;
appendValue(rows(ali2), "CGATAACATTGAAAGA");
appendValue(rows(ali2), "CGATAAACTCAGTACA");
LocalAlignmentFinder<> finder(ali2);
Score<int> scoring(2,-3,-5,-2);
while (localAlignment(ali2, finder, scoring,
8, WatermanEggert())) {
cout << "Score = " << getScore(finder) << endl;
cout << ali2;
cout << "pos1: ";
cout << clippedBeginPosition(row(ali2,0)) << endl;
cout << "pos2: ";
cout << clippedBeginPosition(row(ali2,1)) << endl;

Here is the output that this program produces:
Score = 0
0 .

pos1: 0
pos2: 0

It shows a score of zero, but this should have a score of 12. In fact, the code is supposed to show only those alignments that have scores greater than 8. So internally, the program recognizes that this score is greater than 8, but it reports that the score is zero.

Can you please shed some light on why this is happening?

RynT is offline   Reply With Quote

seqan, waterman-eggert

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 10:54 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO