mbwa_wrapper

mbwa_wrapper aligns short sequencing reads to large reference genomes such as the human genome by leveraging the Burrows-Wheeler Alignment (BWA) algorithm for downstream genomic analyses.

Key Features:

• BWA implementation: Implements the Burrows-Wheeler Alignment (BWA) algorithm that uses the Burrows-Wheeler Transform and backward search for read alignment.
• Speed and efficiency: Provides approximately 10–20× faster alignment than MAQ while maintaining comparable accuracy, enabling large-scale resequencing.
• Read-type support: Supports base space reads from Illumina and color space reads from AB SOLiD.
• Mismatch and gap handling: Accommodates mismatches and gaps (indels) in alignments, suitable for longer reads.
• Output format compatibility: Produces alignments in SAM (Sequence Alignment/Map) format for integration with SAMtools.
• High-throughput suitability: Designed to process large datasets typical of high-throughput next-generation sequencing projects.

Scientific Applications:

• Resequencing projects: Aligns short reads from many individuals against a reference genome for population-scale analyses.
• Variant calling: Generates alignments in SAM format to support downstream variant identification using tools like SAMtools.
• Biomedical research: Enables high-throughput alignment tasks required for genomic studies of disease and biological processes.

Methodology:

Implements the Burrows-Wheeler Alignment (BWA) algorithm employing the Burrows-Wheeler Transform and backward search, permitting mismatches and gaps during alignment.