...

In the ‘2.Setup’ section, you downloaded 3 nfcore/sarek runs to H:/workshop/sarek/runs. These represent sarek variant calling results from the 3 exercises run in session 2:

...

SNPeff creates its own variant analysis reports (html reports) for each sample, which can be seen here:

H:/workshop/sarek/runs/run2_trio/results/reports/snpeff/haplotypecaller

We’ll be creating our of reports, with some additional functional information, using R.

...

Click ‘Knit’ to generate the report. It will ask you to name the file and ask where to save the file. Call it ‘test’ and save it in H:/workshop

...

You will see it generates the test html report.

...

To do this we will use the R clusterProfiler package, which examines enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) terms, based on a set of input genes. In this case, the genes associated with the SNPeff variants.

...

...

Creating the

...

report

The R script and R Markdown scripts have already been created, so you’ll just need to open them in RStudio.

In RStudio, select ‘File’ → “Open File”

Browse to H:/workshop/sarek/runs/run3_liver

Select the ‘modify.R’ file and click ‘Open’.

Open the ‘function_report.Rmd’ file in the same way (it’s in the same location as modify.R)

...

setwd("H:/workshop/sarek/runs/run3_liver")

This sets the working directory you want to run the analysis from. This needs to be the same directory where you ran your sarek script from, i.e. where the sarek ‘results’ directory is.

...

Now just run the rmarkdown::render(.... line to automatically generate the report.

You will now see the HTML report in the working directory you defined in setwd(...). This is a single file, a self-contained report with all the tables and figures you added to it. It can be opened in a web browser, emailed to colleagues or published as a web page.