classify_sequences: Given a set of sequences containing mutations and identical...

Description Usage Arguments Value

View source: R/classify.r

Description

The reads_per_index_summary provides a table of the number of mutant and identical reads associated with every index in the data. Identical reads are (obviously) ones with 0 mutations. Sequencer index/reads are when there is 1 mutant out of n reads for an index. RT index/reads are when all reads for an index are mutant.

Usage

1
2
3
4
5
6
7
classify_sequences(
  chng,
  ident,
  reads_per_index_summary,
  min_sequencer = 10,
  verbose = verbose
)

Arguments

chng

Set of reads with mutants.

ident

Set of reads which are identical.

reads_per_index_summary

Table containing how many reads/index are identical/mutant.

min_sequencer

Define how stringent one must be to deem a given read/index group as a sequencer-based error. Too small and gets too many false positives, too large and one gets insufficient data to play with.

verbose

Print some information while running.

Value

Groups of reads which are RT, identical, and sequencer; along with summary information.


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