diff --git a/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/bamstats/Histogram.scala b/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/bamstats/Histogram.scala
index ce0c7fce8aa8cad2f01b6a899df90a179e5459ec..56644eb9952588eef5684386b812108cbae27828 100644
--- a/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/bamstats/Histogram.scala
+++ b/biopet-tools/src/main/scala/nl/lumc/sasc/biopet/tools/bamstats/Histogram.scala
@@ -82,4 +82,11 @@ class Histogram[T](_counts: Map[T, Long] = Map[T, Long]())(implicit ord: Numeric
     } else Map()
   }
 
+  /** Write histogram to a tsv/count file */
+  def writeAggregateToTsv(file: File): Unit = {
+    val writer = new PrintWriter(file)
+    aggregateStats.foreach(x => writer.println(x._1 + "\t" + x._2))
+    writer.close()
+  }
+
 }