Description Usage Arguments Value See Also Examples
View source: R/helpers_plots.R
Plots mismatches in relation to the amplicon, assumes your reads are relative to the respective amplicon sequences predicted cut sites. Top plot is for the forward reads, middle one shows amplicon sequence, and bottom plot is for reverse reads.
1 2 | plot_mismatches(alignments, config, id, cut_buffer = 5,
xlab_spacing = 4)
|
alignments |
(data.frame) Loaded alignment information from alignments_events.csv file. |
config |
(data.frame) Loaded table from config_summary.csv file. |
id |
(string or vector of strings) Name of the ID column from config file or name of multiple IDs, if it is possible to group them. They have to have the same amplicon, amplicons on the reverse strand will be reverse complemented to match forward strand amplicons. |
cut_buffer |
(numeric) Default is 5, you should specify the same as used in the analysis. |
xlab_spacing |
(numeric) Spacing of the x axis labels. Default is 4. |
(mismatches plot) gtable object of mismatches plot
Other specialized plots: metaplot_deletions
,
metaplot_insertions
,
metaplot_mismatches
,
plot_cuts
, plot_deletions
,
plot_heterogeneity
,
plot_insertions
,
plot_variants
1 2 3 4 5 6 7 8 9 | #example config
config <- read.csv(system.file("extdata", "results", "config_summary.csv",
package = "amplican"))
#example alignments results
alignments_file <- system.file("extdata", "results", "alignments",
"events_filtered_shifted_normalized.csv",
package = "amplican")
alignments <- read.csv(alignments_file)
p <- plot_mismatches(alignments, config, c('ID_1', 'ID_3'))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.