Question

In the context of SAM/BAM file manipulation, which specific index-based query operation allows a tool to retrieve only the reads mapped to a specific chromosomal region without reading the entire file?

Accepted Answer

The specific operation is known as a coordinate-based range query, which utilizes a BAI or CSI index file to perform random access on a sorted BAM file. A BAM file is a binary version of a SAM file, which stores high-throughput DNA sequencing data. Because these files are often gigabytes or terabytes in size, reading the entire file to find reads in a specific location is computationally inefficient. To enable fast retrieval, the BAM file must first be sorted by coordinate, meaning the reads are arranged according to their position on the reference genome. A companion index file, usually ending in .bai, is then generated. This index acts like a book index by storing pointers or byte offsets that map specific genomic intervals to precise locations within the compressed binary file. When a user requests a specific region, such as chromosome 1 from base 1000 to 2000, the tool reads the index to determine exactly which blocks of the compressed BAM file contain the relevant data. The tool then jumps directly to those specific byte offsets in the BAM file and reads only that small portion, bypassing the rest of the data entirely. This process allows software like samtools or bedtools to extract specific reads instantly without loading the entire genomic dataset into memory.

Home → All Courses → Science and Mathematics Courses → Genomic Data Science → Flashcard

In the context of SAM/BAM file manipulation, which specific index-based query operation allows a tool to retrieve only the reads mapped to a specific chromosomal region without reading the entire file?