| Bioinformatics Toolbox | |
Display open reading frames in a sequence
seqshoworfs(SeqNT, 'PropertyName', PropertyValue) seqshoworfs(..., 'Frames', FramesValue) seqshoworfs(..., 'GeneticCode', GeneticCodeValue) seqshoworfs(..., 'MinimumLength', MinimumLengthValue) seqshoworfs(..., 'AlternativeStartCodons', StartCodonsValue) seqshoworfs(..., 'Color', ColorValue) seqshoworfs(..., 'Columns', ColumnsValue)
SeqNT | Nucleotide sequence. Enter either a character string with the characters A, T (U), G, C, and ambiguous characters R, Y, K, M, S, W, B, D, H, V, N, or a vector of integers. You can also enter a structure with the field Sequence. |
FramesValue | Property to select the frame. Enter 1, 2, 3, -1, -2, -3, enter a vector with integers, or 'all'. The default value is the vector [1 2 3]. Frames -1, -2, and -3 correspond to the first, second, and third reading frames for the reverse complement. |
GeneticCodeValue | Genetic code name. Enter a code number or a code name from the table geneticcode. |
MinimumLengthValue | Property to set the minimum number of codons in an ORF. |
StartCodonsValue | Property to control using alternative start codons. Enter either true or false. The default value is false. |
ColorValue | Property to select the color for highlighting the reading frame. Enter either a 1-by-3 RGB vector specifying the intensity (0 to 255) of the red, green, and blue components of the color, or a character from the following list: 'b'—blue, 'g'—green, 'r'—red, 'c'—cyan, 'm'—magenta, or 'y'—yellow. To specify different colors for the three reading frames, use a 1-by-3 cell array of color values. If you are displaying reverse complement reading frames, then COLOR should be a 1-by-6 cell array of color values. |
ColumnsValue | Property to specify the number of columns in the output. |
seqshoworfs identifies and highlights all open reading frames using the standard or an alternative genetic code.
seqshoworfs(SeqNT) displays the sequence with all open reading frames highlighted, and it returns a structure of start and stop positions for each ORF in each reading frame. The standard genetic code is used with start codon 'AUG' and stop codons 'UAA', 'UAG', and 'UGA'.
seqshoworfs(..., 'Frames', FramesValue) specifies the reading frames to display. The default is to display the first, second, and third reading frames with ORFs highlighted in each frame.
seqshoworfs(..., 'GeneticCode', GeneticCodeValue) specifies the genetic code to use for finding open reading frames.
seqshoworfs(..., 'MinimumLength', MinimumLengthValue) sets the minimum number of codons for an ORF to be considered valid. The default value is 10.
seqshoworfs(..., 'AlternativeStartCodons', StartCodonsValue) uses alternative start codons if AlternativeStartCodons is set to true. For example, in the human mitochondrial genetic code, AUA and AUU are known to be alternative start codons. For more details of alternative start codons, see
http://www.ncbi.nlm.nih.gov/Taxonomy/Utils/ wprintgc.cgi?mode=t#SG1
seqshoworfs(..., 'Color', ColorValue) selects the color used to highlight the open reading frames in the output display. The default color scheme is blue for the first reading frame, red for the second, and green for the third frame.
seqshoworfs(..., 'Columns', ColumnsValue) specifies how many columns per line to use in the output. The default value is 64.
Look for the open reading frames in a random nucleotide sequence.
s = randseq(200,'alphabet', 'dna'); seqshoworfs(s);
Identify the open reading frames in a GenBank sequence.
HLA_DQB1 = getgenbank('NM_002123');
seqshoworfs(HLA_DQB1.Sequence);
Bioinformatics Toolbox functions codoncount, geneticcode, seqdisp,seqshowwords, seqwordcount
MATLAB function regexp
| | seqreverse | seqshowwords | |
© 1994-2005 The MathWorks, Inc.