Cigar string and read different length bug fixed.

lbcb-sci · Nov 4, 2019 · 3bd398d · 3bd398d
1 parent 220d611
commit 3bd398d
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/alignment/alignment.cc b/src/alignment/alignment.cc
@@ -297,7 +297,7 @@ int CheckAlignmentSane(std::vector<unsigned char> &alignment, const SingleSequen
           // If there are insertions following deletions (or other way around), something is wrong again.
           if ((last_move == EDLIB_I && alignment_char == EDLIB_D) || (last_move == EDLIB_D && alignment_char == EDLIB_I)) {
             LogSystem::GetInstance().Log(VERBOSE_LEVEL_ALL_DEBUG, true, FormatString("CheckAlignmentSane returned false! return 2. Consecutive I and D operations! last_move = %c, alignment_char = %c. num_same_moves = %ld, qname: '%s', read_length: %ld, ref_length: %ld\n", (last_move == EDLIB_I) ? 'I' : 'D', (alignment_char == EDLIB_I) ? 'I' : 'D', num_same_moves, read->get_header(), read_length, ref_length), "CheckAlignmentSane");
-            return 0;
+            return read_length != read->get_sequence_length();
           }
         }
         if (i < alignment.size()) {

diff --git a/src/graphmap/process_read.cc b/src/graphmap/process_read.cc
@@ -1136,7 +1136,7 @@ int GraphMap::RNAGenerateAlignments_(int order_number, MappingData *mapping_data
 	  cigarExons2.push_back(cigExon);
   }
 
-  double difference = abs(score22-score);
+  double difference = fabs(score22-score);
   betterCigarExons = cigarExons1;
 
   if(aln_result2->rv == is::AlignmentReturnValue::OK) {