align2cigar

Преобразуйте выровненные последовательности в соответствующие подписи в формате CIGAR

Синтаксис

[Cigars,Starts] = align2cigar(Alignment,Ref)

Описание

[Cigars,Starts] = align2cigar(Alignment,Ref) преобразовывает выровненные последовательности, представленные в Alignment, массиве ячеек выровненных векторов символов, вектора строки или символьного массива, в Cigars, массив ячеек соответствующих отформатированных СИГАРОЙ векторов символов или вектора строки, с помощью ссылочной последовательности, заданной Ref, вектором символов или строкой. Это также возвращает Starts, вектор целых чисел, указывающих на положение запуска каждой выровненной последовательности относительно несодержащей разрывы ссылочной последовательности.

Входные параметры

Alignment

Массив ячеек вектора символов, представьте в виде строки вектор или символьный массив, представляющий выровненные последовательности. Мягкие усечения приняты, чтобы быть представленными строчными буквами в выровненных последовательностях. Пропущенные положения приняты, чтобы быть представленными. в выровненных последовательностях.

Ref

Вектор символов или строка, задающая выровненную ссылочную последовательность. Длина Ref должна равняться количеству столбцов в Alignment.

Выходные аргументы

Cigars

Массив ячеек отформатированных СИГАРОЙ векторов символов или вектора строки, соответствующего каждой выровненной последовательности в Alignment.

Starts

Вектор целых чисел, указывающих на положение запуска каждой выровненной последовательности относительно несодержащей разрывы ссылочной последовательности.

Примеры

свернуть все

Этот пример показывает, как преобразовать выровненные строки в строки CIGAR

Создайте массив ячеек выровненных строк, создайте строку, задающую ссылочную последовательность, и затем преобразуйте выравнивание в строки CIGAR:

aln = ['ACG-ATGC'; 'ACGT-TGC'; '  GTAT-C']
aln = 3x8 char array
    'ACG-ATGC'
    'ACGT-TGC'
    '  GTAT-C'

ref =  'ACGTATGC';
[cigar, start] = align2cigar(aln, ref)
cigar = 1x3 cell array
    {'3=1D4='}    {'4=1D3='}    {'4=1D1='}

start = 1×3

     1     1     3

Представленный в R2010b

Для просмотра документации необходимо авторизоваться на сайте