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

    %% Plot Fractions of Reads Stratified by Phread Scores
% 

% Copyright 2015 The MathWorks, Inc.


%%
% 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.
QSObj = BioReadQualityStatistics('SRR005164_1_50.fastq','FilterLength',40,...
                                  'QualityScoreThreshold',5);
%%
% Plot the fractions of all reads with Phred scores in the ranges of 0-10,
% 11-20, 21-30, and 31-40 at each base position.
plotPerPositionCountByQuality(QSObj)