Храните данные о чтении из SAM-отформатированного файла в объекте BioRead.
br =
BioRead with properties:
Quality: [1501x1 File indexed property]
Sequence: [1501x1 File indexed property]
Header: [1501x1 File indexed property]
NSeqs: 1501
Name: ''
Получите последовательности (чтения) из объекта.
Получите первые, третьи, и пятые последовательности из объекта.
seqs2 = 3x1 cell array
{'CACTAGTGGCTCATTGTAAATGTGTGGTTTAACTCG'}
{'AGTGGCTCATTGTAAATGTGTGGTTTAACTCGTCC' }
{'GCTCATTGTAAATGTGTGGTTTAACTCGTCCATGG' }
Получите первые пять положений тех последовательностей.
seqs3 = 3x1 cell array
{'CACTA'}
{'AGTGG'}
{'GCTCA'}
Можно использовать заголовок, чтобы получить соответствующие последовательности с тем заголовком. Если несколько последовательностей имеют тот же заголовок, функция возвращает все те последовательности.
Получите первые пять положений последовательностей с заголовком B7_591:4:96:693:509.
seqs4 = 1x1 cell array
{'CACTA'}
Получите первые, четвертые, и шестые положения первых трех последовательностей.
seq5 = 3x1 cell array
{'CTG'}
{'CGG'}
{'AGC'}