diff --git a/libsequence/src/libsequence/Sequence/PolyTable.hpp b/libsequence/src/libsequence/Sequence/PolyTable.hpp
index 59ee93f..38b03d3 100644
--- a/libsequence/src/libsequence/Sequence/PolyTable.hpp
+++ b/libsequence/src/libsequence/Sequence/PolyTable.hpp
@@ -53,6 +53,7 @@ long with libsequence.  If not, see <http://www.gnu.org/licenses/>.
 #include <type_traits>
 #include <functional>
 #include <memory>
+#include <stdexcept>
 #include <Sequence/polySiteVector.hpp>
 /*! \example PolyTableIterators.cc */
 namespace Sequence
diff --git a/libsequence/src/libsequence/Sequence/Translate.hpp b/libsequence/src/libsequence/Sequence/Translate.hpp
index 6efba67..d08d660 100644
--- a/libsequence/src/libsequence/Sequence/Translate.hpp
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/Sequence/VariantMatrixCapsule.hpp b/libsequence/src/libsequence/Sequence/VariantMatrixCapsule.hpp
index 7d1e6b6..9ebb783 100644
--- a/libsequence/src/libsequence/Sequence/VariantMatrixCapsule.hpp
+++ b/libsequence/src/libsequence/Sequence/VariantMatrixCapsule.hpp
@@ -4,6 +4,7 @@
 #include <cstdint>
 #include <memory>
 #include <vector>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/libsequence/src/libsequence/Sequence/bits/AlignStream.tcc b/libsequence/src/libsequence/Sequence/bits/AlignStream.tcc
index 43db3b8..a9663eb 100644
--- a/libsequence/src/libsequence/Sequence/bits/AlignStream.tcc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/Sequence/bits/Alignment.tcc b/libsequence/src/libsequence/Sequence/bits/Alignment.tcc
index 1ea92e5..7bc6f6a 100644
--- a/libsequence/src/libsequence/Sequence/bits/Alignment.tcc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/Sequence/bits/Clustalw.tcc b/libsequence/src/libsequence/Sequence/bits/Clustalw.tcc
index d56967a..5ca0c95 100644
--- a/libsequence/src/libsequence/Sequence/bits/Clustalw.tcc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/Sequence/bits/phylipData.tcc b/libsequence/src/libsequence/Sequence/bits/phylipData.tcc
index bdb00e7..ea74eda 100644
--- a/libsequence/src/libsequence/Sequence/bits/phylipData.tcc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/src/PolySites.cc b/libsequence/src/libsequence/src/PolySites.cc
index 5682c70..7be9a61 100644
--- a/libsequence/src/libsequence/src/PolySites.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/src/RedundancyCom95.cc b/libsequence/src/libsequence/src/RedundancyCom95.cc
index 5327294..eefcd24 100644
--- a/libsequence/src/libsequence/src/RedundancyCom95.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/src/Seq/fastq.cc b/libsequence/src/libsequence/src/Seq/fastq.cc
index 84cdb8f..7b067c2 100644
--- a/libsequence/src/libsequence/src/Seq/fastq.cc
+++ b/libsequence/src/libsequence/src/Seq/fastq.cc
@@ -1,6 +1,7 @@
 #include <Sequence/fastq.hpp>
 #include <functional>
 #include <iostream>
+#include <stdexcept>
 
 namespace Sequence
 {
diff --git a/libsequence/src/libsequence/src/Specializations.cc b/libsequence/src/libsequence/src/Specializations.cc
index 5a19b93..23d1edc 100644
--- a/libsequence/src/libsequence/src/Specializations.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/src/shortestPath.cc b/libsequence/src/libsequence/src/shortestPath.cc
index 563ce9f..0a9e660 100644
--- a/libsequence/src/libsequence/src/shortestPath.cc
+++ b/libsequence/src/libsequence/src/shortestPath.cc
@@ -21,6 +21,7 @@
 
 */
 
+#include <stdexcept>
 #include <Sequence/shortestPath.hpp>
 #include <Sequence/PathwayHelper.hpp>
 #include <Sequence/Translate.hpp>
diff --git a/libsequence/src/libsequence/src/summstats_deprecated/Recombination.cc b/libsequence/src/libsequence/src/summstats_deprecated/Recombination.cc
index 02a3319..649a607 100644
--- a/libsequence/src/libsequence/src/summstats_deprecated/Recombination.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/src/summstats_deprecated/nSL.cc b/libsequence/src/libsequence/src/summstats_deprecated/nSL.cc
index 54a346f..7e35d43 100644
--- a/libsequence/src/libsequence/src/summstats_deprecated/nSL.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/src/variant_matrix/nonowningcapsules.cc b/libsequence/src/libsequence/src/variant_matrix/nonowningcapsules.cc
index b9f3e85..726e626 100644
--- a/libsequence/src/libsequence/src/variant_matrix/nonowningcapsules.cc
+++ b/libsequence/src/libsequence/src/variant_matrix/nonowningcapsules.cc
@@ -3,6 +3,7 @@
 #include <algorithm>
 #include <iterator>
 #include <exception>
+#include <stdexcept>
 
 namespace
 {
diff --git a/libsequence/src/libsequence/test/AlignStreamTest.cc b/libsequence/src/libsequence/test/AlignStreamTest.cc
index 73123cf..5310de9 100644
--- a/libsequence/src/libsequence/test/AlignStreamTest.cc
+++ b/libsequence/src/libsequence/test/AlignStreamTest.cc
@@ -7,6 +7,7 @@
 #include <fstream>
 #include <iostream>
 #include <unistd.h>
+#include <stdexcept>
 
 struct infilenames
 {
diff --git a/libsequence/src/libsequence/test/ComparisonsTest.cc b/libsequence/src/libsequence/test/ComparisonsTest.cc
index b20daf3..c440e3f 100644
--- a/libsequence/src/libsequence/test/ComparisonsTest.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/test/FastaIO.cc b/libsequence/src/libsequence/test/FastaIO.cc
index 26a3ecf..a215b81 100644
--- a/libsequence/src/libsequence/test/FastaIO.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/test/PolyTableBadBehavior.cc b/libsequence/src/libsequence/test/PolyTableBadBehavior.cc
index 837a7a7..60e65e8 100644
--- a/libsequence/src/libsequence/test/PolyTableBadBehavior.cc
+++ b/libsequence/src/libsequence/test/PolyTableBadBehavior.cc
@@ -18,6 +18,7 @@
 #include <cctype>
 #include <iostream>
 #include <functional>
+#include <stdexcept>
 
 BOOST_AUTO_TEST_SUITE(PolyTableBadBehaviorTest)
 
diff --git a/libsequence/src/libsequence/test/testAlleleCountMatrix.cc b/libsequence/src/libsequence/test/testAlleleCountMatrix.cc
index 19941bd..197f340 100644
--- a/libsequence/src/libsequence/test/testAlleleCountMatrix.cc
+++ b/libsequence/src/libsequence/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/libsequence/src/libsequence/test/testClassicSummstats.cc b/libsequence/src/libsequence/test/testClassicSummstats.cc
index 5680bc2..a113e22 100644
--- a/libsequence/src/libsequence/test/testClassicSummstats.cc
+++ b/libsequence/src/libsequence/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>
diff --git a/libsequence/src/polytable.cc b/libsequence/src/polytable.cc
index 76c6e90..4e9ba12 100644
--- a/libsequence/src/polytable.cc
+++ b/libsequence/src/polytable.cc
@@ -18,6 +18,7 @@
 //
 
 #include <iostream>
+#include <stdexcept>
 #include <pybind11/pybind11.h>
 #include <pybind11/functional.h>
 #include <pybind11/stl.h>
diff --git a/libsequence/src/variant_matrix.cc b/libsequence/src/variant_matrix.cc
index b5b39a5..629dbe1 100644
--- a/libsequence/src/variant_matrix.cc
+++ b/libsequence/src/variant_matrix.cc
@@ -1,4 +1,5 @@
 #include <iostream>
+#include <stdexcept>
 #include <pybind11/pybind11.h>
 #include <pybind11/functional.h>
 #include <pybind11/numpy.h>
