Skip to content

Commit

Permalink
Add hide/show buttons for the example solutions of sbatch scripts
Browse files Browse the repository at this point in the history
  • Loading branch information
malinlarsson committed Nov 11, 2024
1 parent f062fa1 commit e87e202
Showing 1 changed file with 16 additions and 3 deletions.
19 changes: 16 additions & 3 deletions topics/vc/lab_vc.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -766,7 +766,7 @@ To learn more about SLURM and SBATCH scripts please look the [How to Run Jobs](h
## Variant calling in cohort {-}

Below is a skeleton script that can be used as a template for running [variant calling in a cohort](#jointvc). Please modify it to run all the steps in part two of this workshop.

```{r,accordion=TRUE,echo=FALSE,comment="",class.output="sh"}
cat(paste0('#!/bin/bash
#SBATCH -A ',id_project,'
Expand All @@ -792,7 +792,7 @@ done
#Fill in the code for GenotypeGVCFs here
'))
```
Please save the sbatch script in your Dardel folder and call it "joint_genotyping.sbatch" or similar. Make the script executable by this command:
Please save the template in your workspace and call it "joint_genotyping.sbatch" or similar. Make the script executable by this command:

```bash
chmod u+x joint_genotyping.sbatch
Expand All @@ -810,7 +810,10 @@ If you have an active node reservation you can run the script as a normal bash s
./joint_genotyping.sbatch
```

If you would like more help with creating the sbatch script, please look at our example solution:
If you would like more help with creating the sbatch script, please look at our example solution below (click on the link):

<details>
<summary>Example solution</summary>

```{r,accordion=TRUE,echo=FALSE,comment="",class.output="sh"}
cat(paste0('#!/bin/bash
Expand Down Expand Up @@ -854,10 +857,17 @@ gatk --java-options -Xmx2g GenotypeGVCFs \\
-O cohort.vcf'))
```


</details>


## GATK best practices {-}

Now please try to incorporate the additional steps from [GATK's best practices](#gatk_bp) into the workflow. If you run out of time you can sneak peek at our example solution below.

<details>
<summary>Example solution</summary>

```{r,accordion=TRUE,echo=FALSE,comment="",class.output="sh"}
cat(paste0('#!/bin/bash
#SBATCH -A ',id_project,'
Expand Down Expand Up @@ -978,6 +988,9 @@ gatk --java-options -Xmx2g MergeVcfs \\
-O cohort.filtered.vcf'))
```


</details>

# Additional information {-}

* Here is a technical documentation of [Illumina Quality Scores](https://www.illumina.com/documents/products/technotes/technote_understanding_quality_scores.pdf)
Expand Down

0 comments on commit e87e202

Please sign in to comment.