Skip to content

Commit d25a086

Browse files
authored
Update variant_calling_gatk1.md
1 parent 5eae9c3 commit d25a086

File tree

1 file changed

+11
-3
lines changed

1 file changed

+11
-3
lines changed

docs/tutorials/variant_calling_gatk1/variant_calling_gatk1.md

+11-3
Original file line numberDiff line numberDiff line change
@@ -18,7 +18,7 @@ For queries relating to this workshop, contact Melbourne Bioinformatics (bioinfo
1818
Khalid Mahmood
1919
Melbourne Bioinformatics, University of Melbourne
2020
Developed: July 2021
21-
Reviewed: June 2023
21+
Reviewed: August 2024
2222
-------------------------------
2323
## Overview
2424

@@ -629,7 +629,15 @@ gatk VariantsToTable \
629629
```
630630
631631
### 2. HTML report
632-
Another useful method for sharing data is an interactive HTML file. This is suited for sharing a smaller subset of variants along with sequencing data. Here we will go through a simple example using the [jigv](https://github.com/brentp/jigv) tool.
632+
Another useful method for sharing data is an interactive HTML file. This is suited for sharing a smaller subset of variants along with sequencing data. Here we will go through a simple example using the [jigv](https://github.com/brentp/jigv) tool. See below on how to download and setup jigv.
633+
634+
```bash
635+
# download jigv
636+
wget https://github.com/brentp/jigv/releases/download/v0.1.10/jigv
637+
638+
# make jigv executable
639+
chmod +x jigv
640+
```
633641
634642
![fig2](./media/fig2.png)
635643
@@ -645,7 +653,7 @@ tabix -p vcf output/subset.vcf.gz
645653
646654
Now, we will call the jigv tool command to generate the report.
647655
```bash
648-
jigv --sample NA12878 \
656+
./jigv --sample NA12878 \
649657
--sites output/subset.vcf.gz \
650658
--fasta reference/hg38/Homo_sapiens_assembly38.fasta \
651659
output/NA12878.sort.dup.bqsr.bam > output/NA12878.jigv.html

0 commit comments

Comments
 (0)