Function get_cram_reads 

pub fn get_cram_reads(
    config: &BamConfig,
    query: impl Iterator<Item = Result<RecordBuf>>,
    region: &Region,
    lflank: usize,
    rflank: usize,
) -> Result<(Vec<BamRead>, usize)>

Extract subsequences from all CRAM records that overlap a given region.

Works identically to get_bam_reads but consumes any iterator that yields io::Result<sam::alignment::RecordBuf> — the item type produced by cram::io::IndexedReader::query. Returns (reads, candidates_seen); see get_bam_reads for what candidates_seen means.