Plot summary statistics of BioRead object
[H,qsObj] = plotSummary(___)
plotSummary(BRObj) generates a summary statistics figure containing six plots about the average quality score for each base position, average quality score distribution, read count percentage for each base position, percentage of GC content for each base position, GC content distribution, and nucleotide distribution.
[H,qsObj] = plotSummary(___) returns a column vector H of handles to the axes in the generated figure and a BioReadQualityStatistics object qsObj using any of the input arguments from the previous syntaxes.
String specifying the format used for the characters encoding the sequence information and quality scores in a sequence file. Supported formats are 'Sanger', 'Illumina13', 'Illumina15', 'Illumina18', and 'Solexa'.
Positive integer (n) specifying the first n characters of each read to be used. Use an empty array to supress filtering.
Scalar value specifying a minimum average quality threshold for a read to be considered. Any read with an average score of less than the specified threshold value is ignored. The default value is –Inf, which causes all reads to be considered.
Column vector of handles to the axes in the generated figure.
BioReadQualityStatistics object that stores the data represented by the graphs.
Plot Summary Statistics of BioRead Object
Create a BioRead object from a FASTQ file.
BRObj = BioRead('SRR005164_1_50.fastq');
Plot summary statistics of the BioRead object BRObj using the first 40 characters of each read with the minimum average quality score of 5.