Added missing empty check.

dellytools · Jun 3, 2014 · d18123a · d18123a
1 parent 2b23372
commit d18123a
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 9 deletions.
diff --git a/src/delly.cpp b/src/delly.cpp
@@ -1219,24 +1219,28 @@ vcfOutput(TConfig const& c, std::vector<TStructuralVariantRecord> const& svs, TJ
 
       // Always take the minimum to be conservative (and to flag unclear samples as LowQual)
       if ((countMapItLeft!=normalCountMap.end()) && (countMapItRight!=normalCountMap.end())) {
-	if (countMapItLeft->second[0].size() <= countMapItRight->second[0].size()) countMapIt=countMapItLeft;
-	else countMapIt=countMapItRight;
+	if ((!countMapItLeft->second.empty()) && (!countMapItRight->second.empty())) {
+	  if (countMapItLeft->second[0].size() <= countMapItRight->second[0].size()) countMapIt=countMapItLeft;
+	  else countMapIt=countMapItRight;
+	}
       }
       if ((abCountMapItLeft!=abnormalCountMap.end()) && (abCountMapItRight!=abnormalCountMap.end())) {
-	if ((svIter->chr != svIter->chr2) || (abCountMapItLeft->second[0].size() > abCountMapItRight->second[0].size())) abCountMapIt=abCountMapItRight;
-	else abCountMapIt=abCountMapItLeft;
+	if ((!abCountMapItLeft->second.empty()) && (!abCountMapItRight->second.empty())) {
+	  if ((svIter->chr != svIter->chr2) || (abCountMapItLeft->second[0].size() > abCountMapItRight->second[0].size())) abCountMapIt=abCountMapItRight;
+	  else abCountMapIt=abCountMapItLeft;
+	}
       }
 
-      TMapqVector mapqRef;
-      TMapqVector mapqAlt;
+      TMapqVector mapqRef = TMapqVector();
+      TMapqVector mapqAlt = TMapqVector();
       if ((svIter->precise) || (c.junction)) {
 	// Genotyping for precise events uses junction read qualities
 	mapqRef = jctCountMapIt->second.first;
 	mapqAlt = jctCountMapIt->second.second;
       } else {
 	// Genotyping for imprecise events uses paired-end qualities
-	mapqRef = countMapIt->second[0];
-	mapqAlt = abCountMapIt->second[0];
+	if ((countMapIt!=normalCountMap.end()) && (!countMapIt->second.empty())) mapqRef = countMapIt->second[0];
+	if ((abCountMapIt!=abnormalCountMap.end()) && (!abCountMapIt->second.empty())) mapqAlt = abCountMapIt->second[0];
       }
 
       // Compute genotype likelihoods

diff --git a/src/version.h b/src/version.h
@@ -27,7 +27,7 @@ Contact: Tobias Rausch (rausch@embl.de)
 namespace torali
 {
 
-  std::string dellyVersionNumber = "0.5.4";
+  std::string dellyVersionNumber = "0.5.5";
 
   inline 
     void printTitle(std::string const& title)