diff --git a/PyNutil/coordinate_extraction.py b/PyNutil/coordinate_extraction.py
index 3f29d1bac41278d3fe87f452b22879f39bf275c9..9e1c3fd0c13e2b20751ffee7171c19a4c644da46 100644
--- a/PyNutil/coordinate_extraction.py
+++ b/PyNutil/coordinate_extraction.py
@@ -1,6 +1,6 @@
 import numpy as np
 import pandas as pd
-from .read_and_write import load_visualign_json
+from .io.read_and_write import load_visualign_json
 from .counting_and_load import flat_to_dataframe
 from .visualign_deformations import triangulate, transform_vec
 from glob import glob
@@ -8,7 +8,7 @@ import cv2
 from skimage import measure
 import threading
 import re
-from .reconstruct_dzi import reconstruct_dzi
+from .io.reconstruct_dzi import reconstruct_dzi
 
 
 def number_sections(filenames, legacy=False):
diff --git a/PyNutil/io/__init__.py b/PyNutil/io/__init__.py
new file mode 100644
index 0000000000000000000000000000000000000000..b6e690fd59145ce8900fd9ab8d8a996ee7d33834
--- /dev/null
+++ b/PyNutil/io/__init__.py
@@ -0,0 +1 @@
+from . import *
diff --git a/PyNutil/atlas_loader.py b/PyNutil/io/atlas_loader.py
similarity index 100%
rename from PyNutil/atlas_loader.py
rename to PyNutil/io/atlas_loader.py
diff --git a/PyNutil/file_operations.py b/PyNutil/io/file_operations.py
similarity index 100%
rename from PyNutil/file_operations.py
rename to PyNutil/io/file_operations.py
diff --git a/PyNutil/propagation.py b/PyNutil/io/propagation.py
similarity index 100%
rename from PyNutil/propagation.py
rename to PyNutil/io/propagation.py
diff --git a/PyNutil/read_and_write.py b/PyNutil/io/read_and_write.py
similarity index 100%
rename from PyNutil/read_and_write.py
rename to PyNutil/io/read_and_write.py
diff --git a/PyNutil/reconstruct_dzi.py b/PyNutil/io/reconstruct_dzi.py
similarity index 100%
rename from PyNutil/reconstruct_dzi.py
rename to PyNutil/io/reconstruct_dzi.py
diff --git a/PyNutil/main.py b/PyNutil/main.py
index c76bdaf1309e1fccd62506c547b76f7ab2c151f4..b4cd48943012c84bfaa15e47e37478c664b54744 100644
--- a/PyNutil/main.py
+++ b/PyNutil/main.py
@@ -1,10 +1,10 @@
 import json
 import os
 
-from .atlas_loader import load_atlas_data, load_custom_atlas
+from .io.atlas_loader import load_atlas_data, load_custom_atlas
 from .coordinate_processing import extract_coordinates, label_points_group
 from .data_analysis import quantify_labeled_points
-from .file_operations import save_analysis_output
+from .io.file_operations import save_analysis_output
 
 
 class PyNutil: