7 * Created by Sarah Westcott on 5/7/09.
8 * Copyright 2009 Schloss Lab UMASS Amherst. All rights reserved.
15 /**************************************************************************************************/
17 class oneGapDist : public Dist {
22 oneGapDist(const oneGapDist& ddb) {}
24 void calcDist(Sequence A, Sequence B){
32 string seqA = A.getAligned();
33 string seqB = B.getAligned();
34 int alignLength = seqA.length();
36 for(int i=0;i<alignLength;i++){
37 if((seqA[i] != '.' || seqB[i] != '.')){
43 for(int i=start;i<alignLength;i++){
44 if((seqA[i] == '-' && seqB[i] == '-') || (seqA[i] == '.' && seqB[i] == '-') || (seqA[i] == '-' && seqB[i] == '.')){ ; }
45 else if(seqA[i] == '.' && seqB[i] == '.'){
48 else if(seqB[i] != '-' && (seqA[i] == '-' || seqA[i] == '.')){
56 else if(seqA[i] != '-' && (seqB[i] == '-' || seqB[i] == '.')){
64 else if(seqA[i] != '-' && seqB[i] != '-'){
65 if(seqA[i] != seqB[i]){
79 if(minLength == 0) { dist = 1.0000; }
80 else { dist = (double)difference / minLength; }
85 /**************************************************************************************************/