www.gusucode.com > bioinfo 案例源码程序 matlab代码 > bioinfo/PlotPhredScoreDistributionsExample.m

    %% Plot Phred score distributions
% 
%%
% Create a |BioReadQualityStatistics| object from a FASTQ file using only the
% first 40 characters of each read with a minimum average quality score of
% 5.

% Copyright 2015 The MathWorks, Inc.

QSObj = BioReadQualityStatistics('SRR005164_1_50.fastq', ...
                                  'FilterLength', 40, ...
                                  'QualityScoreThreshold', 5);
%%
% Plot Phred quality scores at given base positions.
plotPerPositionQuality(QSObj)
%%
%