Selecting the quality plot from the template display view menu displays a
horizontal plot as seen above. For each base in the consensus a quality is
computed based on the accuracy of data and strand coverage. This information
is then plotted in a colour and height specific to that type according to the
table below.
For example, in the figure we see that the first two hundred bases are mostly
only well determined on the reverse strand. Immediately following this is a
small section of double stranded sequence with many disagreements between
strands.
Note that when displaying many bases the screen resolution implies that the
quality codes for many bases will appear in the same screen pixel. However
the use of varying heights ensures that all problematic regions will be
visible, even when the problem is only with a single base position. Hence when
the quality plot consists of a single grey line all known quality problems
have been resolved, at the current consensus and quality cutoffs. As the
quality for a base position is determined from quality of coverage in each
strand and the agreement between the strands it is dependent on the consensus
and quality cutoffs.
The quality plot appears as "Calculate quality" in the results manager window.
Commands available on the data include "Information", which lists a summary of
the distribution of quality types to the output window, and "List" which lists
the actual quality values for each base to the output window. These quality
values are written in a textual form of single letters per base and are listed
below.
An example of the output using "Information" and "List" follows.
(Click for full size image)
grey 0 to 0 OK on both strands, both agree
blue 0 to 1 OK on plus strand only
green -1 to 0 OK on minus strand only
red -1 to 1 Bad on both strands
black -2 to 2 OK on both strands but they disagree
============================================================
Thu 11:09:29 AM: quality summary
------------------------------------------------------------
81.00 OK on both strands and they agree(a)
3.94 OK on plus strand only(b,d)
11.98 OK on minus strand only(c,e)
1.85 Bad on both strands(f,g,h,j)
1.22 OK on both strands but they disagree(i)
============================================================
Thu 11:09:35 AM: quality listing
------------------------------------------------------------
10 20 30 40 50 60
eeeeeeeeee eeeeeeeeee eeeeeeeeee eeeeeeehee eeeeeeeeee eeeeeeeeee
70 80 90 100 110 120
eeeeeeeeee eeeeeeeeee eeeeeeeeee eeeeeeeeee eeeeeeeeee eeeeeeeeee
130 140 150 160 170 180
eeeeeeeeee eeeeeeeeee eeeeeeeeee eeeeeeeeee eeeeeeeeee eeeeeeeeee
190 200 210 220 230 240
eeeeeeeeee eeeeeeeeee heeeeeeeee eeeeeeeici iiaiaciiia aaaaaaaaac
250 260 270 280 290 300
aaaacaaaaa aaaaaaaiia aaaaaaaaaa aaaaaaaaaa aaaabaaaaa aaaaaaaaaa
310 320 330 340 350 360
aaaaaaaaaa aaaaaaaaaa aaaaaaaaaa aaaaaaaaaa aaaaaaaaaa faaaaaaaaa
[ output removed for brevity ]
This page is maintained by
James Bonfield.
Last generated on 29 April 1996.
URL: http://www.mrc-lmb.cam.ac.uk/pubseq/manual/gap4_69.html