seqcomplement

Complement of nucleotide sequence

Syntax

SeqC = seqcomplement(SeqNT)

Description

SeqC = seqcomplement(SeqNT) returns the complementary sequence of the DNA or RNA nucleotide sequence contained in SeqNT. The format of SeqC matches SeqNT. For example, if SeqNT is a vector of integers, then so is SeqC. For complementary base pairs, see Complementary DNA and RNA Base Pairs. For complementary ambiguous characters, see Complementary IUPAC Ambiguity Codes.

Examples

Find Complement of DNA Nucleotide Sequence

Create a random DNA nucleotide sequence.

SeqNT = randseq(8,Alphabet="dna")

SeqNT = 
'TTATGACG'

Find the complement sequence.

SeqC = seqcomplement(SeqNT)

SeqC = 
'AATACTGC'

Input Arguments

SeqNT — Nucleotide sequence

character vector | string scalar | row vector of integers | structure

Nucleotide sequence, specified as one of the following:

• Character vector or string scalar containing the letters A, C, G, T, and U, and ambiguous characters R, Y, K, M, S, W, B, D, H, V, and N.

• Row vector of integers from Mapping Nucleotide Integers to Letter Codes.

• MATLAB^® structure containing a Sequence field that specifies a nucleotide sequence, such as the output returned by emblread, fastaread, fastqread, genbankread, getembl, or getgenbank.

More About

Complementary DNA and RNA Base Pairs

Nucleotide in Input Sequence	Nucleotide in Complimentary Sequence
A	T or U
C	G
G	C
T or U	A

Complementary IUPAC Ambiguity Codes

IUPAC Code in Input Sequence	IUPAC Code in Complimentary Sequence
M	K
R	Y
W	W
S	S
V	B
N	N

Version History

Introduced before R2006a