Храните данные о чтении из SAM-отформатированного файла в объекте BioRead.
br =
BioRead with properties:
Quality: [1501x1 File indexed property]
Sequence: [1501x1 File indexed property]
Header: [1501x1 File indexed property]
NSeqs: 1501
Name: ''
Получите информацию о качестве последовательности.
Получите информацию о качестве последовательности из первых и третьих элементов в объекте.
seqQuals2 = 2x1 cell array
{'<<<<<<<<<<<<<<<;<<<<<<<<<5<<<<<;:<;7'}
{'<<<<<<<<<<<7;71<<;<;;<7;<<3;);3*8/5' }
Используйте логический вектор, чтобы получить ту же информацию.
seqQuals3 = 2x1 cell array
{'<<<<<<<<<<<<<<<;<<<<<<<<<5<<<<<;:<;7'}
{'<<<<<<<<<<<7;71<<;<;;<7;<<3;);3*8/5' }
Можно использовать заголовок, чтобы получить качество соответствующей последовательности с тем заголовком. Если несколько последовательностей имеют тот же заголовок, функция возвращает информацию о качестве всех тех последовательностей.
Получите качественную информацию последовательностей с заголовком B7_591:4:96:693:509.
seqQuals4 = 1x1 cell array
{'<<<<<<<<<<<<<<<;<<<<<<<<<5<<<<<;:<;7'}
Доступ к каждому свойству объекта с помощью записи через точку.
seqQuals2 = 2x1 cell array
{'<<<<<<<<<<<<<<<;<<<<<<<<<5<<<<<;:<;7'}
{'<<<<<<<<<<<7;71<<;<;;<7;<<3;);3*8/5' }