exponenta event banner

fastqwrite

Запись в файл в формате FASTQ

Синтаксис

fastqwrite(File, FASTQStruct)
fastqwrite(File, Header, Sequence, Qual)

Описание

fastqwrite(File, FASTQStruct) записывает содержимое структуры MATLAB ® или массива структур в файл в формате FASTQ. Если указан существующий файл в формате FASTQ ,fastqwrite добавляет данные в файл вместо перезаписи файла.

fastqwrite(File, Header, Sequence, Qual) записывает данные заголовка, последовательности и качества в файл в формате FASTQ.

Совет

Чтобы добавить данные в формате FASTQ к существующему файлу, просто укажите это имя файла. fastqwrite добавляет данные в конец файла.

Если используется fastqwrite в сценарии можно отключить предупреждающее сообщение добавления, введя следующие командные строки перед fastqwrite команда:

warnState = warning %Save the current warning state
warning('off','Bioinfo:fastqwrite:AppendToFile'); 
Затем введите следующую командную строку после fastqwrite команда:
warning(warnState) %Reset warning state to previous settings

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

File

Символьный вектор или строка, указывающая либо имя файла, либо путь и имя файла для сохранения данных в формате FASTQ. Если указано только имя файла, fastqwrite сохраняет файл в текущей папке MATLAB. Если указан существующий файл, fastqwrite добавляет данные в файл вместо перезаписи файла.

FASTQStruct

Структура MATLAB или массив структур, содержащих поля Header, Sequence, и Quality, например, возвращено fastqread.

Header

Символьный вектор или строка, содержащая информацию заголовка о нуклеотидной последовательности. Этот текст отображается в заголовке файла в формате FASTQ. File.

Sequence

Символьный вектор или строка, содержащая нуклеотидную последовательность с использованием стандартных букв или целочисленных кодов IUB/IUPAC. Список допустимых символов см. в разделе Аминокислотный или нуклеотидный поиск.

Qual

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

Примеры

Запись нескольких последовательностей в файл FASTQ из массива структур:

% Read the contents of a FASTQ-formatted file into
% an array of structures
reads = fastqread('SRR005164_1_50.fastq');
% Create another array of structures for the first five reads
reads5 = reads(1:5);
% Write the first five reads to a separate FASTQ-formatted file
fastqwrite('fiveReads.fastq', reads5)

Запишите одну последовательность в файл FASTQ из отдельных переменных:

% Create separate variables for the header, sequence, and 
% quality information of a nucleotide sequence
h = 'MYSEQ-000_1_1_1_953_493';
s = 'GTTACCATGATGTTATTTCTTCATTTGGAGGTAAAA';
q = ']]]]]]]]]]]]]]]]]]]]]]T]]]]RJRZTQLOA';
% Write the information to a FASTQ-formatted file
fastqwrite('oneRead.fastq', h, s, q)

Подробнее

свернуть все

Формат файла FASTQ

Файл в формате FASTQ содержит нуклеотидную последовательность и информацию о качестве по четырем линиям:

  • Строка 1 - Информация заголовка с префиксом @ символ

  • Линия 2 - Нуклеотидная последовательность

  • Строка 3 - Информация заголовка с префиксом + символ

  • Линия 4 - ASCII представление показателей качества по основаниям для нуклеотидной последовательности с использованием кодирования Phred или Solexa

Представлен в R2009b