extract_mismatch_positions: Extract mismatch positions from BAM data

In JakobPedersenLab/dreams: Deep Read-level Error Model for Sequencing data

Extract mismatch positions from BAM data

Description

This function filters a dataframe for reads that contain mismatches when compared to the reference genome, as indicated by the MD tag within the BAM data. It then extracts the genomic positions of these mismatches, providing a detailed account of where the sequencing reads do not match the reference sequence.

Usage

extract_mismatch_positions(bam_df)

Arguments

bam_df

A dataframe obtained from the load_BAM function.

Value

A dataframe of mismatch positions, which includes the genomic positions of all identified mismatches across the reads.

JakobPedersenLab/dreams documentation built on Feb. 2, 2024, 3:14 p.m.