sequencer_error: Calculate a lower limit error rate for the sequencer.

In abelew/Rerrrt: Examine RT error rates through sequencing.

Description

I say 'lower limit' because I assume my set of sequencer-based errors did not detect the full set of error actually from the sequencer. I chose to take only the set with >= 10 reads / index for which only n (1) read(s) is a mutant. My hope is that this will avoid false positives, but it will also limit the perceived sequencer error rate.

Usage

sequencer_error(data_summary)

Arguments

data_summary

Result from create_matrices().

Value

Estimate of errors/nucleotide deemed to originate from the sequencer.

abelew/Rerrrt documentation built on Jan. 15, 2022, 8 a.m.