#sequence length refers to sequence length AFTER removal of Primers, Barcodes and trimming. this ensures that downstream analyis tools will have appropiate sequence information
minSeqLength	170
maxSeqLength	1000
minAvgQuality	28
#*minSeqLength	170
#*minAvgQuality	20
#truncate total Sequence length to X (length after Barcode, Adapter and Primer removals)
TruncateSequenceLength	170

#Ambiguous bases in Sequence - uclust only supports 0 ambiguous nucleotides
maxAmbiguousNT	0
#*maxAmbiguousNT	1

#Binomial error model of expected errors per sequence (see https://github.com/fpusan/moira), to deactivate, set BinErrorModelAlpha to -1
BinErrorModelMaxExpError	2.5
BinErrorModelAlpha	0.01

#Homonucleotide Runs.. this should normally be filtered by sequencer software
maxHomonucleotide	12

#Filter whole sequence if one window of quality scores is below average
QualWindowWidth	50
QualWindowThreshhold	25

#Trim the end of a sequence if a window falls below quality threshhold. Useful for removing low qulaity trailing ends of sequence
TrimWindowWidth	20
TrimWindowThreshhold	25

#Max number of accumulated P for a mismatch. After this length, the rest of the sequence will be deleted. Complimentary to TrimWindowThreshhold. (-1) deactivates this option.
maxAccumulatedError	1
#not implemented
#*maxAccumulatedError	-1


#Barcode Errors - currently this can only be 0; 
maxBarcodeErrs	1
maxPrimerErrs	1

#keep Barcode / Primer Sequence in the output fasta file - in a normal 16S analysis this should be deactivated (0) for Barcode and de-activated (0) for primer
keepBarcodeSeq	0
keepPrimerSeq	0


#set fastqVersion to 1 if you use Sanger, Illumina 1.8+ or NCBI SRA files. Set fastqVersion to 2, if you use Illumina 1.3+ - 1.7+ or Solexa fastq files.
fastqVersion	1

#if one or more files have a technical adapter still included (e.g. TCAG 454) this can be removed by setting this option
TechnicalAdapter	

#delete X NTs (e.g. if the first 5 bases are known to have strange biases)
TrimStartNTs	0

#correct PE header format (1/2) this is to accomodate the illumina miSeq paired end annotations 2="@XXX 1:0:4" insteand of 1="@XXX/1". Note that the format will be automatically detected
PEheaderPairFmt	1

#sets if sequences without match to reverse primer (ReversePrimer) will be accepted (T=reject ; F=accept all); default=F
RejectSeqWithoutRevPrim	F
#*RejectSeqWithoutRevPrim	F
#sets if sequences without a forward (LinkerPrimerSequence) primer will be accepted (T=reject ; F=accept all); default=F
RejectSeqWithoutFwdPrim	F
#*RejectSeqWithoutFwdPrim	T

#checks if the whole amplicon was reverse-transcribed sequenced 
CheckForReversedSeqs	T
CheckForMixedPairs	T


#this option should be "T" if your amplicons are possibly shorter than a read in a paired end sequencing run (e.g. amplicon of 300 in 250x2 miSeq is "T")
AmpliconShortPE	T

