Usage
processAmplicons(readfile, readfile2=NULL, barcodefile, hairpinfile,
barcodeStart=1, barcodeEnd=5,
barcode2Start=NULL, barcode2End=NULL,
barcodeStartRev=NULL, barcodeEndRev=NULL,
hairpinStart=37, hairpinEnd=57,
allowShifting=FALSE, shiftingBase=3,
allowMismatch=FALSE, barcodeMismatchBase=1,
hairpinMismatchBase=2, allowShiftedMismatch=FALSE,
verbose=FALSE)
Arguments
readfile
character vector giving one or more fastq filenames
readfile2
character vector giving one or more fastq filenames for reverse read, default to NULL
barcodefile
filename containing sample-specific barcode ids and sequences
hairpinfile
filename containing hairpin/sgRNA-specific ids and sequences
barcodeStart
numeric value, starting position (inclusive) of barcode sequence in reads
barcodeEnd
numeric value, ending position (inclusive) of barcode sequence in reads
barcode2Start
numeric value, starting position (inclusive) of second barcode sequence in forward reads
barcode2End
numeric value, ending position (inclusive) of second barcode sequence in forward reads
barcodeStartRev
numeric value, starting position (inclusive) of barcode sequence in reverse reads, default to NULL
barcodeEndRev
numeric value, ending position (inclusive) of barcode sequence in reverse reads, default to NULL
hairpinStart
numeric value, starting position (inclusive) of hairpin/sgRNA sequence in reads
hairpinEnd
numeric value, ending position (inclusive) of hairpin/sgRNA sequence in reads
allowShifting
logical, indicates whether a given hairpin/sgRNA can be matched to a neighbouring position
shiftingBase
numeric value of maximum number of shifted bases from input hairpinStart
and hairpinEnd
should the program check for a hairpin/sgRNA match when allowShifting
is TRUE
allowMismatch
logical, indicates whether sequence mismatch is allowed
barcodeMismatchBase
numeric value of maximum number of base sequence mismatches allowed in a barcode sequence when allowShifting
is TRUE
hairpinMismatchBase
numeric value of maximum number of base sequence mismatches allowed in a hairpin/sgRNA sequence when allowShifting
is TRUE
allowShiftedMismatch
logical, effective when allowShifting
and allowMismatch
are both TRUE
. It indicates whether we check for sequence mismatches at a shifted position.
verbose
if TRUE
, output program progress