seqdisp - Format long sequence output for easy viewing

Syntax

seqdisp(Seq)

seqdisp(Seq, ...'Row', RowValue, ...)
seqdisp(Seq, ...'Column', ColumnValue, ...)
seqdisp(Seq, ...'ShowNumbers', ShowNumbersValue, ...)

Arguments

Seq	Nucleotide or amino acid sequence represented by any of the following: Character array FASTA file name MATLAB structure with the field Sequence Multiply aligned sequences are allowed. FASTA files can have the file extension fa, fasta, fas, fsa, or fst.
RowValue	Integer that specifies the length of each row. Default is 60.
ColumnValue	Integer that specifies the column width or number of symbols before displaying a space. Default is 10.
ShowNumbersValue	Controls the display of numbers at the start of each row. Choices are true (default) to show numbers, or false to hide numbers.

Description

seqdisp(Seq) displays a sequence in rows, with a default row length of 60 and a default column width of 10.

seqdisp(Seq, ...'PropertyName', PropertyValue, ...) calls seqdisp 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:

seqdisp(Seq, ...'Row', RowValue, ...) specifies the length of each row for the displayed sequence.

seqdisp(Seq, ...'Column', ColumnValue, ...) specifies the number of letters to display before adding a space. RowValue must be larger than and evenly divisible by ColumnValue.

seqdisp(Seq, ...'ShowNumbers', ShowNumbersValue, ...) controls the display of numbers at the start of each row. Choices are true (default) to show numbers, or false to hide numbers.

Examples

Read sequence information from the GenBank database. Display the sequence in rows with 50 letters, and within a row, separate every 10 letters with a space.

mouseHEXA = getgenbank('AK080777');
seqdisp(mouseHEXA, 'Row', 50, 'Column', 10)

Create and save a FASTA file with two sequences, and then display it.

hdr = ['Sequence A'; 'Sequence B'];
seq = ['TAGCTGRCCAAGGCCAAGCGAGCTTN';'ATCGACYGGTTCCGGTTCGCTCGAAN']
fastawrite('local.fa', hdr, seq);
seqdisp('local.fa', 'ShowNumbers', false')

ans =
>Sequence A
 1  TAGCTGRCCA AGGCCAAGCG AGCTTN
>Sequence B
 1  ATCGACYGGT TCCGGTTCGC TCGAAN

See Also

getgenbank | multialignread | multialignwrite | seqconsensus | seqlogo | seqprofile | seqshoworfs | seqshowwords | seqtool

© 1984-2012- The MathWorks, Inc.    -   Site Help   -   Patents   -   Trademarks   -   Privacy Policy   -   Preventing Piracy   -   RSS