Main Content

ntdensity

Plot density of nucleotides along sequence

Syntax

ntdensity(SeqNT)
Density = ntdensity(SeqNT)
... = ntdensity(..., 'Window', WindowValue, ...)
[Density, HighCG] = ntdensity(..., 'CGThreshold', CGThresholdValue, ...)

Arguments

SeqNT

One of the following:

Note

Although you can submit a sequence with nucleotides other than A, C, G, and T, ntdensity plots only A, C, G, and T.

WindowValue

Value that specifies the window length for the density calculation. Default is length(SeqNT)/20.

CGThresholdValue

Controls the return of indices for regions where the CG content of SeqNT is greater than CGThresholdValue. Default is 5.

Description

ntdensity(SeqNT) plots the density of nucleotides A, C, G, and T in sequence SeqNT.

Density = ntdensity(SeqNT) returns a MATLAB structure with the density of nucleotides A, C, G, and T.

... = ntdensity(SeqNT, ...'PropertyName', PropertyValue, ...) calls ntdensity with optional properties that use property name/property value pairs. You can specify one or more properties in any order. Each PropertyName must be enclosed in single quotation marks and is case insensitive. These property name/property value pairs are as follows:

... = ntdensity(..., 'Window', WindowValue, ...) uses a window of length WindowValue for the density calculation. Default WindowValue is length(SeqNT)/20.

[Density, HighCG] = ntdensity(..., 'CGThreshold', CGThresholdValue, ...) returns indices for regions where the CG content of SeqNT is greater than CGThresholdValue. Default CGThresholdValue is 5.

Examples

  1. Create a random character vector to represent a nucleotide sequence.

    s = randseq(1000, 'alphabet', 'dna');
  2. Plot the density of nucleotides along the sequence.

    ntdensity(s)

Version History

Introduced before R2006a