exponenta event banner

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, a 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 которые выравниваются в пределах определенной области ссылочной последовательности.

Примеры

Построить 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) не отображаются в последовательностях или выравнивании.