add snp, which was missing by accident

pdimens · Apr 25, 2024 · 749c205 · 749c205
1 parent fc23288
commit 749c205
Show file tree

Hide file tree

Showing 2 changed files with 69 additions and 0 deletions.
diff --git a/.github/workflows/tests.yml b/.github/workflows/tests.yml
@@ -229,6 +229,74 @@ jobs:
         shell: micromamba-shell {0}
         run: harpy align minimap -g test/genome/genome.fasta.gz -s "--show-failed-logs" -x "--seed 13" test/fastq
 
+  test_mpileup:
+    needs: [changes, pkgbuild]
+    if: ${{ needs.changes.outputs.mpileup == 'true' && needs.pkgbuild.result == 'success' }}
+    name: test mpileup
+    runs-on: ubuntu-latest
+    steps:
+      - name: checkout
+        uses: actions/checkout@v4
+        with:
+          fetch-depth: 1
+      - name: setup mamba
+        uses: mamba-org/setup-micromamba@v1
+        env:
+          ACTIONS_STEP_DEBUG: true
+        with:
+          init-shell: bash
+          generate-run-shell: true
+          environment-file: workflow/envs/harpy.yaml
+          cache-environment: true
+          post-cleanup: 'all'
+          log-level: error
+      - name: Install harpy
+        shell: micromamba-shell {0}
+        run: |
+          python3 -m pip install --upgrade build && python3 -m build
+          pip install dist/*.whl
+          resources/buildforCI.sh
+      - name: snp mpileup
+        shell: micromamba-shell {0}
+        run: harpy snp mpileup -w 150000 -g test/genome/genome.fasta.gz -q -x "--ignore-RG" test/bam
+      - name: snp mpileup-pop
+        shell: micromamba-shell {0}
+        run: harpy snp mpileup -w 150000 -g test/genome/genome.fasta.gz -p test/samples.groups -q test/bam
+
+  test_freebayes:
+    needs: [changes, pkgbuild]
+    if: ${{ needs.changes.outputs.freebayes == 'true' && needs.pkgbuild.result == 'success' }}
+    name: test freebayes
+    runs-on: ubuntu-latest
+    steps:
+      - name: checkout
+        uses: actions/checkout@v4
+        with:
+          fetch-depth: 1
+      - name: setup mamba
+        uses: mamba-org/setup-micromamba@v1
+        env:
+          ACTIONS_STEP_DEBUG: true
+        with:
+          init-shell: bash
+          generate-run-shell: true
+          environment-file: workflow/envs/harpy.yaml
+          cache-environment: true
+          post-cleanup: 'all'
+          log-level: error
+      - name: Install harpy
+        shell: micromamba-shell {0}
+        run: |
+          python3 -m pip install --upgrade build && python3 -m build
+          pip install dist/*.whl
+          resources/buildforCI.sh
+      - name: snp freebayes
+        shell: micromamba-shell {0}
+        run: harpy snp freebayes -w 150000 -g test/genome/genome.fasta.gz -q -x "-g 200" test/bam
+      - name: snp freebayes-pop
+        shell: micromamba-shell {0}
+        run: harpy snp freebayes -w 150000 -g test/genome/genome.fasta.gz -p test/samples.groups -q test/bam
+
   test_impute:
     needs: [changes, pkgbuild]
     if: ${{ needs.changes.outputs.impute == 'true' && needs.pkgbuild.result == 'success' }}

diff --git a/src/harpy/rules/demultiplex.gen1.smk b/src/harpy/rules/demultiplex.gen1.smk
@@ -212,6 +212,7 @@ rule log_workflow:
     message:
         "Summarizing the workflow: {output}"
     run:
+        os.makedirs(f"{outdir}workflow/", exist_ok= True)
         with open(outdir + "workflow/demux.gen1.summary", "w") as f:
             _ = f.write("The harpy demultiplex module ran using these parameters:\n\n")
             _ = f.write("Haplotag technology: Generation I\n")