somAgg code book¶
This code book provides some sample snippets to help you use somAgg in your analyses. These include using BEDtools to find the correct chunk file to use, and using BCFtools to query the aggregate files themselves.
Overview¶
The code snippets assume that you are working in the HPC environment and that you submit jobs to the cluster. Please see In-Depth Guide to HPC Usage for more information.
Feedback and requests
For any feedback and requests to the somAgg code book, or if you encounter issues running one of the examples, please reach out via the Genomics England Service Desk including "somAgg" in the title/description of your inquiry.
Applications¶
The majority of queries to aggV2 can be implemented using the applications below:
| Application | Description |
|---|---|
| bcftools | A set of utilities that manipulate variant calls in the Variant Call Format (VCF). Use version 1.16 via module load bcftools/1.16 |
| split-vep | A bcftools plug-in to parse VEP annotation (comes with bcftools version 1.16. |
| LabKey APIs | The LabKey client libraries (APIs) provide programmatic access to the clinical/phenotype data. |
| R / Python | For downstream processing. |
| bedtools | To intersect, merge, count, complement, and shuffle genomic intervals. Use version 2.31.0 via module load bedtools/2.31.1 |
Code book structure¶
We have divided the code book into the following sections: