| Reads were mapped to the reference human genome build 36 |
| (hg18) using BLAT. Reads were realigned by cross_match to a local 40 |
| kb fragment of the reference genome, and all variation from each |
| read was saved. Variant positions were aggregated by genome |
| coordinate and filtered to remove false positives based on the |
| following criteria: Variant quality score, fraction of bases at each |
| position that were variant and proximity to a homopolymer run of 5 bp |
| or longer. See wheeler et al. Nature 2008 for details. |