N50 statistic

In computational biology, the N50 statistic is a statistic of a set of contig or scaffold lengths. The N50 is similar to a mean or median, but has greater weight given to the longer contigs. It is used widely in genome assembly, especially in reference to contig lengths within a draft assembly. There are also the related N90, NG50, and D50 statistics.

Definition

N50

Given a set of contigs, each with its own length, the N50 length is defined as the length for which the collection of all contigs of that length or longer contains at least half of the sum of the lengths of all contigs, and for which the collection of all contigs of that length or shorter also contains at least half of the sum of the lengths of all contigs. (When more than one value of length meets both these criteria then the N50 is the average of the longest and shortest lengths that meet the criteria.) This can be thought of as the point of half of the mass of the distribution; the number of bases from all contigs shorter than the N50 will be close to the number of bases from all contigs longer than the N50.

N90

The N90 statistic is less than or equal to the N50 statistic; it is the length for which the collection of all contigs of that length or longer contains at least 90% of the sum of the lengths of all contigs, and for which the collection of all contigs of that length or shorter contains at least 10% of the sum of the lengths of all contigs.

NG50

Note that N50 is calculated in the context of the assembly size rather than the genome size. Therefore, comparisons of N50 values derived from assemblies of significantly different lengths are usually not informative, even if for the same genome. To address this, the authors of the Assemblathon competition derived a new measure called NG50. The NG50 statistic is the same as N50 except that it is 50% of the known or estimated genome size that must be achieved. This allows for meaningful comparisons between different assemblies.

D50

The D50 statistic (also termed D50 test) is similar to the N50 statistic in definition though it is generally not used to describe genome assemblies. The D50 statistic is the lowest value d for which the sum of the lengths of the largest d lengths is at least 50% of the sum of all of the lengths.

Examples

Consider two fictional, and highly simplified, genome assemblies (A & B) that are derived from two different species. Assembly A contains six contigs of lengths 80 kbp, 70 kbp, 50 kbp, 40 kbp, 30 kbp, and 20 kbp. The sum size of assembly A would be 290 kbp and so the N50 contig length would be 70 kbp (because 80 + 70 is greater than 50% of 290 kbp). Now lets assume that the contig lengths of assembly B are the same as assembly A except for the presence of two additional contigs of 10 kbp and 5 kbp. Thus the size of assembly B is 305 kbp, and now the N50 contig length drops to 50 kbp (80 + 70 + 50 is greater than 50% of 305 kbp). This example illustrates that one can sometimes increase the N50 length simply by removing some of the shortest contigs or scaffolds from an assembly.

If the estimated or known size of the genome from the fictional species A was 500 kbp, then the NG50 contig length would be 30 kbp (80 + 70 + 50 + 40 + 30 is greater than 50% of 500 kbp). In contrast, if the estimated or known size of the genome from species B was 350 kbp, then it would have an NG50 contig length of 50 kbp (80 + 70 + 50 is greater than 50% of 350 kbp).

Alternate computation

N50 can be found mathematically for a list L of positive integers as follows:

  1. Create another list L' , which is identical to L, except that every element n in L has been replaced with n copies of itself.
  2. The median of L' is the N50 of L. (The 10% quantile of L' is the N90 statistic.)

For example: If L = (2, 2, 2, 3, 3, 4, 8, 8), then L' consists of six 2's, six 3's, four 4's, and sixteen 8's. That is, L' has twice as many 2s as L; it has three times as many 3s as L; it has four times as many 4s; etc. The median of the 32-element set L' is the average of the 16th smallest element, 4, and 17th smallest element, 8, so the N50 is 6. We can see that the sum of all values in the list L that are smaller than or equal to the N50 of 6 is 16 = 2+2+2+3+3+4 and the sum of all values in the list L that are larger than or equal to 6 is also 16 = 8+8. For comparison with the N50 of 6, note that the mean of the list L is 4 while the median is 3.

Contradictory definitions

Some contradictions in the definition(s) of the N50 value have been identified, as discussed in a thread on the SEQ Answers forum. Also see

References

See also