Utility description: Variant Call Acceptance
Variant Call Acceptance Checker — Independent Raw Variant Call Verification
ℹ️ Utility performs multi-parameter verification of raw variant calls according to CAP/CLIA, FDA NGS guidance and GATK Best Practices:
• Depth Filtering: Context-dependent depth thresholds (≥20× germline, ≥100× somatic).
• Allele Balance Check: Verification of expected allele ratios for heterozygous and somatic calls.
• Artifact Detection: Identification of strand bias (FS), low QualByDepth (QD), and other positional artifacts.
• Mapping & Base Quality: Control of alignment and sequencing quality at individual read level.
• Control Sample Check: Automatic filtering of variants detected in negative controls.
• Problematic Regions: Flagging of calls in homopolymers and segmental duplications.
⚠️ IMPORTANT:
• Variant caller filters are INSUFFICIENT. Independent verification is mandatory for clinical tests.
• Strand bias (FS > 60) is the most frequent indicator of false-positive SNV.
• Indels in homopolymers have FP-rate up to 50% without additional filters.
• Variant in negative control = systemic issue. All such calls must be rejected.
Usage:
VariantCallAcceptanceChecker.exe → demo mode (console output)
VariantCallAcceptanceChecker.exe input.csv output.json → evaluate your data
Input format:
SampleID,VariantID,VariantType,Context,TotalDepth,MinDepth_Germline,MinDepth_Somatic,AlleleBalance,MinAlleleBalance_Het,MaxAlleleBalance_Het,MappingQuality,MinMappingQuality,BaseQuality,MinBaseQuality,StrandBias_FS,MaxStrandBias_FS,ReadPositionBias_QD,MinQD,GenotypeQuality_GQ,MinGenotypeQuality,InControlSample,InGNomAD_PassFilter,CallerFilter,IsHomopolymerRegion,IsSegmentalDuplication
Example:
WES-001,chr7:55249063:G:T,SNV,Germline,85,20,100,0.48,0.30,0.70,60,40,35,20,2.1,60,18.5,2.0,99,20,false,true,PASS,false,false
📍 Scope of Application (Usage Where):
• Clinical NGS Laboratories: Gate between variant calling and annotation.
• Bioinformatics Pipeline Validation: Confirmation of filter effectiveness.
• Liquid Biopsy: Strict QC of somatic calls at low VAF.
• Research Projects: Ensuring high-quality datasets for publication.
— WHY IS THIS NEEDED?
Variant callers are optimized for sensitivity, not specificity.
Without additional filters, 5-20% of calls are technical artifacts.
Every false-positive variant entering annotation wastes expert time and may lead to erroneous conclusion.
Independent acceptance check is the second safety barrier after the caller itself.
⚠️ CRITICAL:
• Depth ≥ Context-Specific Min: Germline ≥20×, somatic ≥100× (for VAF <5%).
• Strand Bias FS ≤ 60 (germline) / ≤ 200 (somatic): Exceedance = artifact.
• QD ≥ 2.0: Low QD = variant not supported by base quality.
• GQ ≥ 20: Below = uncertain genotype.
• In Control Sample = False: ANY variant in negative control = FILTERED.
• Homopolymer Indels: Require manual review or orthogonal confirmation.
Key features:
• Nine-parameter call quality assessment
• Context-dependent thresholds (somatic vs germline)
• Automatic problematic genomic region detection
• Three-tier classification (Accepted / Review Required / Filtered)
• Compliance with GATK Best Practices and CAP/CLIA
Critical parameters:
• Total Depth: ≥ Context-specific minimum
• Allele Balance: Within expected range
• Mapping Quality: ≥ 40
• Strand Bias (FS): ≤ Threshold
• QualByDepth (QD): ≥ 2.0
• Genotype Quality (GQ): ≥ 20
• Negative Control: Not present
• Caller Filter: PASS (or reviewed)
💡 Usage tips:
1. Variant-Type Thresholds: Indels require stricter FS/QD than SNVs. Configure separately.
2. Panel of Normals (PoN): Use PoN from ≥40 samples to filter recurrent artifacts.
3. Orthogonal Validation: Confirm somatic calls with VAF <5% by ddPCR or amplicon-seq.
4. Artifact Trends: Monitor % filtered by type to detect reagent/instrument degradation.
5. Integration: Run checker automatically after variant calling, before annotation.
⚠️ Note: This utility verifies TECHNICAL QUALITY of the call, not biological significance. A technically quality call may be a benign polymorphism. Biological interpretation is performed at subsequent stages (VariantClassificationGateChecker, VariantReportabilityRulesChecker).
input.csv
SampleID,VariantID,VariantType,Context,TotalDepth,MinDepth_Germline,MinDepth_Somatic,AlleleBalance,MinAlleleBalance_Het,MaxAlleleBalance_Het,MappingQuality,MinMappingQuality,BaseQuality,MinBaseQuality,StrandBias_FS,MaxStrandBias_FS,ReadPositionBias_QD,MinQD,GenotypeQuality_GQ,MinGenotypeQuality,InControlSample,InGNomAD_PassFilter,CallerFilter,IsHomopolymerRegion,IsSegmentalDuplication
WES-2026-001,chr7:55249063:G:T,SNV,Germline,85,20,100,0.48,0.30,0.70,60,40,35,20,2.1,60,18.5,2.0,99,20,false,true,PASS,false,false
WES-2026-001,chr1:115256530:A:AT,Indel,Germline,45,20,100,0.35,0.30,0.70,38,40,22,20,125.0,60,1.2,2.0,15,20,false,false,LowQual,true,false
TUMOR-2026-001,chr12:25398284:C:A,SNV,Somatic,65,20,100,0.08,0.02,0.40,55,40,30,20,8.5,200,12.0,2.0,45,20,false,true,PASS,false,false
WES-2026-002,chr3:37090354:G:A,SNV,Germline,120,20,100,0.50,0.30,0.70,60,40,38,20,1.0,60,22.0,2.0,99,20,true,true,PASS,false,false
URS & FS — User Requirements and Functional Specification
This document describes the controlled interface and behaviour of VariantCallAcceptanceChecker for Variant Call Acceptance Checker.
Domain limits and critical parameters
Key fragments from the source description are shown below. Before production use, limits must be verified against the approved specification, registration dossier and local SOPs.
- • Depth Filtering: Context-dependent depth thresholds (≥20× germline, ≥100× somatic).
- • Strand bias (FS > 60) is the most frequent indicator of false-positive SNV.
- • Variant in negative control = systemic issue. All such calls must be rejected.
- ⚠️ CRITICAL:
- • Depth ≥ Context-Specific Min: Germline ≥20×, somatic ≥100× (for VAF <5%).
- • Strand Bias FS ≤ 60 (germline) / ≤ 200 (somatic): Exceedance = artifact.
- • QD ≥ 2.0: Low QD = variant not supported by base quality.
- • GQ ≥ 20: Below = uncertain genotype.
- Critical parameters:
- • Total Depth: ≥ Context-specific minimum
- • Mapping Quality: ≥ 40
- • Strand Bias (FS): ≤ Threshold
- • QualByDepth (QD): ≥ 2.0
- • Genotype Quality (GQ): ≥ 20
- 2. Panel of Normals (PoN): Use PoN from ≥40 samples to filter recurrent artifacts.
- 3. Orthogonal Validation: Confirm somatic calls with VAF <5% by ddPCR or amplicon-seq.
URS — User Requirements Specification
| ID | Requirement | Criticality | Acceptance criterion |
|---|
| URS-001 | The utility shall accept an input.csv file for Variant Call Acceptance Checker with headers defined in the data contract. | High | The file is processed without manual header editing. |
| URS-002 | The utility shall perform deterministic QC evaluation without machine learning and without probabilistic conformance decisions. | High | Identical input data, rule version and configuration produce reproducible results. |
| URS-003 | The utility shall validate mandatory fields, data types, ranges, units and domain plausibility. | High | Schema, conversion and range errors are explicitly reported. |
| URS-004 | The utility shall apply domain limits and rules from the description, approved specification, registration dossier and local SOPs. | High | Each check has PASS/WARNING/FAIL and a clear message. |
| URS-005 | The utility shall generate output.json with machine-readable results, source values, warnings, failures and critical findings. | High | JSON is suitable for LIMS/ELN/MES integration and QA/QC review. |
| URS-006 | The utility shall preserve traceability between batch/sample, input file, applied rules and final status. | High | Output contains identifiers, checked parameters and audit metadata. |
| URS-007 | The documentation shall support IQ/OQ/PQ, CSV/CSA and review by internal QA or inspectors. | Medium | URS, FS, input/output contract and test scenarios are supplied with the utility. |
| URS-008 | The utility shall be used as a QC decision-support tool and not as a substitute for approved specifications and QA/QP release decision. | Medium | Documentation states change control and limit-verification expectations. |
input.csv contract
| # | Field | Type | Sample | Purpose |
|---|
| 1 | SampleID | string / controlled vocabulary | WES-2026-001 | Sample or laboratory specimen identifier. |
| 2 | VariantID | string / controlled vocabulary | chr7:55249063:G:T | Controlled input parameter for deterministic QC rules. |
| 3 | VariantType | string / controlled vocabulary | SNV | Controlled input parameter for deterministic QC rules. |
| 4 | Context | string / controlled vocabulary | Germline | Controlled input parameter for deterministic QC rules. |
| 5 | TotalDepth | decimal | 85 | Controlled input parameter for deterministic QC rules. |
| 6 | MinDepth_Germline | decimal | 20 | Controlled input parameter for deterministic QC rules. |
| 7 | MinDepth_Somatic | decimal | 100 | Controlled input parameter for deterministic QC rules. |
| 8 | AlleleBalance | decimal | 0.48 | Controlled input parameter for deterministic QC rules. |
| 9 | MinAlleleBalance_Het | decimal | 0.30 | Controlled input parameter for deterministic QC rules. |
| 10 | MaxAlleleBalance_Het | decimal | 0.70 | Controlled input parameter for deterministic QC rules. |
| 11 | MappingQuality | decimal | 60 | Controlled input parameter for deterministic QC rules. |
| 12 | MinMappingQuality | decimal | 40 | Controlled input parameter for deterministic QC rules. |
| 13 | BaseQuality | decimal | 35 | Controlled input parameter for deterministic QC rules. |
| 14 | MinBaseQuality | decimal | 20 | Controlled input parameter for deterministic QC rules. |
| 15 | StrandBias_FS | decimal | 2.1 | Controlled input parameter for deterministic QC rules. |
| 16 | MaxStrandBias_FS | decimal | 60 | Controlled input parameter for deterministic QC rules. |
| 17 | ReadPositionBias_QD | decimal | 18.5 | Controlled input parameter for deterministic QC rules. |
| 18 | MinQD | decimal | 2.0 | Controlled input parameter for deterministic QC rules. |
| 19 | GenotypeQuality_GQ | string / controlled vocabulary | 99 | Controlled input parameter for deterministic QC rules. |
| 20 | MinGenotypeQuality | decimal | 20 | Controlled input parameter for deterministic QC rules. |
| 21 | InControlSample | string / controlled vocabulary | false | Sample or laboratory specimen identifier. |
| 22 | InGNomAD_PassFilter | string / controlled vocabulary | true | Controlled input parameter for deterministic QC rules. |
| 23 | CallerFilter | string / controlled vocabulary | PASS | Controlled input parameter for deterministic QC rules. |
| 24 | IsHomopolymerRegion | string / controlled vocabulary | false | Controlled input parameter for deterministic QC rules. |
| 25 | IsSegmentalDuplication | string / controlled vocabulary | false | Controlled input parameter for deterministic QC rules. |
SampleID,VariantID,VariantType,Context,TotalDepth,MinDepth_Germline,MinDepth_Somatic,AlleleBalance,MinAlleleBalance_Het,MaxAlleleBalance_Het,MappingQuality,MinMappingQuality,BaseQuality,MinBaseQuality,StrandBias_FS,MaxStrandBias_FS,ReadPositionBias_QD,MinQD,GenotypeQuality_GQ,MinGenotypeQuality,InControlSample,InGNomAD_PassFilter,CallerFilter,IsHomopolymerRegion,IsSegmentalDuplication
WES-2026-001,chr7:55249063:G:T,SNV,Germline,85,20,100,0.48,0.30,0.70,60,40,35,20,2.1,60,18.5,2.0,99,20,false,true,PASS,false,false
WES-2026-001,chr1:115256530:A:AT,Indel,Germline,45,20,100,0.35,0.30,0.70,38,40,22,20,125.0,60,1.2,2.0,15,20,false,false,LowQual,true,false
TUMOR-2026-001,chr12:25398284:C:A,SNV,Somatic,65,20,100,0.08,0.02,0.40,55,40,30,20,8.5,200,12.0,2.0,45,20,false,true,PASS,false,false
Input validation rules
| ID | Field | Rule | Criticality |
|---|
| VR-001 | SampleID | The field shall match an approved dictionary or accepted string representation. | High |
| VR-002 | VariantID | The field shall match an approved dictionary or accepted string representation. | High |
| VR-003 | VariantType | The field shall match an approved dictionary or accepted string representation. | High |
| VR-004 | Context | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-005 | TotalDepth | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-006 | MinDepth_Germline | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-007 | MinDepth_Somatic | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-008 | AlleleBalance | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-009 | MinAlleleBalance_Het | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-010 | MaxAlleleBalance_Het | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-011 | MappingQuality | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-012 | MinMappingQuality | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-013 | BaseQuality | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-014 | MinBaseQuality | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-015 | StrandBias_FS | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-016 | MaxStrandBias_FS | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-017 | ReadPositionBias_QD | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-018 | MinQD | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-019 | GenotypeQuality_GQ | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-020 | MinGenotypeQuality | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-021 | InControlSample | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-022 | InGNomAD_PassFilter | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-023 | CallerFilter | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-024 | IsHomopolymerRegion | The field shall match an approved dictionary or accepted string representation. | Medium |
| VR-025 | IsSegmentalDuplication | The field shall match an approved dictionary or accepted string representation. | Medium |
FS — Functional Specification
| ID | Function | Implementation |
|---|
| FS-001 | CLI execution | Support execution modes: demo mode without arguments and production mode input.csv output.json. |
| FS-002 | CSV import | Read input.csv in UTF-8/CSV-compatible format and validate header and expected columns. |
| FS-003 | Schema validation | Check mandatory fields, column count, unknown key fields and empty mandatory values. |
| FS-004 | Type conversion | Convert numeric, flag and text values; invalid format is recorded as a row-level error. |
| FS-005 | Domain rule engine | Apply rules for Variant Call Acceptance Checker, including critical limits from the description and approved specification. |
| FS-006 | Status aggregation | Produce final status: FAIL for critical failure, WARNING for non-critical deviation, PASS for conformance. |
| FS-007 | JSON export | Write output.json with detailed checks, source values, warnings, failures and critical findings. |
| FS-008 | Audit support | Keep result structure suitable for review, deviation investigation and calculation reproduction. |
| FS-009 | Integration contract | Support the scenario LIMS/ELN/MES → input.csv → utility → output.json → portal/admin review. |
| FS-010 | Error handling | Return explicit messages for missing file, empty CSV, invalid schema, output write failure and invalid format. |
Example output.json
{
"utilityId": "variantcallacceptancechecker",
"utilityFolder": "VariantCallAcceptanceChecker",
"package": "LiquidBiopsy",
"overallStatus": "PASS|WARNING|FAIL",
"sourceFile": "input.csv",
"processedAtUtc": "2026-06-10T00:00:00Z",
"checks": [
{
"parameter": "SampleID",
"value": "WES-2026-001",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-001"
},
{
"parameter": "VariantID",
"value": "chr7:55249063:G:T",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-002"
},
{
"parameter": "VariantType",
"value": "SNV",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-003"
},
{
"parameter": "Context",
"value": "Germline",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-004"
},
{
"parameter": "TotalDepth",
"value": "85",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-005"
},
{
"parameter": "MinDepth_Germline",
"value": "20",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-006"
},
{
"parameter": "MinDepth_Somatic",
"value": "100",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-007"
},
{
"parameter": "AlleleBalance",
"value": "0.48",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-008"
},
{
"parameter": "MinAlleleBalance_Het",
"value": "0.30",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-009"
},
{
"parameter": "MaxAlleleBalance_Het",
"value": "0.70",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-010"
},
{
"parameter": "MappingQuality",
"value": "60",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-011"
},
{
"parameter": "MinMappingQuality",
"value": "40",
"status": "PASS|WARNING|FAIL",
"message": "Deterministic rule-based check result",
"ruleReference": "FS-RULE-012"
}
],
"criticalFindings": [],
"warnings": [],
"audit": {
"inputHash": "sha256:<calculated at runtime>",
"rulesVersion": "<utility executable version>",
"documentation": "VariantCallAcceptanceChecker.documentation.html"
}
}
Traceability matrix
| URS | FS | Test | Evidence |
|---|
| URS-001 | FS-001, FS-002 | OQ-001 | Verify execution and import of valid input.csv. |
| URS-002 | FS-005, FS-006 | OQ-004 | Repeat the same dataset and compare output.json. |
| URS-003 | FS-003, FS-004, FS-010 | OQ-002, OQ-003 | Verify missing columns and invalid types. |
| URS-004 | FS-005, FS-006 | OQ-004, PQ-001 | Verify critical deviations on real/boundary data. |
| URS-005 | FS-007, FS-009 | OQ-005 | Verify JSON schema and downstream-system suitability. |
| URS-006 | FS-008 | OQ-006 | Verify identifiers and audit metadata. |
| URS-007 | FS-008, FS-010 | IQ-001, OQ-007 | Verify documentation completeness and control evidence. |
| URS-008 | FS-005, FS-008 | PQ-002 | Verify review workflow and no replacement of QA decision. |
IQ/OQ/PQ test scenarios
| ID | Scenario | Expected result |
|---|
| IQ-001 | Verify executable, input.csv, documentation and checksum availability. | Delivery set is complete; version is recorded. |
| OQ-001 | Valid sample row from input.csv. | PASS or acceptable WARNING according to rules. |
| OQ-002 | Remove a mandatory CSV column. | Schema error or FAIL with missing-column reference. |
| OQ-003 | Place a non-numeric value into a numeric field. | Type-conversion error with row/field reference. |
| OQ-004 | Set a critical parameter outside the limit. | FAIL and critical finding. |
| OQ-005 | Verify output.json structure. | All mandatory sections are present and JSON is valid. |
| OQ-006 | Verify batch/sample traceability. | Input and result identifiers match. |
| PQ-001 | Verify 3–5 real user batches/samples. | Result is confirmed by QC/QA review. |
| PQ-002 | Verify deviation workflow and manual QA decision. | Utility supports review but does not replace approved decision. |
QA/QC and change control
- Do not rename columns without updating validator, documentation and test set.
- Retain
input.csv, output.json, executable version and checksum. - Before production use, perform IQ/OQ/PQ or equivalent CSV/CSA verification.
- Critical limits shall be verified against the approved specification, registration dossier and local SOPs.
- The utility provides structured QC decision support; final release decision remains with QA/QP and approved procedures.