From 4984ff265bb03f99203d771a4209d13be89d6be2 Mon Sep 17 00:00:00 2001
From: Tobias Rausch <rauschtobi@gmail.com>
Date: Fri, 6 Dec 2024 10:57:07 +0100
Subject: [PATCH] v1.3.2

---
 src/modvcf.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/modvcf.h b/src/modvcf.h
index c578df1..755d2e1 100644
--- a/src/modvcf.h
+++ b/src/modvcf.h
@@ -319,7 +319,7 @@ vcfParse(TConfig const& c, bam_hdr_t* hd, std::vector<TStructuralVariantRecord>&
 
 template<typename TConfig, typename TStructuralVariantRecord, typename TJunctionCountMap, typename TReadCountMap, typename TCountMap>
 inline void
-vcfOutput(TConfig const& c, std::vector<TStructuralVariantRecord>& svs, TJunctionCountMap const& jctCountMap, TReadCountMap const& readCountMap, TCountMap const& spanCountMap)
+vcfOutput(TConfig const& c, std::vector<TStructuralVariantRecord> const& svs, TJunctionCountMap const& jctCountMap, TReadCountMap const& readCountMap, TCountMap const& spanCountMap)
 {
   // BoLog class
   BoLog<double> bl;
@@ -416,7 +416,7 @@ vcfOutput(TConfig const& c, std::vector<TStructuralVariantRecord>& svs, TJunctio
     now = boost::posix_time::second_clock::local_time();
     std::cerr << '[' << boost::posix_time::to_simple_string(now) << "] " << "Genotyping" << std::endl;
     bcf1_t *rec = bcf_init();
-    for(typename TSVs::iterator svIter = svs.begin(); svIter!=svs.end(); ++svIter) {
+    for(typename TSVs::const_iterator svIter = svs.begin(); svIter!=svs.end(); ++svIter) {
       if ((svIter->srSupport == 0) && (svIter->peSupport == 0)) continue;
       // In discovery mode, skip SVs that have less than 2 reads support after genotyping
       if (!c.hasVcfFile) {