makeAlignment: Make alignments helper.

In valenlab/amplican: Automated analysis of CRISPR experiments

Make alignments helper.

Description

Aligning reads to the amplicons for each ID in this barcode, constructing amplicanAlignment. Assume that all IDs here belong to the same barcode.

Usage

makeAlignment(
  cfgT,
  average_quality,
  min_quality,
  filter_n,
  batch_size,
  scoring_matrix,
  gap_opening,
  gap_extension,
  fastqfiles,
  primer_mismatch,
  donor_mismatch,
  donor_strict
)

Arguments

cfgT	config file as data table
average_quality	(numeric) The FASTQ file have a quality for each nucleotide, depending on sequencing technology there exist many formats. This package uses readFastq to parse the reads. If the average quality of the reads fall below value of average_quality then sequence is filtered. Default is 0.
min_quality	(numeric) Similar as in average_quality, but depicts the minimum quality for ALL nucleotides in given read. If one of nucleotides has quality BELLOW min_quality, then the sequence is filtered. Default is 20.
filter_n	(boolean) Whether to filter out reads containing N base.
batch_size	(numeric) How many reads to analyze at a time? Needed for filtering of large fastq files.
scoring_matrix	(matrix) Default is 'NUC44'. Pass desired matrix using nucleotideSubstitutionMatrix.
gap_opening	(numeric) The opening gap score.
gap_extension	(numeric) The gap extension score.
fastqfiles	(numeric) Normally you want to use both FASTQ files. But in some special cases, you may want to use only the forward file, or only the reverse file. Possible options: 0 Use both FASTQ files. 0.5 Use both FASTQ files, but only for one of the reads (forward or reverse) is required to have primer perfectly matched to sequence - eg. use when reverse reads are trimmed of primers, but forward reads have forward primer in the sequence. 1 Use only the forward FASTQ file. 2 Use only the reverse FASTQ file.
primer_mismatch	(numeric) Decide how many mismatches are allowed during primer matching of the reads, that groups reads by experiments. When primer_mismatch = 0 no mismatches are allowed, which can increase number of unasssigned read.
donor_mismatch	(numeric) How many events of length 1 (mismatches, deletions and insertions of length 1) are allowed when aligning toward the donor template. This parameter is only used when donor template is specified. The higher the parameter the less strict will be algorithm accepting read as HDR. Set to 0 if only perfect alignments to the donor template marked as HDR, unadvised due to error rate of the sequencers.
donor_strict	(logical) Applies more strict algorithm for HDR detection. Only these reads that have all of the donor events will count as HDR. Tolerates 'donor_mismatch' level of noise, but no indels are allowed. Use this when your reads should span over the whole window of the donor events. Might be more time consuming.

Value

amplicanAlignment object for this barcode experiments

valenlab/amplican documentation built on Jan. 28, 2024, 5:10 a.m.