Megadepth

Megadepth quantifies sequencing coverage across genomic regions from BigWig, BAM, and CRAM files for interval-level and base-pair coverage analyses.

Key Features:

• Efficiency in data processing: Implements high-performance processing with reduced memory usage and parallel threading, demonstrated by summarizing coverage for all disjoint intervals of the Gencode V35 annotation across over 19,000 GTExV8 BigWig files in approximately one hour using 32 threads.
• Versatile input compatibility: Supports BigWig and BAM/CRAM file formats for unified quantification of alignments and coverage.
• Comprehensive coverage analysis: Computes coverage for specific genomic regions or annotations and can produce base-pair coverage summaries for downstream analyses.

Scientific Applications:

• Interval-level coverage summarization: Summarizes sequencing coverage across annotated genomic intervals such as gene models from Gencode V35 for large cohorts and annotation sets.
• Base-pair coverage matrix generation: Produces base-pair resolution coverage summaries suitable for matrix-based downstream analyses.
• Support for downstream genomics analyses: Facilitates downstream tasks including differential expression studies, variant calling, and functional genomics investigations.

Methodology:

Employs optimized algorithms with efficient memory management and parallel processing capabilities to perform fast and accurate coverage quantification.