list_mismatches: Lists mismatches between two DNAstrings

In rocioespci/single: Accurate consensus sequence from nanopore reads of a gene library

Lists mismatches between two DNAstrings

Description

This is an auxiliary function in single package, to list the mutations of two DNAstrings.

Usage

list_mismatches(ref, seq)

Arguments

ref	DNAString, reference sequence.
seq	DNAString, target sequence, same length as ref.

Value

Character vector containing Nucleotide in ref Position Nucleotide in seq. If ref and seq are equal, it returns NA.

Examples

ref = Biostrings::DNAString("AAAA")
seq = Biostrings::DNAString("AGAT")
list_mismatches(ref,seq)
list_mismatches(ref,ref)

rocioespci/single documentation built on April 18, 2023, 8:48 p.m.