Creating an external database VCF file
Prerequisites
Format: Database file must follow VCF 4.2 specifications.
Tools:
Required:
awkDocumentationbgzipDocumentationtabixDocumentation
Optional:
vcftoolsUseful for population frequency calculations. Documentation
How to create a noise or historic database file
Calculate population statistics
General: Allele Number (AN): Calculate the total number of alleles in your population by multiplying the number of individuals () by 2: .
For each variant within the dataset:
Allele Count (AC): Determine the number of times alternate allele of a variant appears across all individuals. This is inferred from genotype counts (). In case of a biallelic variant where allele A is reference and allele B is alternate, the allele count for alternate allele is calculated as follows:
Allele Frequency (AF): Calculate as the ratio of Allele Count to Allele Number: .
Create a VCF file with your variants
In the
INFOfield, include theAFsub-field.Optionally, include
AC,AN, andTEN(a list of up to 10 samples carrying the variant). You may add other fields, provided the field names do not contain underscores or hyphens.Specify the exact format of each
INFOsub-field in ##INFO meta-information lines.
See format example below.
Sort variants in the VCF based on chromosome and position with awk
awk '$1 ~ /^#/ {print $0;next} {print $0 | "sort -k1,1 -k2,2n"}'Compress the VCF with bgzip
bgzip <your_db>.vcfCreate a TBI index file with tabix
tabix -p vcf <your_db>.vcf.gzExample historic DB VCF header and variant line
How to create a curated database file
Create a VCF file with your variants
In the
INFOfield, include thesignificancesub-field and assign its value to each variant based on Table 1. Only onesignificancevalue is allowed per variant. If a variant has multiple interpretations, list the variant in separate rows, each with a differentsignificancevalue.Table 1. Mapping of significance values to pathogenicity classes.
0
Unknown
1
Benign
2
Likely Benign
3
VUS
4
Likely Pathogenic
5
Pathogenic
Optionally, include
comment,category, or other fields to capture text or numerical values that are relevant to classification. You may add other fields, provided the field names do not contain underscores or hyphens.Specify the exact format of each
INFOsub-field in ##INFO meta-information lines.
See format example below.
Sort variants in the VCF based on chromosome and position with awk
Compress the VCF with bgzip
Create a TBI index file with tabix
Example curated DB VCF variant lines
Small variant
Copy number variant
Next steps
Reach out to Illumina support
Provide VCF and TBI files to Illumina support to upload to your organization's dedicated storage bucket, along with information:
Database name. Underscores ("_") in a database name are not allowed.
Database type (Noise, Historic, Curated)
Variant type (SNV, CNV)
Genome reference (GRCh37, GRCh38)
Register the database
Once the database is uploaded, the user with appropriate permissions can register the database by selecting it from their bucket in Settings.
Last updated
Was this helpful?

