plotPerSequenceQuality

Class: BioReadQualityStatistics

Plot distribution of average quality scores

Syntax

plotPerSequenceQuality(QSObj)
H = plotPerSequenceQuality(QSObj)

Description

plotPerSequenceQuality(QSObj) plots the distribution of average quality scores for short-read sequences of the BioReadQualityStatistics object QSObj.

H = plotPerSequenceQuality(QSObj) returns the handle H to the axes object containing the generated plot.

Input Arguments

QSObj

BioReadQualityStatistics object.

Output Arguments

H

Handle to axes object containing generated plot.

Examples

expand all

Plot the Distribution of Average Quality Scores

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);

Display the distribution of average quality scores.

plotPerSequenceQuality(QSObj)
ans = 

  Axes (Quality Distribution) with properties:

             XLim: [0 40]
             YLim: [0 50]
           XScale: 'linear'
           YScale: 'linear'
    GridLineStyle: '-'
         Position: [0.1300 0.1100 0.7750 0.8150]
            Units: 'normalized'

  Use GET to show all properties

Was this topic helpful?