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)