diff --git a/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/Stats.scala b/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/Stats.scala
index 1023e280237ccb699de7f77f99634b9e947b9e82..d7f0961fc4dcaf2dd03f7773a9826908e9d00794 100644
--- a/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/Stats.scala
+++ b/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/Stats.scala
@@ -17,11 +17,10 @@ package nl.lumc.sasc.biopet.tools.vcfstats
 import java.io.{File, FileOutputStream, IOException, PrintWriter}
 
 import nl.lumc.sasc.biopet.tools.vcfstats.Stats.plotHeatmap
-import nl.lumc.sasc.biopet.tools.vcfstats.VcfStats.{getClass, logger, sampleDistributions}
-
-import scala.collection.mutable
+import nl.lumc.sasc.biopet.tools.vcfstats.VcfStats.logger
 import nl.lumc.sasc.biopet.utils.{ConfigUtils, sortAnyAny}
 
+import scala.collection.mutable
 import scala.sys.process.{Process, ProcessLogger}
 
 /**
@@ -161,26 +160,32 @@ case class Stats(generalStats: mutable.Map[String, mutable.Map[Any, Int]] = muta
 
   def writeAllOutput(outputDir: File,
                      samples: List[String],
-                     genotypeFields: List[String] = Nil,
-                     infoFields: List[String] = Nil,
-                     sampleDistributions: List[String] = Nil): Unit = {
+                     genotypeFields: List[String],
+                     infoFields: List[String],
+                     sampleDistributions: List[String],
+                     contig: Option[String]): Unit = {
     outputDir.mkdirs()
     this.writeToFile(new File(outputDir, "stats.json"),
                      samples,
                      genotypeFields,
                      infoFields,
-                     sampleDistributions)
+                     sampleDistributions,
+                     contig)
     writeOverlap(outputDir, samples)
   }
 
   def writeToFile(outputFile: File,
                   samples: List[String],
-                  genotypeFields: List[String] = Nil,
-                  infoFields: List[String] = Nil,
-                  sampleDistributions: List[String] = Nil): Unit = {
+                  genotypeFields: List[String],
+                  infoFields: List[String],
+                  sampleDistributions: List[String],
+                  contig: Option[String]): Unit = {
     val allWriter = new PrintWriter(outputFile)
-    val json = ConfigUtils.mapToJson(
-      this.getStatsAsMap(samples, genotypeFields, infoFields, sampleDistributions))
+    val map = this.getStatsAsMap(samples, genotypeFields, infoFields, sampleDistributions)
+    val json = contig match {
+      case Some(c) => ConfigUtils.mapToJson(Map("contigs" -> Map(c -> map)))
+      case _ => ConfigUtils.mapToJson(Map("total" -> map))
+    }
     allWriter.println(json.nospaces)
     allWriter.close()
   }
diff --git a/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/VcfStats.scala b/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/VcfStats.scala
index d2cafb767514562a2809d06889522a577fa2821d..cd49082ae5110fc70e6700d7a65de73cbd90f7a5 100644
--- a/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/VcfStats.scala
+++ b/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/vcfstats/VcfStats.scala
@@ -89,7 +89,8 @@ object VcfStats extends ToolCommand {
                          samples,
                          adGenotypeTags,
                          adInfoTags,
-                         sampleDistributions))
+                         sampleDistributions,
+                         None))
 
     regionStats
       .flatMap(_._2)
@@ -100,7 +101,8 @@ object VcfStats extends ToolCommand {
                            samples,
                            adGenotypeTags,
                            adInfoTags,
-                           sampleDistributions)
+                           sampleDistributions,
+                           Some(k))
       }
     regionStats.unpersist()
 
@@ -159,7 +161,8 @@ object VcfStats extends ToolCommand {
             samples,
             adGenotypeTags,
             adInfoTags,
-            sampleDistributions)
+            sampleDistributions,
+            Some(bedRecord.chr))
           None
         }
       } else Future.successful(Some(bedRecord.chr -> stats))