seqRFLP (version 1.0.1)

clipprobe: Finding the sequences that could be clipped given two primers.

Description

Finding the sequences that could be clipped given the forward and backward primers.

Usage

clipprobe(fil, forProbe, bacProbe, tol = 3, clipped.only = TRUE)

Arguments

fil

Sequences in fasta format to be analyzed. Must be in "fasta" format.

forProbe

The forward primer.

bacProbe

The backward primer.

tol

Mumber of DNA base that can not be matched by the primers.

clipped.only

Whether only show the sequences that can be clipped.

Value

if clipped.only = TRUE, the function will return to the sequence(s) that could be obtained by PCR.

Details

This function can be used for predicting whether the sequences can be obtained by Polymerase chain reaction (PCR) given the forward and backward primers. Users may adjust the precision of DNA sites matching by the given primers on the sequences.

References

None.

See Also

See Also frag.dat for restriction enzyme clippling pattern.

Examples

### clipprobe() example

## Step1 Specify the forward and backward primer.
#clip the sequence between the backword and forward primer.
forProbe = ITS1F = 'CTTGGTCATTTAGAGGAAGTAA' # forward primer should be from the 5' to 3' end.
bacProbe = ITS4 = 'GCATATCAATAAGCGGAGGA'    # backward primer
#only sequence with two probes found could be clipped.

### Step2 reading sequences from external data in package.
directory <- system.file("extdata", package = "seqRFLP")
path <- file.path(directory, "seqs.fasta")
fas <- read.fasta(path)

## Step3 Clipping. Find clipped sequences, this usually takes less than 1 minute.
## please wait for a moment.
clipped <- clipprobe(fas, forProbe, bacProbe, tol = 2, clipped.only = TRUE)

## Step4 Checking the results.
## 368 selected sequences that could be clipped.
length(gnames.fas(clipped)) 
## 393 original sequences.
length(gnames.fas(fas))
## Sequences that can not be clipped.
setdiff(gnames.fas(fas), gnames.fas(clipped))