deBGA-VARA

deBGA-VARA performs variation-aware high-throughput sequencing read alignment by incorporating known genomic variants to improve mapping sensitivity and accuracy.

Key Features:

• Variation-Aware Alignment: Integrates the VARA variation-aware global alignment algorithm into deBGA to align reads against multiple genomic sequences using known genetic variants to improve alignment quality and variant calling accuracy.
• Efficient Memory Usage: Dynamically extracts regional variants on-the-fly and constructs a pseudo tree-based structure for seed extension, avoiding loading the entire genome variation set into memory to reduce RAM usage.
• Algorithmic Foundation: Implements VARA based on the Landau-Vishkin algorithm to enable efficient global alignment by dynamically handling genomic variations during read mapping.
• Benchmarking and Performance: Evaluated with simulated reads and real datasets (e.g., NA12878), demonstrating improved sensitivity and accuracy in read alignment while maintaining manageable memory footprints for variant-calling workflows.

Scientific Applications:

• Variation-aware read mapping: Accurate mapping of high-throughput sequencing reads to references that incorporate known variants.
• Variant calling: Improves sensitivity and accuracy of downstream variant calling from aligned reads.
• Genome-wide association studies (GWAS): Provides higher-quality alignments that can enhance genotype accuracy for GWAS analyses.
• Comparative and evolutionary genomics: Enables more precise identification of genetic variants for evolutionary biology studies.
• Personalized medicine: Facilitates precise identification of genetic variants relevant to personalized medicine analyses.

Methodology:

VARA was integrated into deBGA; deBGA-VARA dynamically extracts regional variants on-the-fly, constructs a pseudo tree-based structure for seed extension, and applies a Landau-Vishkin-based global alignment algorithm; evaluation used simulated reads and the NA12878 dataset.