Effective Genome Size
A number of tools can accept an “effective genome size”. This is defined as the length of the “mappable” genome. There are two common alternative ways to calculate this:
1. The number of non-N bases in the genome.
2. The number of regions (of some size) in the genome that are uniquely mappable (possibly given some maximal edit distance).
Option 1 can be computed using faCount from Kents tools.
The effective genome size for a number of genomes using this method is given below:
Genome |
Effective size |
|---|---|
GRCh37 |
2864785220 |
GRCh38 |
3130250755 |
T2T/CHM13CAT_v2 |
3117275501 |
GRCm37 (mm9) |
1423050485 |
GRCm38 (mm10) |
2677331709 |
GRCm39 |
2654621783 |
dm3 |
162367812 |
dm6 |
142573024 |
GRCz10 |
1369631918 |
GRCz11 |
1368780147 |
WBcel235 |
100286401 |
TAIR10 |
119481543 |
These values only appropriate if multimapping reads are included. If they are excluded (or there’s any MAPQ filter applied),
then values derived from option 2 are more appropriate.
These are then based on the read length.
We can approximate these values for various read lengths using the khmer program program and unique-kmers.py in particular.
A table of effective genome sizes given a read length using this method is provided below:
Read length |
GRCh37 |
GRCh38 |
T2T/CHM13CAT_v2 |
GRCm37 (mm9) |
GRCm38 (mm10) |
GRCm39 |
dm3 |
dm6 |
GRCz10 |
GRCz11 |
WBcel235 |
TAIR10 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
50 |
2685511454 |
2708135512 |
2725240337 |
1476879138 |
2308840607 |
2309746861 |
130428510 |
125494459 |
1195445541 |
1197575653 |
95159402 |
114339094 |
75 |
2736124898 |
2765133945 |
2786136059 |
1547472870 |
2409741014 |
2410055689 |
135004387 |
127324557 |
1251132611 |
1250812288 |
96945370 |
115317469 |
100 |
2776919708 |
2820945813 |
2814334875 |
1639646400 |
2469638445 |
2468088461 |
139647132 |
129789773 |
1280188944 |
1280354977 |
98259898 |
118459858 |
150 |
2827436883 |
2884242992 |
2931551487 |
1750762288 |
2497610741 |
2495461690 |
144307658 |
129946309 |
1312207019 |
1311832909 |
98721103 |
118504138 |
200 |
2855463800 |
2909040054 |
2936403235 |
1856314701 |
2523838356 |
2521902382 |
148523810 |
132533877 |
1321355041 |
1322366338 |
98672558 |
117723393 |
250 |
2855044784 |
2919128283 |
2960856300 |
1924599658 |
2540234125 |
2538633971 |
151901455 |
132912060 |
1339205109 |
1342093482 |
101271756 |
119585546 |