pcr.seqs

The pcr.seqs will trim inputted sequences based on a variety of user-defined options.

Default Settings

The fasta parameter is required.

mothur > pcr.seqs(fasta=silva.bacteria.fasta)


Options

oligos

The oligos file allows you to provide primer information.

mothur > pcr.seqs(fasta=silva.bacteria.fasta, oligos=pcrTest.oligos)


pdiffs

The pdiffs parameter is used to allow differences in the forward primers. pdiffs is maximum number of differences to the forward primer sequence, default=0.

rdiffs

The rdiffs parameter is used to allow differences in the reverse primers. rdiffs is maximum number of differences to the reverse primer sequence, default=0.

checkorient

The checkorient parameter will look for the reoriented and reverse compliment of the primers if the primers are not found. The default is true.

name && group && count && taxonomy

If you apply pcr.seqs as we have in the previous sections and then attempt to run any downstream commands that need a name, group and taxonomy file, the files will be incompatible since those files will still contain information for sequences that you culled. To get around this you can use the group and name or count option:

mothur > pcr.seqs(fasta=silva.bacteria.fasta, taxonomy=silva.bacteria.rdp.tax)


ecoli

The ecoli parameter is used to provide a fasta file containing a single reference sequence (e.g. for e. coli) this must be aligned. mothur will trim to the start and end positions of the reference sequence.

start && end

The start parameter allows you to provide a starting position to trim to. The end parameter allows you to provide a ending position to trim from.

mothur > pcr.seqs(fasta=silva.bacteria.fasta, start=1044, end=6500)


nomatch

The nomatch parameter allows you to decide what to do with sequences where the primer is not found. Default=reject, meaning remove from fasta file. if nomatch=keep, then do nothing to sequence.

keepprimer

The keepprimer parameter allows you to keep the primer, default=false. Note that even if keepprimer=TRUE, mothur will not remove primers if it finds multiple potential matches to a given primer sequence.

keepdots

The keepdots parameter allows you to keep the leading and trailing .’s, default=true.

processors

The processors parameter allows you to use multiple processors. Default processors=Autodetect number of available processors and use all available.

Revisions

• 1.25.0 - First Introduced
• 1.28.0 - Added count option
• 1.29.2 - Bug Fix: If file was aligned and both forward and reverse primers were given, reverse primer was not trimmed properly.
• 1.30.0 - Added pdiffs parameter
• 1.30.2 - Bug Fix: For forward primer trimming with aligned sequences and keepdots=t. If the character before the first primer base was a base and not a gap the base was not trimmed. https://forum.mothur.org/viewtopic.php?f=4&t=2209