diff --git a/subworkflows/local/pairalign_m2m/main.nf b/subworkflows/local/pairalign_m2m/main.nf
index 39e40f8..479a7fd 100644
--- a/subworkflows/local/pairalign_m2m/main.nf
+++ b/subworkflows/local/pairalign_m2m/main.nf
@@ -122,6 +122,7 @@ workflow PAIRALIGN_M2M {
     m2o = LAST_SPLIT_M2O.out.maf
     o2m = LAST_SPLIT_O2M.out.maf
     o2o = LAST_SPLIT_O2O.out.maf
+    versions = LAST_LASTDB.out.versions
 }
 
 /*
diff --git a/subworkflows/local/pairalign_m2o/main.nf b/subworkflows/local/pairalign_m2o/main.nf
index 24e90fa..ff244b1 100644
--- a/subworkflows/local/pairalign_m2o/main.nf
+++ b/subworkflows/local/pairalign_m2o/main.nf
@@ -82,6 +82,7 @@ workflow PAIRALIGN_M2O {
         .mix(LAST_SPLIT_O2O.out.multiqc.collect{ it[1]} )
     m2o = LAST_LASTAL_M2O.out.maf
     o2o = LAST_SPLIT_O2O.out.maf
+    versions = LAST_LASTDB.out.versions
 }
 
 /*
diff --git a/workflows/pairgenomealign.nf b/workflows/pairgenomealign.nf
index a7740e5..f84ae3d 100644
--- a/workflows/pairgenomealign.nf
+++ b/workflows/pairgenomealign.nf
@@ -53,7 +53,6 @@ workflow PAIRGENOMEALIGN {
     ASSEMBLYSCAN (
         ch_samplesheet
     )
-    ch_versions = ch_versions.mix(ASSEMBLYSCAN.out.versions.first())
 
     // Prefix query ids with target genome name before producing alignment files
     ch_samplesheet = ch_samplesheet
@@ -88,6 +87,12 @@ workflow PAIRGENOMEALIGN {
 
     // Collate and save software versions
     //
+
+    ch_versions = ch_versions
+        .mix(SEQTK_CUTN_TARGET.out.versions)
+        .mix(     ASSEMBLYSCAN.out.versions)
+        .mix(        pairalign_out.versions)
+
     softwareVersionsToYAML(ch_versions)
         .collectFile(
             storeDir: "${params.outdir}/pipeline_info",