diff --git a/Sequence/Translate.hpp b/Sequence/Translate.hpp
index 6efba67..d08d660 100644
--- a/Sequence/Translate.hpp
+++ b/Sequence/Translate.hpp
@@ -24,6 +24,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 #ifndef __TRANSLATE_HPP__
 #define __TRANSLATE_HPP__
 #include <string>
+#include <stdexcept>
 #include <Sequence/SeqEnums.hpp>
 /*! \file Translate.hpp
   @brief declares Sequence::Translate,a function to translate CDS sequences into peptide sequences
diff --git a/Sequence/VariantMatrixCapsule.hpp b/Sequence/VariantMatrixCapsule.hpp
index 7d1e6b6..9ebb783 100644
--- a/Sequence/VariantMatrixCapsule.hpp
+++ b/Sequence/VariantMatrixCapsule.hpp
@@ -4,6 +4,7 @@
 #include <cstdint>
 #include <memory>
 #include <vector>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/Sequence/bits/AlignStream.tcc b/Sequence/bits/AlignStream.tcc
index 43db3b8..a9663eb 100644
--- a/Sequence/bits/AlignStream.tcc
+++ b/Sequence/bits/AlignStream.tcc
@@ -28,6 +28,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 */
 
 #include <Sequence/AlignStream.hpp>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/Sequence/bits/Alignment.tcc b/Sequence/bits/Alignment.tcc
index 1ea92e5..7bc6f6a 100644
--- a/Sequence/bits/Alignment.tcc
+++ b/Sequence/bits/Alignment.tcc
@@ -37,6 +37,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 #include <iterator>
 #include <algorithm>
 #include <iostream>
+#include <stdexcept>
 namespace Sequence
 {
     namespace Alignment
diff --git a/Sequence/bits/Clustalw.tcc b/Sequence/bits/Clustalw.tcc
index d56967a..5ca0c95 100644
--- a/Sequence/bits/Clustalw.tcc
+++ b/Sequence/bits/Clustalw.tcc
@@ -31,6 +31,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 #include <iterator>
 #include <algorithm>
 #include <cassert>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/Sequence/bits/PolyTableSlice.tcc b/Sequence/bits/PolyTableSlice.tcc
index 9017c16..5e2cc70 100644
--- a/Sequence/bits/PolyTableSlice.tcc
+++ b/Sequence/bits/PolyTableSlice.tcc
@@ -27,6 +27,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 #include <iostream>
 #include <stdexcept>
 #include <cmath>
+#include <stdexcept>
 namespace Sequence
 {
   template<typename T>
diff --git a/Sequence/bits/phylipData.tcc b/Sequence/bits/phylipData.tcc
index bdb00e7..ea74eda 100644
--- a/Sequence/bits/phylipData.tcc
+++ b/Sequence/bits/phylipData.tcc
@@ -29,6 +29,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 #include <utility>
 #include <algorithm>
 #include <cctype>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/src/PolySites.cc b/src/PolySites.cc
index 5682c70..7be9a61 100644
--- a/src/PolySites.cc
+++ b/src/PolySites.cc
@@ -27,6 +27,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 #include <iterator>
 #include <algorithm>
 #include <iostream>
+#include <stdexcept>
 
 /*! \class Sequence::PolySites Sequence/PolySites.hpp
   \ingroup polytables
diff --git a/src/RedundancyCom95.cc b/src/RedundancyCom95.cc
index 5327294..eefcd24 100644
--- a/src/RedundancyCom95.cc
+++ b/src/RedundancyCom95.cc
@@ -27,6 +27,7 @@
 #include <cstring>
 #include <algorithm>
 #include <functional>
+#include <stdexcept>
 #include <Sequence/SeqEnums.hpp>
 #include <Sequence/Translate.hpp>
 #include <Sequence/RedundancyCom95.hpp>
diff --git a/src/Seq/fastq.cc b/src/Seq/fastq.cc
index 84cdb8f..7b067c2 100644
--- a/src/Seq/fastq.cc
+++ b/src/Seq/fastq.cc
@@ -1,6 +1,7 @@
 #include <Sequence/fastq.hpp>
 #include <functional>
 #include <iostream>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/src/Specializations.cc b/src/Specializations.cc
index 5a19b93..23d1edc 100644
--- a/src/Specializations.cc
+++ b/src/Specializations.cc
@@ -23,6 +23,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 
 //Code for the -*- C++ -*- Template Specializations for libsequence
 #include <Sequence/Alignment.hpp>
+#include <stdexcept>
 
 /*! \file Specializations.cc
   @brief Definitions of template specializations for library functions
diff --git a/src/shortestPath.cc b/src/shortestPath.cc
index 563ce9f..8e60efc 100644
--- a/src/shortestPath.cc
+++ b/src/shortestPath.cc
@@ -29,6 +29,7 @@
 #include <Sequence/SeqAlphabets.hpp>
 #include <Sequence/Comparisons.hpp>
 #include <algorithm>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/src/summstats_deprecated/Recombination.cc b/src/summstats_deprecated/Recombination.cc
index 02a3319..649a607 100644
--- a/src/summstats_deprecated/Recombination.cc
+++ b/src/summstats_deprecated/Recombination.cc
@@ -24,6 +24,7 @@
 #include <algorithm>
 #include <limits>
 #include <cctype>
+#include <stdexcept>
 #include <Sequence/PolyTable.hpp>
 #include <Sequence/SimData.hpp>
 #include <Sequence/Recombination.hpp>
diff --git a/src/summstats_deprecated/nSL.cc b/src/summstats_deprecated/nSL.cc
index 54a346f..7e35d43 100644
--- a/src/summstats_deprecated/nSL.cc
+++ b/src/summstats_deprecated/nSL.cc
@@ -6,6 +6,7 @@
 #include <cmath>
 #include <limits>
 #include <unordered_map>
+#include <stdexcept>
 // For parallelizing nSL:
 #include <functional>
 #include <iostream>
diff --git a/src/variant_matrix/nonowningcapsules.cc b/src/variant_matrix/nonowningcapsules.cc
index b9f3e85..726e626 100644
--- a/src/variant_matrix/nonowningcapsules.cc
+++ b/src/variant_matrix/nonowningcapsules.cc
@@ -3,6 +3,7 @@
 #include <algorithm>
 #include <iterator>
 #include <exception>
+#include <stdexcept>
 
 namespace
 {
diff --git a/test/AlignStreamTest.cc b/test/AlignStreamTest.cc
index 73123cf..5310de9 100644
--- a/test/AlignStreamTest.cc
+++ b/test/AlignStreamTest.cc
@@ -7,6 +7,7 @@
 #include <fstream>
 #include <iostream>
 #include <unistd.h>
+#include <stdexcept>
 
 struct infilenames
 {
diff --git a/test/ComparisonsTest.cc b/test/ComparisonsTest.cc
index b20daf3..c440e3f 100644
--- a/test/ComparisonsTest.cc
+++ b/test/ComparisonsTest.cc
@@ -5,6 +5,7 @@
 #include <boost/test/unit_test.hpp>
 #include <algorithm>
 #include <iterator>
+#include <stdexcept>
 BOOST_AUTO_TEST_SUITE(ComparisonsTest)
 
 BOOST_AUTO_TEST_CASE( notagap ) //Silly!
diff --git a/test/FastaIO.cc b/test/FastaIO.cc
index 26a3ecf..a215b81 100644
--- a/test/FastaIO.cc
+++ b/test/FastaIO.cc
@@ -4,6 +4,7 @@
 #include <boost/test/unit_test.hpp>
 #include <unistd.h>
 #include <iterator>
+#include <stdexcept>
 
 struct fasta_io_fixture
 {
diff --git a/test/PolyTableBadBehavior.cc b/test/PolyTableBadBehavior.cc
index 837a7a7..60e65e8 100644
--- a/test/PolyTableBadBehavior.cc
+++ b/test/PolyTableBadBehavior.cc
@@ -18,6 +18,7 @@
 #include <cctype>
 #include <iostream>
 #include <functional>
+#include <stdexcept>
 
 BOOST_AUTO_TEST_SUITE(PolyTableBadBehaviorTest)
 
diff --git a/test/testAlleleCountMatrix.cc b/test/testAlleleCountMatrix.cc
index 19941bd..197f340 100644
--- a/test/testAlleleCountMatrix.cc
+++ b/test/testAlleleCountMatrix.cc
@@ -6,6 +6,7 @@
 #include <algorithm>
 #include <numeric> //for std::iota
 #include <iterator>
+#include <stdexcept>
 
 BOOST_FIXTURE_TEST_SUITE(test_allele_count_matrix, vmatrix_from_msprime)
 
diff --git a/test/testClassicSummstats.cc b/test/testClassicSummstats.cc
index 5680bc2..a113e22 100644
--- a/test/testClassicSummstats.cc
+++ b/test/testClassicSummstats.cc
@@ -7,6 +7,7 @@
 #include <string>
 #include <vector>
 #include <iostream>
+#include <stdexcept>
 #include <Sequence/VariantMatrix.hpp>
 #include <Sequence/VariantMatrixViews.hpp>
 #include <Sequence/summstats/classics.hpp>
