Graphs were generated using ggplot2 (64) and R base packages. Box plots of gene expression data show values corrected for confounding variance using gene expression PC as per eQTL analysis. Circular plots were created using the ggbio R package (65), and local association plots were generated with LocusZoom (66).
