From e5ceb0bcef00f5a2637e5cf9b36132aff3f83621 Mon Sep 17 00:00:00 2001 From: tobiasrausch Date: Wed, 15 Jul 2015 12:48:36 +0200 Subject: [PATCH] fixed translocation genotyping bug --- src/spanning.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/spanning.h b/src/spanning.h index 2745c3b..8e974c7 100644 --- a/src/spanning.h +++ b/src/spanning.h @@ -248,7 +248,7 @@ namespace torali { if (std::abs(midPoint - itSV->svStart) < std::abs(itSV->svEnd - midPoint)) { if ((itSV->chr==rec->core.tid) && (itSV->svStart>=sPos) && (itSV->svStart<=ePos)) leftIt->second.first.push_back(pairQuality); } else { - if ((itSV->chr==rec->core.tid) && (itSV->svEnd>=sPos) && (itSV->svEnd<=ePos)) rightIt->second.first.push_back(pairQuality); + if ((itSV->chr2==rec->core.tid) && (itSV->svEnd>=sPos) && (itSV->svEnd<=ePos)) rightIt->second.first.push_back(pairQuality); } } else if ((getStrandIndependentOrientation(rec->core) != libIt->second.defaultOrient) || (outerISize < libIt->second.minNormalISize) || (outerISize > libIt->second.maxNormalISize) || (rec->core.tid!=rec->core.mtid)) { // Missing spanning coverage