LevioSAM

LevioSAM translates SAM/BAM read alignments between reference coordinate systems using Variant Call Format (VCF) files to apply population-specific variants for coordinate lifting.

Key Features:

• VCF-Based Lifting: Uses VCF files containing population-specific variant information instead of chain files to lift alignments between references.
• Efficient Multithreading: Employs multiple threads and succinct data structures to reduce post-alignment computational bottlenecks and scale with increased threading.
• Performance Optimization: When run downstream of a read aligner on 16 threads, completes the lifting process in less than 13% of the time required by the aligner under similar conditions.

Scientific Applications:

• Population genetics: Translates alignments between different reference genomes to support population genetics analyses.
• Comparative genomics: Enables comparison of alignments across reference assemblies for comparative genomics studies.
• Variant analysis across populations: Facilitates analysis of variants across populations by remapping alignments to population-specific references.
• Integration of diverse genetic studies: Supports integrating alignment data from diverse genetic studies mapped to different references.

Methodology:

Reads SAM/BAM alignment files and uses population-specific variants encoded in VCF files to map alignments onto a reference genome, employing succinct/advanced data structures and multithreading for scalability and speed.