getAlignment

Класс: BioMap

Создайте выравнивание, представленную в BioMap объект

Синтаксис

Alignment = getAlignment(BioObj, StartPos, EndPos)
Alignment = getAlignment(BioObj, StartPos, EndPos, R)
Alignment = getAlignment(..., 'ParameterName', ParameterValue)
[Alignment, Indices] = getAlignment(...)

Описание

Alignment = getAlignment(BioObj, StartPos, EndPos) возвращает Alignment, символьный массив, содержащий выровненные последовательности чтения из BioObj, а BioMap объект. Последовательности чтения должны выровняться внутри определенной области ссылочной последовательности, которая определяется StartPos и EndPos, два положительных целого числа, таких что StartPos меньше EndPos, и оба меньше, чем длина ссылки последовательности.

Alignment = getAlignment(BioObj, StartPos, EndPos, R) выбирает ссылку где getAlignment восстанавливает выравнивание.

Alignment = getAlignment(..., 'ParameterName', ParameterValue) принимает одну или несколько пар имя/значение параметра, разделенных запятыми. Задайте ParameterName внутри одинарные кавычки.

[Alignment, Indices] = getAlignment(...) возвращает Indicesвектор индексов, задающих считанные последовательности, которые выравниваются внутри определенной области ссылки последовательности.

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

BioObj

Объект BioMap класс.

StartPos

Положительное целое число, которое задает начало области ссылки последовательности. StartPos должно быть меньше EndPosи меньше, чем общая длина последовательности ссылки.

EndPos

Положительное целое число, которое задает конец области ссылки последовательности. EndPos должно быть больше StartPosи меньше, чем общая длина последовательности ссылки.

R

Положительное целое число, индексирующее SequenceDictionary свойство BioObj, или вектор символов или строка, задающая фактическое имя ссылки.

Аргументы в виде пар имя-значение

'OffsetPad'

Определяет, добавляются ли заполняющие заготовки в начале каждой выровненной последовательности, чтобы представлять смещение начального положения каждой выровненной последовательности относительно ссылки. Варианты true или false (по умолчанию).

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

Alignment

Символьный массив, содержащий выровненные последовательности чтения из BioObj которые выравниваются внутри определенной области ссылки последовательности. Каждая строка символьного массива содержит одну выровненную последовательность, то есть положения последовательности, которые попадают в заданную область ссылки последовательности. Каждая выровненная последовательность может включать погрешности.

Indices

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

Примеры

Создайте a BioMap объект, а затем восстановите выравнивание между положениями 10 и 25 последовательности ссылки:

% Construct a BioMap object from a SAM file 
BMObj1 = BioMap('ex1.sam');
% Construct the alignment between positions 10 and 25 of the
% reference sequence. 
Alignment = getAlignment(BMObj1, 10, 25)
Alignment =

CTCATTGTAAATGTGT
CTCATTGTAAATGTGT
CTCATTGTAAATGTGT
CTCATTGTAATTTTTT
CTCATTGTAAATGTGT
   ATTGTAAATGTGT
   ATTGTAAATGTGT
     TGTAAATGTGT
        AAATGTGT
            GTGT
            GTGT
              GT

Алгоритмы

getAlignment предполагает, что опорная последовательность не имеет погрешностей. Поэтому положения в показаниях, соответствующих вставкам (I) и заполнению (P), не появляются в выравнивании.

Поскольку мягкие обрезанные положения (S) не связаны с положениями, которые выравниваются по ссылочной последовательности, они не появляются в выравнивании.

Пропущенное положение (N) появляется как . (период) в выравнивании.

Жесткие отсеченные положения (H) не появляются в последовательностях или выравнивании.