###Homework for Lab 2: DUE Friday, October 2nd
We have already discussed Hardy-Weinberg Equilibrium in class, but we haven’t really linked it to population genetics yet. Today’s lab will focus on the uses of Hardy-Weinberg Equilibrium in modern population genetics as we use R to look at Hardy-Weinberg Equilibrium in our assigned populations.
In modern population genetics studies, Hardy-Weinberg Equilibrium has a slightly different meaning than what we have talked about so far in class. In class, we learned that a population that is in true Hardy-Weinberg Equilibrium has to fit five conditions:
Obviously, holding a population of any species to this standard is pretty much impossible, so we will use Hardy-Weinberg differently to understand real populations (like humans) in our own population genetics study.
You will often see Hardy-Weinberg Equilibrium being used in modern population genetics to simply measure whether the observed genotype frequencies of a specific SNP or microsatellite are the same as the estimated genotype frequencies. Population geneticists will often use this measure of deviance as a “starting point,” so to speak, for their studies. They use Hardy-Weinberg as their starting point for a few reasons: to check the genotype quality of their data, to get a general sense of the amount of allelic variability and selection happening at specific loci, or, most commonly, as an assumption about the frequency of disease causing genotypes in their population. The assigned Benetti et al. (2020) and Maekawa et al. (2014) papers use Hardy-Weinberg in the latter manner, as a check to make sure there was no selection happening at their loci of interest.
In order to determine how much the observed genotype deviates from the expected genotype, we’ve been using a Chi-Squared Goodness-of-Fit test, the output of which is a P-value that will tell you how much the observed genotype frequency deviates from the expected frequency. However, we learned from the Wigginton reading that the Chi-Squared test has a tendency to give false positives, meaning it will tell us something is significantly different when it’s actually not; this is also called a Type I Error). In their paper, they propose a new mathematical method of calculating the “True” Hardy-Weinberg value, which fixes the problem of false positives. We will not need to get into the specifics of the equation to calculate the “True” Hardy-Weinberg Equilibrium, but you will be calculating your population’s “True” Hardy-Weinberg Equilibrium as well as using a traditional Chi-Squared test.
Now that we have learned a little more about the applications of Hardy-Weinberg Equilibrium in modern population genetics, let’s get to the module activities!
Use the SCC and R coding language to observe and understand population differences in ACE2 variation, as well as vcftools to understand differences in TMPRSS2 variation.
Calculate Hardy-Weinberg Equilibrium for all ACE2 OR TMPRSS2 SNPs in individual populations using a traditional Chi-Squared test.
Perform a check on all SNPs not in Hardy-Weinberg Equilibrium by calculating “True” Hardy-Weinberg with the built-in Shiny App.
Research the consequence types of these SNPs in order to better understand how these SNPs might affect the function of the gene itself, and also how they might affect subsequent phenotypes.
Calculate the “True” Hardy-Weinberg Equilibrium using the Shiny App for EITHER: ACE2 SNP 2074192, an exonic SNP that has known phenotypic consequences (according to Fan et al., 2011), and determine what the Hardy-Weinberg Equilibrium of that SNP says about the population.
OR
TMPRSS2 SNP rs12329760, an exonic SNP that has known phenotypic consequences (according to Maekawa et al., 2014), and determine what the Hardy-Weinberg Equilibrium of that SNP says about the population.
IMPORTANT NOTE: Some students will find that their ACE2 gene region will not work for the HWE exercises. This is due to a lack of variation data in their specific population (in other words, all individuals have the same allelic variants at all ACE2 loci). If you find that’s the case (e.g., your HWE table is empty!), then run the following exercises using TMPRSS2.
The first thing we need to do is access R Studio within the SCC environment. This is a little different from just logging in to the SCC. To use R Studio, we actually have to log in to the SCC On Demand interface:
Use this link to log in to SCC On Demand using your Kerberos password.
Once logged in, navigate to the Interactive Apps menu item.
Within that menu item, choose the RStudio Server
When that launches, your screen should look like this:
Click on the Launch
button.
It may take a moment, but once your SCC-based R Studio session is ready, you’ll see this on your screen:
Click on Connect to RStudio Server
. Once your R Studio interface launches, it will look like this:
The pane on the left is your console, which is where you will enter code to run. The pane on the upper right is your (environment), which is where anything running, or any datasets or items running in R will be shown. Finally, the pane on the lower right is your working directory. It should be your SCC home space when you first log in.
If you want to change your working directory to your anth333
project space, you may nagivate there by clicking on the ...
button on the upper right side of the pane and entering your project space name:
Finally, once you’ve navigated to the correct folder, you’ll want to set this as your working directory in R by clicking on the More
menu item in the pane, and cliking on Set As Working Directory
:
Conveniently enough, you can also access your space SCC via a Terminal window in this R Studio session using the Terminal tab next to your console. Please note that you may need to set your proper working directory in the Terminal, as well, using the cd
command (e.g., cd /project/anth333/caschmit
).
Now that we are in R Studio and have our data in place in our working directory (from Lab 1), we can get started. The first thing we need to do here is install two packages.
In R space, packages are essentially bundles of functions that R users have written to perform analyses that do not come with the base program. In this class we will be using packages that were designed specifically for population genetics analysis. To install these packages, run these two lines of code in the console (this may take a few minutes!):
install.packages("vcfR")
install.packages("pegas")
We only need to install these packages once, but we need to load them in to our workspace every time we open the program. To do this, use the following code:
library(vcfR)
library(pegas)
Now that we have all our packages loaded, we can get to the analysis! In the “files” tab in the bottom left quadrant, you should see the VCF files you downloaded in your folder (give you put them there properly).
Try these exercises first with your ACE2 VCF file.
For the YRI population, however, there is no variation in the ACE2 gene region. Because of this, the demonstration data here will have to be done with TMPRSS2. If you find that you also have no variation within ACE2 in your population, please also use TMPRSS2 for the remainder of the module.
We will need to load our VCF file into our R workspace much like we did with our packages. Unlike the packages, however, we need to give a name to the file. To make things simple, it would be a good idea to name this file in accordance with the acronym of your population. I will use the YRI population as an example. Use one of the two options below to enter your VCF data into the R space:
#If your file is zipped (e.g., ends in '.vcf.gz') then use this code:
TMPRSS2_YRI <- read.vcfR("TMPRSS2_YRI.vcf.gz", verbose = TRUE)
## Scanning file to determine attributes.
## File attributes:
## meta lines: 131
## header_line: 132
## variant count: 2246
## column count: 117
##
Meta line 131 read in.
## All meta lines processed.
## gt matrix initialized.
## Character matrix gt created.
## Character matrix gt rows: 2246
## Character matrix gt cols: 117
## skip: 0
## nrows: 2246
## row_num: 0
##
Processed variant 1000
Processed variant 2000
Processed variant: 2246
## All variants processed
#If your file is unzipped (e.g., ends in just '.vcf') then use this code:
TMPRSS2_YRI <- read.vcfR("TMPRSS2_YRI.vcf", verbose = TRUE)
This will give you a file in your workspace called a vcfR file. Most genetics-based packages have their own way of formatting and working with genetic data, and this is a format R can use. The next line of code will seem a little counterintuitive, because we have to change the format of the file again to work with functions in the pegas package. Doing this will not delete the original file (as long as you name the file something different), so we can use the vcfR file again in later analyses. The type of file we have to convert to for pegas is called a genind format, so we will do that here:
pegas.TMPRSS2_YRI <- vcfR2genind(TMPRSS2_YRI, sep = "[|/]")
Finally, now that we have our data in the format that we need, we can do our Hardy-Weinberg tests. This function will test whether all of your population’s SNPs in the coding region of ACE2 or TMPRSS2 are in Hardy-Weinberg Equilibrium. As an example, I will run the funtion with my pegas.TMPRSS2_YRI file. Remember to save your results by assigning them a name:
HWE <- hw.test(pegas.TMPRSS2_YRI, B = 0) #B MUST equal 0 here, changing this value will affect the result!
HWE #Running the name of your results file will print the results
## chi^2 df Pr(chi^2 >)
## rs575760704 0.000000e+00 0 1.000000e+00
## rs544440280 2.336398e-03 1 9.614482e-01
## rs564273360 0.000000e+00 0 1.000000e+00
## rs559811756 9.433138e-03 1 9.226276e-01
## rs553400003 0.000000e+00 0 1.000000e+00
## rs560183935 0.000000e+00 0 1.000000e+00
## rs181162134 0.000000e+00 0 1.000000e+00
## rs462471 2.466520e+00 1 1.162952e-01
## rs76000363 1.494810e+00 1 2.214717e-01
## rs143680939 1.494810e+00 1 2.214717e-01
## rs456142 2.466520e+00 1 1.162952e-01
## rs551045243 0.000000e+00 0 1.000000e+00
## rs571050655 0.000000e+00 0 1.000000e+00
## rs539639766 0.000000e+00 0 1.000000e+00
## rs560898303 0.000000e+00 0 1.000000e+00
## rs185078457 0.000000e+00 0 1.000000e+00
## rs112657409 7.881619e-01 1 3.746565e-01
## rs2838038 2.545009e-01 1 6.139235e-01
## rs555769183 0.000000e+00 0 1.000000e+00
## rs462574 3.266541e+00 1 7.070658e-02
## rs456298 2.466520e+00 1 1.162952e-01
## rs564403172 0.000000e+00 0 1.000000e+00
## rs142444476 2.336398e-03 1 9.614482e-01
## rs540046047 0.000000e+00 0 1.000000e+00
## rs560408119 0.000000e+00 0 1.000000e+00
## rs117652812 0.000000e+00 0 1.000000e+00
## rs548592256 0.000000e+00 0 1.000000e+00
## rs140532244 0.000000e+00 0 1.000000e+00
## rs531183955 0.000000e+00 0 1.000000e+00
## rs551106657 0.000000e+00 0 1.000000e+00
## rs17001042 5.439334e-01 1 4.608072e-01
## rs11910678 2.809573e+00 1 9.370329e-02
## rs189546747 0.000000e+00 0 1.000000e+00
## rs567175997 0.000000e+00 0 1.000000e+00
## rs183418223 9.433138e-03 1 9.226276e-01
## rs556079659 0.000000e+00 0 1.000000e+00
## rs569303172 0.000000e+00 0 1.000000e+00
## rs538240292 0.000000e+00 0 1.000000e+00
## rs558272207 0.000000e+00 0 1.000000e+00
## rs578005531 0.000000e+00 0 1.000000e+00
## rs187742976 0.000000e+00 0 1.000000e+00
## rs377348479 0.000000e+00 0 1.000000e+00
## rs573615979 2.142432e-02 1 8.836289e-01
## rs542538936 0.000000e+00 0 1.000000e+00
## rs75317729 9.433138e-03 1 9.226276e-01
## rs77675406 2.545009e-01 1 6.139235e-01
## rs544635611 0.000000e+00 0 1.000000e+00
## rs542458473 0.000000e+00 0 1.000000e+00
## rs12627374 0.000000e+00 0 1.000000e+00
## rs547186777 0.000000e+00 0 1.000000e+00
## rs567103273 0.000000e+00 0 1.000000e+00
## rs372405355 0.000000e+00 0 1.000000e+00
## rs549745562 0.000000e+00 0 1.000000e+00
## rs569561312 0.000000e+00 0 1.000000e+00
## rs538252674 0.000000e+00 0 1.000000e+00
## rs557822084 0.000000e+00 0 1.000000e+00
## rs191697229 2.336398e-03 1 9.614482e-01
## rs150014829 0.000000e+00 0 1.000000e+00
## rs62217525 0.000000e+00 0 1.000000e+00
## rs146564124 0.000000e+00 0 1.000000e+00
## rs542402063 0.000000e+00 0 1.000000e+00
## rs77996454 2.545009e-01 1 6.139235e-01
## rs149695119 5.883901e-02 1 8.083404e-01
## rs141230106 0.000000e+00 0 1.000000e+00
## rs200164183 0.000000e+00 0 1.000000e+00
## rs564691729 0.000000e+00 0 1.000000e+00
## rs527261679 0.000000e+00 0 1.000000e+00
## rs370347248 0.000000e+00 0 1.000000e+00
## rs374666536 0.000000e+00 0 1.000000e+00
## rs142194573 0.000000e+00 0 1.000000e+00
## rs118028230 0.000000e+00 0 1.000000e+00
## rs73905370 2.545009e-01 1 6.139235e-01
## rs569577322 0.000000e+00 0 1.000000e+00
## rs569583276 2.336398e-03 1 9.614482e-01
## rs55896064 2.545009e-01 1 6.139235e-01
## rs571587918 0.000000e+00 0 1.000000e+00
## rs534192391 0.000000e+00 0 1.000000e+00
## rs547523412 0.000000e+00 0 1.000000e+00
## rs57474639 3.737024e-01 1 5.409930e-01
## rs79971314 6.629730e-01 1 4.155123e-01
## rs555911696 9.433138e-03 1 9.226276e-01
## rs73372161 3.419130e-03 1 9.533716e-01
## rs538588854 0.000000e+00 0 1.000000e+00
## rs558791980 0.000000e+00 0 1.000000e+00
## rs571993620 0.000000e+00 0 1.000000e+00
## rs61459778 7.626580e-03 1 9.304090e-01
## rs145024812 9.433138e-03 1 9.226276e-01
## rs574538362 0.000000e+00 0 1.000000e+00
## rs138812485 0.000000e+00 0 1.000000e+00
## rs377591443 0.000000e+00 0 1.000000e+00
## rs531920590 0.000000e+00 0 1.000000e+00
## rs551817968 0.000000e+00 0 1.000000e+00
## rs566012796 0.000000e+00 0 1.000000e+00
## rs527693716 0.000000e+00 0 1.000000e+00
## rs148910500 0.000000e+00 0 1.000000e+00
## rs115975538 1.004612e+00 1 3.161970e-01
## rs536311940 0.000000e+00 0 1.000000e+00
## rs118133613 0.000000e+00 0 1.000000e+00
## rs143595083 9.433138e-03 1 9.226276e-01
## rs148038688 9.433138e-03 1 9.226276e-01
## rs546913068 0.000000e+00 0 1.000000e+00
## rs572265353 0.000000e+00 0 1.000000e+00
## rs534326511 0.000000e+00 0 1.000000e+00
## rs73905371 3.737024e-01 1 5.409930e-01
## rs73372163 1.088502e-02 1 9.169064e-01
## rs543578752 0.000000e+00 0 1.000000e+00
## rs563352095 0.000000e+00 0 1.000000e+00
## rs73372166 1.088502e-02 1 9.169064e-01
## rs113562865 1.887624e-01 1 6.639487e-01
## rs455281 7.305144e-01 1 3.927165e-01
## rs457909 0.000000e+00 0 1.000000e+00
## rs547401659 0.000000e+00 0 1.000000e+00
## rs461194 0.000000e+00 0 1.000000e+00
## rs73372168 1.211511e-01 1 7.277894e-01
## rs188235654 0.000000e+00 0 1.000000e+00
## rs570126729 0.000000e+00 0 1.000000e+00
## rs532216261 0.000000e+00 0 1.000000e+00
## rs552303675 0.000000e+00 0 1.000000e+00
## rs548979551 0.000000e+00 0 1.000000e+00
## rs192955773 9.433138e-03 1 9.226276e-01
## rs554890599 0.000000e+00 0 1.000000e+00
## rs540024706 0.000000e+00 0 1.000000e+00
## rs138765307 0.000000e+00 0 1.000000e+00
## rs556787864 0.000000e+00 0 1.000000e+00
## rs200395836 0.000000e+00 0 1.000000e+00
## rs183984610 0.000000e+00 0 1.000000e+00
## rs142750000 0.000000e+00 0 1.000000e+00
## rs572530227 0.000000e+00 0 1.000000e+00
## rs541351488 0.000000e+00 0 1.000000e+00
## rs374377315 0.000000e+00 0 1.000000e+00
## rs147359020 0.000000e+00 0 1.000000e+00
## rs186826600 0.000000e+00 0 1.000000e+00
## rs563633617 0.000000e+00 0 1.000000e+00
## rs532610867 0.000000e+00 0 1.000000e+00
## rs552105586 0.000000e+00 0 1.000000e+00
## rs565903681 0.000000e+00 0 1.000000e+00
## rs139015396 9.433138e-03 1 9.226276e-01
## rs548552823 0.000000e+00 0 1.000000e+00
## rs149855493 9.433138e-03 1 9.226276e-01
## rs546295824 0.000000e+00 0 1.000000e+00
## rs144823388 9.433138e-03 1 9.226276e-01
## rs570467504 0.000000e+00 0 1.000000e+00
## rs191394761 0.000000e+00 0 1.000000e+00
## rs148561341 0.000000e+00 0 1.000000e+00
## rs572351113 0.000000e+00 0 1.000000e+00
## rs76315847 3.844785e-02 1 8.445466e-01
## rs554852594 2.336398e-03 1 9.614482e-01
## rs185018051 0.000000e+00 0 1.000000e+00
## rs144359794 0.000000e+00 0 1.000000e+00
## rs146605032 0.000000e+00 0 1.000000e+00
## rs743542 5.216761e-04 1 9.817777e-01
## rs190618812 0.000000e+00 0 1.000000e+00
## rs181872728 0.000000e+00 0 1.000000e+00
## rs528337775 0.000000e+00 0 1.000000e+00
## rs548354821 0.000000e+00 0 1.000000e+00
## rs561982326 0.000000e+00 0 1.000000e+00
## rs113438095 0.000000e+00 0 1.000000e+00
## rs148988435 0.000000e+00 0 1.000000e+00
## rs570454392 0.000000e+00 0 1.000000e+00
## rs371012018 0.000000e+00 0 1.000000e+00
## rs185312677 2.336398e-03 1 9.614482e-01
## rs73372170 3.324100e-02 1 8.553307e-01
## rs534892317 0.000000e+00 0 1.000000e+00
## rs462321 3.803969e+00 1 5.113126e-02
## rs462326 2.989617e+00 1 8.380000e-02
## rs141323355 0.000000e+00 0 1.000000e+00
## rs143672898 0.000000e+00 0 1.000000e+00
## rs577217013 0.000000e+00 0 1.000000e+00
## rs545870269 0.000000e+00 0 1.000000e+00
## rs146797606 9.433138e-03 1 9.226276e-01
## rs573343651 0.000000e+00 0 1.000000e+00
## rs542060428 0.000000e+00 0 1.000000e+00
## rs565902164 0.000000e+00 0 1.000000e+00
## rs140428704 0.000000e+00 0 1.000000e+00
## rs530701233 0.000000e+00 0 1.000000e+00
## rs7283324 2.729747e+00 1 9.849486e-02
## rs462448 0.000000e+00 0 1.000000e+00
## rs532950488 0.000000e+00 0 1.000000e+00
## rs546232086 0.000000e+00 0 1.000000e+00
## rs575517837 0.000000e+00 0 1.000000e+00
## rs189793246 0.000000e+00 0 1.000000e+00
## rs548653178 0.000000e+00 0 1.000000e+00
## rs568346732 0.000000e+00 0 1.000000e+00
## rs150382508 9.433138e-03 1 9.226276e-01
## rs115774576 9.433138e-03 1 9.226276e-01
## rs376074355 0.000000e+00 0 1.000000e+00
## rs114549926 9.433138e-03 1 9.226276e-01
## rs181461974 2.336398e-03 1 9.614482e-01
## rs376403654 0.000000e+00 0 1.000000e+00
## rs542368768 0.000000e+00 0 1.000000e+00
## rs138064996 0.000000e+00 0 1.000000e+00
## rs7279603 1.185306e+00 1 2.762784e-01
## rs544036587 0.000000e+00 0 1.000000e+00
## rs67617179 4.089941e-01 1 5.224805e-01
## rs7278739 4.295801e+00 1 3.820660e-02
## rs533050705 0.000000e+00 0 1.000000e+00
## rs539858204 0.000000e+00 0 1.000000e+00
## rs7278627 1.887624e-01 1 6.639487e-01
## rs149109132 9.433138e-03 1 9.226276e-01
## rs548718178 0.000000e+00 0 1.000000e+00
## rs568605816 0.000000e+00 0 1.000000e+00
## rs531148133 0.000000e+00 0 1.000000e+00
## rs185416185 0.000000e+00 0 1.000000e+00
## rs571137450 0.000000e+00 0 1.000000e+00
## rs539777782 0.000000e+00 0 1.000000e+00
## rs146052428 2.142432e-02 1 8.836289e-01
## rs566918037 0.000000e+00 0 1.000000e+00
## rs370663590 0.000000e+00 0 1.000000e+00
## rs2070793 1.185306e+00 1 2.762784e-01
## rs117898838 2.336398e-03 1 9.614482e-01
## rs544417878 0.000000e+00 0 1.000000e+00
## rs2070792 1.185306e+00 1 2.762784e-01
## rs577753477 0.000000e+00 0 1.000000e+00
## rs540165785 0.000000e+00 0 1.000000e+00
## rs2070791 2.336398e-03 1 9.614482e-01
## rs2070790 4.047374e+00 1 4.424011e-02
## rs142769034 0.000000e+00 0 1.000000e+00
## rs147368850 0.000000e+00 0 1.000000e+00
## rs111572592 0.000000e+00 0 1.000000e+00
## rs2070789 2.729747e+00 1 9.849486e-02
## rs147986574 0.000000e+00 0 1.000000e+00
## rs141603473 1.211511e-01 1 7.277894e-01
## rs547356228 0.000000e+00 0 1.000000e+00
## rs542946711 0.000000e+00 0 1.000000e+00
## rs113564116 0.000000e+00 0 1.000000e+00
## rs555771010 0.000000e+00 0 1.000000e+00
## rs2070788 2.325824e+00 1 1.272428e-01
## rs537989379 0.000000e+00 0 1.000000e+00
## rs150907799 0.000000e+00 0 1.000000e+00
## rs139432971 2.336398e-03 1 9.614482e-01
## rs189832305 0.000000e+00 0 1.000000e+00
## rs9974589 1.454682e+00 1 2.277781e-01
## rs2070787 4.047374e+00 1 4.424011e-02
## rs542228195 0.000000e+00 0 1.000000e+00
## rs2070786 4.227761e+00 1 3.976788e-02
## rs575968857 0.000000e+00 0 1.000000e+00
## rs544822322 0.000000e+00 0 1.000000e+00
## rs564981706 0.000000e+00 0 1.000000e+00
## rs143060022 0.000000e+00 0 1.000000e+00
## rs547168890 0.000000e+00 0 1.000000e+00
## rs560788815 0.000000e+00 0 1.000000e+00
## rs369619813 0.000000e+00 0 1.000000e+00
## rs78503214 8.816327e-02 1 7.665253e-01
## rs569186960 0.000000e+00 0 1.000000e+00
## rs147233451 0.000000e+00 0 1.000000e+00
## rs146132415 0.000000e+00 0 1.000000e+00
## rs61735794 0.000000e+00 0 1.000000e+00
## rs140121827 0.000000e+00 0 1.000000e+00
## rs143712818 9.433138e-03 1 9.226276e-01
## rs112132031 3.419130e-03 1 9.533716e-01
## rs573482152 0.000000e+00 0 1.000000e+00
## rs180831387 0.000000e+00 0 1.000000e+00
## rs28524972 4.535967e+00 1 3.318970e-02
## rs455922 0.000000e+00 0 1.000000e+00
## rs544946928 0.000000e+00 0 1.000000e+00
## rs465576 9.461447e-02 1 7.583908e-01
## rs572143040 0.000000e+00 0 1.000000e+00
## rs540844961 0.000000e+00 0 1.000000e+00
## rs456016 2.336398e-03 1 9.614482e-01
## rs186429734 0.000000e+00 0 1.000000e+00
## rs549459413 0.000000e+00 0 1.000000e+00
## rs563188125 0.000000e+00 0 1.000000e+00
## rs9305744 2.989617e+00 1 8.380000e-02
## rs35050484 0.000000e+00 0 1.000000e+00
## rs192334131 0.000000e+00 0 1.000000e+00
## rs148509204 0.000000e+00 0 1.000000e+00
## rs372665499 0.000000e+00 0 1.000000e+00
## rs568737536 0.000000e+00 0 1.000000e+00
## rs536026892 0.000000e+00 0 1.000000e+00
## rs556192760 0.000000e+00 0 1.000000e+00
## rs115207130 1.597633e-01 1 6.893745e-01
## rs538276300 0.000000e+00 0 1.000000e+00
## rs558492934 0.000000e+00 0 1.000000e+00
## rs142834250 0.000000e+00 0 1.000000e+00
## rs541155775 0.000000e+00 0 1.000000e+00
## rs554702767 0.000000e+00 0 1.000000e+00
## rs183650725 0.000000e+00 0 1.000000e+00
## rs543168691 0.000000e+00 0 1.000000e+00
## rs370896169 0.000000e+00 0 1.000000e+00
## rs532005855 0.000000e+00 0 1.000000e+00
## rs151035593 1.211511e-01 1 7.277894e-01
## rs34624090 1.388883e+00 1 2.385938e-01
## rs187052796 9.433138e-03 1 9.226276e-01
## rs467375 1.252544e+00 1 2.630671e-01
## rs142303004 2.142432e-02 1 8.836289e-01
## rs567380110 0.000000e+00 0 1.000000e+00
## rs146116431 0.000000e+00 0 1.000000e+00
## rs551100726 0.000000e+00 0 1.000000e+00
## rs569866330 0.000000e+00 0 1.000000e+00
## rs116511699 3.844785e-02 1 8.445466e-01
## rs191441931 0.000000e+00 0 1.000000e+00
## rs376752614 0.000000e+00 0 1.000000e+00
## rs534591508 0.000000e+00 0 1.000000e+00
## rs182796345 0.000000e+00 0 1.000000e+00
## rs370927492 0.000000e+00 0 1.000000e+00
## rs542865477 0.000000e+00 0 1.000000e+00
## rs139144487 0.000000e+00 0 1.000000e+00
## rs147827602 0.000000e+00 0 1.000000e+00
## rs545726689 0.000000e+00 0 1.000000e+00
## rs565237319 0.000000e+00 0 1.000000e+00
## rs527684898 0.000000e+00 0 1.000000e+00
## rs540987630 0.000000e+00 0 1.000000e+00
## rs561063944 0.000000e+00 0 1.000000e+00
## rs187078345 9.433138e-03 1 9.226276e-01
## rs533556786 0.000000e+00 0 1.000000e+00
## rs464431 2.336398e-03 1 9.614482e-01
## rs458213 1.533786e+00 1 2.155451e-01
## rs552550910 0.000000e+00 0 1.000000e+00
## rs565933739 2.336398e-03 1 9.614482e-01
## rs141301979 3.844785e-02 1 8.445466e-01
## rs458280 0.000000e+00 0 1.000000e+00
## rs7364083 2.597386e+00 1 1.070401e-01
## rs111671182 2.336398e-03 1 9.614482e-01
## rs191600263 0.000000e+00 0 1.000000e+00
## rs543253134 0.000000e+00 0 1.000000e+00
## rs538803792 0.000000e+00 0 1.000000e+00
## rs7364088 1.252544e+00 1 2.630671e-01
## rs182254633 9.433138e-03 1 9.226276e-01
## rs146957681 0.000000e+00 0 1.000000e+00
## rs187452528 0.000000e+00 0 1.000000e+00
## rs76135088 1.124531e+00 1 2.889448e-01
## rs117562633 0.000000e+00 0 1.000000e+00
## rs28548447 1.486659e-01 1 6.998135e-01
## rs114844880 3.844785e-02 1 8.445466e-01
## rs193253965 0.000000e+00 0 1.000000e+00
## rs559699558 0.000000e+00 0 1.000000e+00
## rs137962614 0.000000e+00 0 1.000000e+00
## rs79065748 5.187727e-01 1 4.713657e-01
## rs190385097 0.000000e+00 0 1.000000e+00
## rs536642840 0.000000e+00 0 1.000000e+00
## rs192011400 0.000000e+00 0 1.000000e+00
## rs570413089 0.000000e+00 0 1.000000e+00
## rs538864203 0.000000e+00 0 1.000000e+00
## rs116606027 9.433138e-03 1 9.226276e-01
## rs117333005 0.000000e+00 0 1.000000e+00
## rs535019926 0.000000e+00 0 1.000000e+00
## rs112980967 8.816327e-02 1 7.665253e-01
## rs184164037 9.433138e-03 1 9.226276e-01
## rs543404937 0.000000e+00 0 1.000000e+00
## rs563812695 0.000000e+00 0 1.000000e+00
## rs569259137 0.000000e+00 0 1.000000e+00
## rs73372182 9.367695e-02 1 7.595537e-01
## rs143097165 9.433138e-03 1 9.226276e-01
## rs528276471 0.000000e+00 0 1.000000e+00
## rs548552862 2.336398e-03 1 9.614482e-01
## rs561789442 0.000000e+00 0 1.000000e+00
## rs530540353 0.000000e+00 0 1.000000e+00
## rs146142989 0.000000e+00 0 1.000000e+00
## rs117827300 0.000000e+00 0 1.000000e+00
## rs188514624 9.433138e-03 1 9.226276e-01
## rs546273967 0.000000e+00 0 1.000000e+00
## rs565973484 0.000000e+00 0 1.000000e+00
## rs534885992 0.000000e+00 0 1.000000e+00
## rs73230068 0.000000e+00 0 1.000000e+00
## rs568619174 0.000000e+00 0 1.000000e+00
## rs181058683 0.000000e+00 0 1.000000e+00
## rs75168613 4.829224e-01 1 4.871016e-01
## rs577285910 0.000000e+00 0 1.000000e+00
## rs113288437 1.942675e-03 1 9.648440e-01
## rs140141551 0.000000e+00 0 1.000000e+00
## rs201984814 0.000000e+00 0 1.000000e+00
## rs143291395 0.000000e+00 0 1.000000e+00
## rs562033822 0.000000e+00 0 1.000000e+00
## rs530689404 0.000000e+00 0 1.000000e+00
## rs544474510 0.000000e+00 0 1.000000e+00
## rs151338439 0.000000e+00 0 1.000000e+00
## rs532756922 0.000000e+00 0 1.000000e+00
## rs546335233 0.000000e+00 0 1.000000e+00
## rs2298658 0.000000e+00 0 1.000000e+00
## rs150445636 0.000000e+00 0 1.000000e+00
## rs145355824 0.000000e+00 0 1.000000e+00
## rs2298659 3.507074e-01 1 5.537130e-01
## rs537370123 0.000000e+00 0 1.000000e+00
## rs17854725 7.550251e-01 1 3.848899e-01
## rs368812287 0.000000e+00 0 1.000000e+00
## rs199704563 0.000000e+00 0 1.000000e+00
## rs553516093 0.000000e+00 0 1.000000e+00
## rs377496737 0.000000e+00 0 1.000000e+00
## rs542031508 0.000000e+00 0 1.000000e+00
## rs555302570 0.000000e+00 0 1.000000e+00
## rs575981127 0.000000e+00 0 1.000000e+00
## rs113928389 9.433138e-03 1 9.226276e-01
## rs544221548 0.000000e+00 0 1.000000e+00
## rs2298660 3.078428e+00 1 7.933717e-02
## rs55964536 1.687500e+00 1 1.939309e-01
## rs2298661 2.620787e-03 1 9.591712e-01
## rs559830930 0.000000e+00 0 1.000000e+00
## rs77511690 5.187727e-01 1 4.713657e-01
## rs549038772 0.000000e+00 0 1.000000e+00
## rs562177468 0.000000e+00 0 1.000000e+00
## rs531121505 0.000000e+00 0 1.000000e+00
## rs368994585 0.000000e+00 0 1.000000e+00
## rs371053759 0.000000e+00 0 1.000000e+00
## rs144192191 5.169581e+00 1 2.298575e-02
## rs559550182 2.142432e-02 1 8.836289e-01
## rs539860607 0.000000e+00 0 1.000000e+00
## rs386519 2.336398e-03 1 9.614482e-01
## rs150875127 4.630485e-02 1 8.296226e-01
## rs535611212 9.433138e-03 1 9.226276e-01
## rs555820997 2.336398e-03 1 9.614482e-01
## rs111220492 9.367695e-02 1 7.595537e-01
## rs537997349 2.142432e-02 1 8.836289e-01
## rs557639119 2.142432e-02 1 8.836289e-01
## rs548463596 3.161565e-02 1 8.588738e-01
## rs577480005 2.336398e-03 1 9.614482e-01
## rs539864066 9.433138e-03 1 9.226276e-01
## rs559741904 9.433138e-03 1 9.226276e-01
## rs573564326 9.433138e-03 1 9.226276e-01
## rs542395007 0.000000e+00 0 1.000000e+00
## rs562598515 0.000000e+00 0 1.000000e+00
## rs531024324 0.000000e+00 0 1.000000e+00
## rs530918384 0.000000e+00 0 1.000000e+00
## rs551164492 0.000000e+00 0 1.000000e+00
## rs564793560 9.433138e-03 1 9.226276e-01
## rs61728255 0.000000e+00 0 1.000000e+00
## rs141400887 4.628476e+00 1 3.144548e-02
## rs111220509 2.370370e+00 1 1.236577e-01
## rs111220533 1.273219e-01 1 7.212252e-01
## rs114837856 9.854388e-02 1 7.535841e-01
## rs549212787 0.000000e+00 0 1.000000e+00
## rs569292495 0.000000e+00 0 1.000000e+00
## rs115720411 9.367695e-02 1 7.595537e-01
## rs557500448 0.000000e+00 0 1.000000e+00
## rs577393745 0.000000e+00 0 1.000000e+00
## rs533867263 9.433138e-03 1 9.226276e-01
## rs140793040 0.000000e+00 0 1.000000e+00
## rs3787946 1.991770e+00 1 1.581560e-01
## rs374510753 0.000000e+00 0 1.000000e+00
## rs180784757 9.433138e-03 1 9.226276e-01
## rs371965055 0.000000e+00 0 1.000000e+00
## rs117888036 0.000000e+00 0 1.000000e+00
## rs564705715 0.000000e+00 0 1.000000e+00
## rs368735421 0.000000e+00 0 1.000000e+00
## rs66575656 1.486659e-01 1 6.998135e-01
## rs187662898 0.000000e+00 0 1.000000e+00
## rs529280337 0.000000e+00 0 1.000000e+00
## rs189425119 0.000000e+00 0 1.000000e+00
## rs80027429 0.000000e+00 0 1.000000e+00
## rs2298662 9.433138e-03 1 9.226276e-01
## rs551442787 0.000000e+00 0 1.000000e+00
## rs367879274 0.000000e+00 0 1.000000e+00
## rs8129582 0.000000e+00 0 1.000000e+00
## rs2298663 3.003602e+00 1 8.307961e-02
## rs567520083 0.000000e+00 0 1.000000e+00
## rs62217527 2.142432e-02 1 8.836289e-01
## rs181793632 0.000000e+00 0 1.000000e+00
## rs186097675 0.000000e+00 0 1.000000e+00
## rs545165511 0.000000e+00 0 1.000000e+00
## rs112213575 2.142432e-02 1 8.836289e-01
## rs190516010 0.000000e+00 0 1.000000e+00
## rs527284833 9.433138e-03 1 9.226276e-01
## rs113034290 0.000000e+00 0 1.000000e+00
## rs138995130 0.000000e+00 0 1.000000e+00
## rs111383922 0.000000e+00 0 1.000000e+00
## rs562934488 0.000000e+00 0 1.000000e+00
## rs75756279 0.000000e+00 0 1.000000e+00
## rs201599220 0.000000e+00 0 1.000000e+00
## rs200615061 0.000000e+00 0 1.000000e+00
## rs181414852 0.000000e+00 0 1.000000e+00
## rs547544037 0.000000e+00 0 1.000000e+00
## rs140530035 1.520650e-01 1 6.965697e-01
## rs199999942 0.000000e+00 0 1.000000e+00
## rs74423429 0.000000e+00 0 1.000000e+00
## rs144800717 0.000000e+00 0 1.000000e+00
## rs569681994 0.000000e+00 0 1.000000e+00
## rs538591883 0.000000e+00 0 1.000000e+00
## rs558695627 0.000000e+00 0 1.000000e+00
## rs572136846 0.000000e+00 0 1.000000e+00
## rs9985159 2.698733e+00 1 1.004280e-01
## rs148136016 2.336398e-03 1 9.614482e-01
## rs113506821 0.000000e+00 0 1.000000e+00
## rs543344542 0.000000e+00 0 1.000000e+00
## rs118108663 0.000000e+00 0 1.000000e+00
## rs76973757 2.545009e-01 1 6.139235e-01
## rs2094881 3.719985e+00 1 5.376486e-02
## rs565063101 0.000000e+00 0 1.000000e+00
## rs527811791 0.000000e+00 0 1.000000e+00
## rs547558669 0.000000e+00 0 1.000000e+00
## rs141898436 6.064554e-02 1 8.054785e-01
## rs186275240 0.000000e+00 0 1.000000e+00
## rs191223345 0.000000e+00 0 1.000000e+00
## rs182543256 2.336398e-03 1 9.614482e-01
## rs538655114 0.000000e+00 0 1.000000e+00
## rs552174617 0.000000e+00 0 1.000000e+00
## rs2104810 3.032520e+00 1 8.161117e-02
## rs8131648 1.983646e+00 1 1.590068e-01
## rs554357472 0.000000e+00 0 1.000000e+00
## rs8131649 1.570073e+00 1 2.101964e-01
## rs150633108 0.000000e+00 0 1.000000e+00
## rs140715097 0.000000e+00 0 1.000000e+00
## rs576292358 0.000000e+00 0 1.000000e+00
## rs58978895 3.419130e-03 1 9.533716e-01
## rs8134203 3.032520e+00 1 8.161117e-02
## rs8134216 3.032520e+00 1 8.161117e-02
## rs541563671 0.000000e+00 0 1.000000e+00
## rs561179495 0.000000e+00 0 1.000000e+00
## rs529909772 0.000000e+00 0 1.000000e+00
## rs549676267 0.000000e+00 0 1.000000e+00
## rs563595766 0.000000e+00 0 1.000000e+00
## rs144576889 9.433138e-03 1 9.226276e-01
## rs76855393 0.000000e+00 0 1.000000e+00
## rs565720190 0.000000e+00 0 1.000000e+00
## rs534317621 0.000000e+00 0 1.000000e+00
## rs548267325 0.000000e+00 0 1.000000e+00
## rs61325328 2.336398e-03 1 9.614482e-01
## rs541793866 0.000000e+00 0 1.000000e+00
## rs187723495 0.000000e+00 0 1.000000e+00
## rs576215861 0.000000e+00 0 1.000000e+00
## rs538990932 9.433138e-03 1 9.226276e-01
## rs558824762 0.000000e+00 0 1.000000e+00
## rs115967323 2.336398e-03 1 9.614482e-01
## rs9636988 4.674515e-01 1 4.941619e-01
## rs4818239 1.162369e+00 1 2.809748e-01
## rs574811036 0.000000e+00 0 1.000000e+00
## rs116170128 1.091771e+00 1 2.960791e-01
## rs563507800 0.000000e+00 0 1.000000e+00
## rs532335955 0.000000e+00 0 1.000000e+00
## rs546165791 0.000000e+00 0 1.000000e+00
## rs559637785 0.000000e+00 0 1.000000e+00
## rs139816990 2.142432e-02 1 8.836289e-01
## rs9983330 1.030573e+00 1 3.100241e-01
## rs568167564 0.000000e+00 0 1.000000e+00
## rs116577479 2.545009e-01 1 6.139235e-01
## rs550253934 0.000000e+00 0 1.000000e+00
## rs190906841 0.000000e+00 0 1.000000e+00
## rs149952403 0.000000e+00 0 1.000000e+00
## rs558834486 0.000000e+00 0 1.000000e+00
## rs566208761 0.000000e+00 0 1.000000e+00
## rs535041396 0.000000e+00 0 1.000000e+00
## rs554692214 0.000000e+00 0 1.000000e+00
## rs574759417 0.000000e+00 0 1.000000e+00
## rs145026350 0.000000e+00 0 1.000000e+00
## rs149076631 0.000000e+00 0 1.000000e+00
## rs577113395 0.000000e+00 0 1.000000e+00
## rs543171759 0.000000e+00 0 1.000000e+00
## rs559698510 0.000000e+00 0 1.000000e+00
## rs142125685 0.000000e+00 0 1.000000e+00
## rs183147330 0.000000e+00 0 1.000000e+00
## rs561570726 0.000000e+00 0 1.000000e+00
## rs530527676 0.000000e+00 0 1.000000e+00
## rs73372191 3.443878e-02 1 8.527766e-01
## rs9974995 4.674515e-01 1 4.941619e-01
## rs187432464 0.000000e+00 0 1.000000e+00
## rs192026800 0.000000e+00 0 1.000000e+00
## rs73372193 3.443878e-02 1 8.527766e-01
## rs9974933 4.674515e-01 1 4.941619e-01
## rs184993135 0.000000e+00 0 1.000000e+00
## rs9975014 4.674515e-01 1 4.941619e-01
## rs141788162 0.000000e+00 0 1.000000e+00
## rs201661208 0.000000e+00 0 1.000000e+00
## rs150554820 0.000000e+00 0 1.000000e+00
## rs201109436 0.000000e+00 0 1.000000e+00
## rs146681599 0.000000e+00 0 1.000000e+00
## rs74564819 0.000000e+00 0 1.000000e+00
## rs112985398 3.844785e-02 1 8.445466e-01
## rs145292327 0.000000e+00 0 1.000000e+00
## rs569742981 0.000000e+00 0 1.000000e+00
## rs538022064 0.000000e+00 0 1.000000e+00
## rs117432315 0.000000e+00 0 1.000000e+00
## rs137871202 2.336398e-03 1 9.614482e-01
## rs140636987 0.000000e+00 0 1.000000e+00
## rs199575615 0.000000e+00 0 1.000000e+00
## rs117941520 0.000000e+00 0 1.000000e+00
## rs528681261 0.000000e+00 0 1.000000e+00
## rs915823 1.850278e+00 1 1.737510e-01
## rs138365638 2.142432e-02 1 8.836289e-01
## rs553596811 0.000000e+00 0 1.000000e+00
## rs192793047 0.000000e+00 0 1.000000e+00
## rs184380117 0.000000e+00 0 1.000000e+00
## rs570835062 0.000000e+00 0 1.000000e+00
## rs539981354 0.000000e+00 0 1.000000e+00
## rs553431242 0.000000e+00 0 1.000000e+00
## rs28403625 0.000000e+00 0 1.000000e+00
## rs535599066 0.000000e+00 0 1.000000e+00
## rs181979537 0.000000e+00 0 1.000000e+00
## rs575167859 0.000000e+00 0 1.000000e+00
## rs875393 0.000000e+00 0 1.000000e+00
## rs184764113 0.000000e+00 0 1.000000e+00
## rs150473664 0.000000e+00 0 1.000000e+00
## rs138282462 0.000000e+00 0 1.000000e+00
## rs149601802 0.000000e+00 0 1.000000e+00
## rs62217529 0.000000e+00 0 1.000000e+00
## rs548725701 0.000000e+00 0 1.000000e+00
## rs139181746 0.000000e+00 0 1.000000e+00
## rs376897988 0.000000e+00 0 1.000000e+00
## rs551277697 0.000000e+00 0 1.000000e+00
## rs734055 0.000000e+00 0 1.000000e+00
## rs143109187 0.000000e+00 0 1.000000e+00
## rs546790067 0.000000e+00 0 1.000000e+00
## rs566925007 0.000000e+00 0 1.000000e+00
## rs34561135 0.000000e+00 0 1.000000e+00
## rs734056 1.637992e+00 1 2.006012e-01
## rs186734573 0.000000e+00 0 1.000000e+00
## rs79397218 0.000000e+00 0 1.000000e+00
## rs61735789 0.000000e+00 0 1.000000e+00
## rs12329760 1.030573e+00 1 3.100241e-01
## rs371583288 0.000000e+00 0 1.000000e+00
## rs147465180 0.000000e+00 0 1.000000e+00
## rs576916247 2.142432e-02 1 8.836289e-01
## rs378501 0.000000e+00 0 1.000000e+00
## rs542333608 0.000000e+00 0 1.000000e+00
## rs142710225 2.336398e-03 1 9.614482e-01
## rs150581606 0.000000e+00 0 1.000000e+00
## rs370654859 9.433138e-03 1 9.226276e-01
## rs191228684 0.000000e+00 0 1.000000e+00
## rs181107522 2.336398e-03 1 9.614482e-01
## rs185166990 0.000000e+00 0 1.000000e+00
## rs566732470 0.000000e+00 0 1.000000e+00
## rs529424146 0.000000e+00 0 1.000000e+00
## rs112209215 2.142432e-02 1 8.836289e-01
## rs75833467 0.000000e+00 0 1.000000e+00
## rs370641046 0.000000e+00 0 1.000000e+00
## rs75655573 2.142432e-02 1 8.836289e-01
## rs4290734 2.545009e-01 1 6.139235e-01
## rs533571483 0.000000e+00 0 1.000000e+00
## rs553755767 0.000000e+00 0 1.000000e+00
## rs189582211 0.000000e+00 0 1.000000e+00
## rs542746297 0.000000e+00 0 1.000000e+00
## rs139374762 1.388883e+00 1 2.385938e-01
## rs556249043 0.000000e+00 0 1.000000e+00
## rs149708827 0.000000e+00 0 1.000000e+00
## rs544784789 0.000000e+00 0 1.000000e+00
## rs564795317 0.000000e+00 0 1.000000e+00
## rs533708336 9.433138e-03 1 9.226276e-01
## rs183405398 0.000000e+00 0 1.000000e+00
## rs2410428 7.974953e-02 1 7.776371e-01
## rs116568213 3.844785e-02 1 8.445466e-01
## rs549067055 0.000000e+00 0 1.000000e+00
## rs146723217 3.109578e-01 1 5.770931e-01
## rs531418415 0.000000e+00 0 1.000000e+00
## rs190682490 0.000000e+00 0 1.000000e+00
## rs555600211 0.000000e+00 0 1.000000e+00
## rs183385443 0.000000e+00 0 1.000000e+00
## rs188495261 0.000000e+00 0 1.000000e+00
## rs567258163 0.000000e+00 0 1.000000e+00
## rs530689318 0.000000e+00 0 1.000000e+00
## rs535834505 0.000000e+00 0 1.000000e+00
## rs200072801 0.000000e+00 0 1.000000e+00
## rs576181035 0.000000e+00 0 1.000000e+00
## rs111620846 2.142432e-02 1 8.836289e-01
## rs2276205 1.091771e+00 1 2.960791e-01
## rs578103765 0.000000e+00 0 1.000000e+00
## rs373311004 0.000000e+00 0 1.000000e+00
## rs560575827 0.000000e+00 0 1.000000e+00
## rs573986970 0.000000e+00 0 1.000000e+00
## rs542577849 0.000000e+00 0 1.000000e+00
## rs562664177 0.000000e+00 0 1.000000e+00
## rs144033559 3.844785e-02 1 8.445466e-01
## rs551442751 0.000000e+00 0 1.000000e+00
## rs118134524 0.000000e+00 0 1.000000e+00
## rs527489879 0.000000e+00 0 1.000000e+00
## rs541239376 0.000000e+00 0 1.000000e+00
## rs112753686 2.142432e-02 1 8.836289e-01
## rs536278410 0.000000e+00 0 1.000000e+00
## rs549833072 0.000000e+00 0 1.000000e+00
## rs192709500 0.000000e+00 0 1.000000e+00
## rs34783969 1.687500e+00 1 1.939309e-01
## rs558683527 0.000000e+00 0 1.000000e+00
## rs149024920 2.336398e-03 1 9.614482e-01
## rs184767756 0.000000e+00 0 1.000000e+00
## rs554462605 0.000000e+00 0 1.000000e+00
## rs574223231 0.000000e+00 0 1.000000e+00
## rs542888854 2.336398e-03 1 9.614482e-01
## rs562574628 0.000000e+00 0 1.000000e+00
## rs576016640 0.000000e+00 0 1.000000e+00
## rs191353873 0.000000e+00 0 1.000000e+00
## rs143818732 0.000000e+00 0 1.000000e+00
## rs527499027 0.000000e+00 0 1.000000e+00
## rs540843636 0.000000e+00 0 1.000000e+00
## rs9976780 2.636719e+00 1 1.044190e-01
## rs530051850 0.000000e+00 0 1.000000e+00
## rs550091217 0.000000e+00 0 1.000000e+00
## rs569716669 0.000000e+00 0 1.000000e+00
## rs185114529 0.000000e+00 0 1.000000e+00
## rs552256264 0.000000e+00 0 1.000000e+00
## rs565603164 0.000000e+00 0 1.000000e+00
## rs534423780 0.000000e+00 0 1.000000e+00
## rs554097496 0.000000e+00 0 1.000000e+00
## rs567504414 0.000000e+00 0 1.000000e+00
## rs9977234 1.564176e-01 1 6.924764e-01
## rs556435793 0.000000e+00 0 1.000000e+00
## rs576030223 2.336398e-03 1 9.614482e-01
## rs9979311 1.936974e+00 1 1.639975e-01
## rs564092110 2.581501e+00 1 1.081191e-01
## rs9980225 0.000000e+00 0 1.000000e+00
## rs75430506 0.000000e+00 0 1.000000e+00
## rs541499151 0.000000e+00 0 1.000000e+00
## rs189436980 0.000000e+00 0 1.000000e+00
## rs2410429 5.675584e-01 1 4.512308e-01
## rs182231135 0.000000e+00 0 1.000000e+00
## rs148155433 2.336398e-03 1 9.614482e-01
## rs532159311 2.336398e-03 1 9.614482e-01
## rs141091900 0.000000e+00 0 1.000000e+00
## rs2410430 2.336398e-03 1 9.614482e-01
## rs528206157 0.000000e+00 0 1.000000e+00
## rs185188053 0.000000e+00 0 1.000000e+00
## rs146865760 0.000000e+00 0 1.000000e+00
## rs536399456 0.000000e+00 0 1.000000e+00
## rs556456416 0.000000e+00 0 1.000000e+00
## rs569988060 0.000000e+00 0 1.000000e+00
## rs139222305 0.000000e+00 0 1.000000e+00
## rs558715322 0.000000e+00 0 1.000000e+00
## rs143148203 2.336398e-03 1 9.614482e-01
## rs541215881 0.000000e+00 0 1.000000e+00
## rs555056677 0.000000e+00 0 1.000000e+00
## rs56066678 2.251166e-01 1 6.351686e-01
## rs543607180 0.000000e+00 0 1.000000e+00
## rs563595464 0.000000e+00 0 1.000000e+00
## rs190134102 0.000000e+00 0 1.000000e+00
## rs2838039 3.432309e+00 1 6.393265e-02
## rs559125514 0.000000e+00 0 1.000000e+00
## rs181332845 2.336398e-03 1 9.614482e-01
## rs147054075 0.000000e+00 0 1.000000e+00
## rs561404433 0.000000e+00 0 1.000000e+00
## rs530253044 0.000000e+00 0 1.000000e+00
## rs549990870 0.000000e+00 0 1.000000e+00
## rs569899862 0.000000e+00 0 1.000000e+00
## rs138164084 2.336398e-03 1 9.614482e-01
## rs2156300 0.000000e+00 0 1.000000e+00
## rs565946159 0.000000e+00 0 1.000000e+00
## rs534978728 0.000000e+00 0 1.000000e+00
## rs79243099 8.816327e-02 1 7.665253e-01
## rs574957059 0.000000e+00 0 1.000000e+00
## rs185406661 0.000000e+00 0 1.000000e+00
## rs557224603 0.000000e+00 0 1.000000e+00
## rs576955441 0.000000e+00 0 1.000000e+00
## rs545825254 0.000000e+00 0 1.000000e+00
## rs2156301 9.433138e-03 1 9.226276e-01
## rs188263047 9.433138e-03 1 9.226276e-01
## rs541782164 0.000000e+00 0 1.000000e+00
## rs561327241 0.000000e+00 0 1.000000e+00
## rs530553979 0.000000e+00 0 1.000000e+00
## rs550252022 0.000000e+00 0 1.000000e+00
## rs144154504 0.000000e+00 0 1.000000e+00
## rs532694622 0.000000e+00 0 1.000000e+00
## rs552607028 0.000000e+00 0 1.000000e+00
## rs180826598 0.000000e+00 0 1.000000e+00
## rs186422083 0.000000e+00 0 1.000000e+00
## rs559934921 0.000000e+00 0 1.000000e+00
## rs148701953 0.000000e+00 0 1.000000e+00
## rs141324737 0.000000e+00 0 1.000000e+00
## rs146252393 0.000000e+00 0 1.000000e+00
## rs4818240 3.109578e-01 1 5.770931e-01
## rs9983252 4.121729e+00 1 4.233578e-02
## rs4818241 0.000000e+00 0 1.000000e+00
## rs376454804 0.000000e+00 0 1.000000e+00
## rs573121598 0.000000e+00 0 1.000000e+00
## rs4816720 3.109578e-01 1 5.770931e-01
## rs554986094 0.000000e+00 0 1.000000e+00
## rs574890845 0.000000e+00 0 1.000000e+00
## rs543908797 0.000000e+00 0 1.000000e+00
## rs4818242 0.000000e+00 0 1.000000e+00
## rs567039032 0.000000e+00 0 1.000000e+00
## rs11702475 1.764378e-01 1 6.744526e-01
## rs35871560 2.232638e-02 1 8.812221e-01
## rs559668949 0.000000e+00 0 1.000000e+00
## rs186847776 0.000000e+00 0 1.000000e+00
## rs381179 0.000000e+00 0 1.000000e+00
## rs2257202 6.808736e-01 1 4.092861e-01
## rs118108194 0.000000e+00 0 1.000000e+00
## rs551103429 0.000000e+00 0 1.000000e+00
## rs570889046 2.336398e-03 1 9.614482e-01
## rs539528622 0.000000e+00 0 1.000000e+00
## rs553167698 0.000000e+00 0 1.000000e+00
## rs62217531 9.854388e-02 1 7.535841e-01
## rs535407146 0.000000e+00 0 1.000000e+00
## rs555198424 0.000000e+00 0 1.000000e+00
## rs144052153 0.000000e+00 0 1.000000e+00
## rs2298857 4.573750e-02 1 8.306537e-01
## rs557592331 0.000000e+00 0 1.000000e+00
## rs554175681 0.000000e+00 0 1.000000e+00
## rs138928680 2.336398e-03 1 9.614482e-01
## rs182791541 0.000000e+00 0 1.000000e+00
## rs539934676 0.000000e+00 0 1.000000e+00
## rs534648425 9.433138e-03 1 9.226276e-01
## rs559934709 9.433138e-03 1 9.226276e-01
## rs528691645 0.000000e+00 0 1.000000e+00
## rs187175193 0.000000e+00 0 1.000000e+00
## rs562327565 0.000000e+00 0 1.000000e+00
## rs530942388 0.000000e+00 0 1.000000e+00
## rs383510 3.478387e-01 1 5.553391e-01
## rs564834372 0.000000e+00 0 1.000000e+00
## rs139467735 7.881619e-01 1 3.746565e-01
## rs546775580 0.000000e+00 0 1.000000e+00
## rs566662405 0.000000e+00 0 1.000000e+00
## rs535235411 0.000000e+00 0 1.000000e+00
## rs191595416 0.000000e+00 0 1.000000e+00
## rs149676870 0.000000e+00 0 1.000000e+00
## rs537728755 2.336398e-03 1 9.614482e-01
## rs557766326 0.000000e+00 0 1.000000e+00
## rs577554253 0.000000e+00 0 1.000000e+00
## rs533931558 0.000000e+00 0 1.000000e+00
## rs193122175 0.000000e+00 0 1.000000e+00
## rs573484758 0.000000e+00 0 1.000000e+00
## rs542407537 0.000000e+00 0 1.000000e+00
## rs562387776 0.000000e+00 0 1.000000e+00
## rs576042280 0.000000e+00 0 1.000000e+00
## rs364289 3.889116e-01 1 5.328721e-01
## rs370363307 0.000000e+00 0 1.000000e+00
## rs149226336 0.000000e+00 0 1.000000e+00
## rs547082861 0.000000e+00 0 1.000000e+00
## rs560182574 0.000000e+00 0 1.000000e+00
## rs79468500 8.816327e-02 1 7.665253e-01
## rs528477904 0.000000e+00 0 1.000000e+00
## rs185555838 2.336398e-03 1 9.614482e-01
## rs537787922 0.000000e+00 0 1.000000e+00
## rs80275470 0.000000e+00 0 1.000000e+00
## rs365025 7.934694e-01 1 3.730531e-01
## rs533513206 0.000000e+00 0 1.000000e+00
## rs553847722 0.000000e+00 0 1.000000e+00
## rs190125027 0.000000e+00 0 1.000000e+00
## rs79512425 0.000000e+00 0 1.000000e+00
## rs556037160 0.000000e+00 0 1.000000e+00
## rs79517809 9.433138e-03 1 9.226276e-01
## rs193067129 0.000000e+00 0 1.000000e+00
## rs184149889 0.000000e+00 0 1.000000e+00
## rs365724 8.975764e-01 1 3.434324e-01
## rs56097233 8.975764e-01 1 3.434324e-01
## rs554139014 0.000000e+00 0 1.000000e+00
## rs563305848 0.000000e+00 0 1.000000e+00
## rs188508398 0.000000e+00 0 1.000000e+00
## rs147934552 0.000000e+00 0 1.000000e+00
## rs141722242 2.336398e-03 1 9.614482e-01
## rs181051004 0.000000e+00 0 1.000000e+00
## rs372422423 2.336398e-03 1 9.614482e-01
## rs185086499 0.000000e+00 0 1.000000e+00
## rs113670863 0.000000e+00 0 1.000000e+00
## rs571417292 0.000000e+00 0 1.000000e+00
## rs533805020 0.000000e+00 0 1.000000e+00
## rs375408 0.000000e+00 0 1.000000e+00
## rs402197 9.433138e-03 1 9.226276e-01
## rs375760 2.692645e-30 1 1.000000e+00
## rs182175876 0.000000e+00 0 1.000000e+00
## rs147063811 0.000000e+00 0 1.000000e+00
## rs538022362 0.000000e+00 0 1.000000e+00
## rs558449333 0.000000e+00 0 1.000000e+00
## rs570264808 2.336398e-03 1 9.614482e-01
## rs540908584 0.000000e+00 0 1.000000e+00
## rs553788077 0.000000e+00 0 1.000000e+00
## rs415918 1.501186e+00 1 2.204889e-01
## rs542712985 0.000000e+00 0 1.000000e+00
## rs562599408 0.000000e+00 0 1.000000e+00
## rs531548362 0.000000e+00 0 1.000000e+00
## rs544847277 0.000000e+00 0 1.000000e+00
## rs141888586 1.091771e+00 1 2.960791e-01
## rs527623581 0.000000e+00 0 1.000000e+00
## rs150205743 0.000000e+00 0 1.000000e+00
## rs567275769 0.000000e+00 0 1.000000e+00
## rs561575371 0.000000e+00 0 1.000000e+00
## rs528451616 2.336398e-03 1 9.614482e-01
## rs56136037 0.000000e+00 0 1.000000e+00
## rs138705806 0.000000e+00 0 1.000000e+00
## rs558104085 0.000000e+00 0 1.000000e+00
## rs571805265 0.000000e+00 0 1.000000e+00
## rs115265507 6.064554e-02 1 8.054785e-01
## rs573704321 0.000000e+00 0 1.000000e+00
## rs542575868 0.000000e+00 0 1.000000e+00
## rs556260609 0.000000e+00 0 1.000000e+00
## rs371531071 0.000000e+00 0 1.000000e+00
## rs422471 1.119920e+00 1 2.899357e-01
## rs372286621 0.000000e+00 0 1.000000e+00
## rs527435310 0.000000e+00 0 1.000000e+00
## rs386416 1.119920e+00 1 2.899357e-01
## rs187290362 0.000000e+00 0 1.000000e+00
## rs3819138 0.000000e+00 0 1.000000e+00
## rs549832756 0.000000e+00 0 1.000000e+00
## rs536039173 0.000000e+00 0 1.000000e+00
## rs555995855 0.000000e+00 0 1.000000e+00
## rs117696554 0.000000e+00 0 1.000000e+00
## rs191693032 0.000000e+00 0 1.000000e+00
## rs528452128 2.336398e-03 1 9.614482e-01
## rs141864728 0.000000e+00 0 1.000000e+00
## rs3787947 1.119920e+00 1 2.899357e-01
## rs139305247 0.000000e+00 0 1.000000e+00
## rs536525396 9.433138e-03 1 9.226276e-01
## rs142518005 2.336398e-03 1 9.614482e-01
## rs183016576 0.000000e+00 0 1.000000e+00
## rs423596 0.000000e+00 0 1.000000e+00
## rs145996484 0.000000e+00 0 1.000000e+00
## rs572164476 0.000000e+00 0 1.000000e+00
## rs541200538 0.000000e+00 0 1.000000e+00
## rs561043908 0.000000e+00 0 1.000000e+00
## rs185946013 0.000000e+00 0 1.000000e+00
## rs543734531 0.000000e+00 0 1.000000e+00
## rs563186329 0.000000e+00 0 1.000000e+00
## rs532082574 0.000000e+00 0 1.000000e+00
## rs2298664 1.119920e+00 1 2.899357e-01
## rs559363556 0.000000e+00 0 1.000000e+00
## rs527898967 0.000000e+00 0 1.000000e+00
## rs142103907 9.433138e-03 1 9.226276e-01
## rs2298665 0.000000e+00 0 1.000000e+00
## rs536589154 0.000000e+00 0 1.000000e+00
## rs429442 5.232693e-01 1 4.694504e-01
## rs569959837 0.000000e+00 0 1.000000e+00
## rs562131885 9.433138e-03 1 9.226276e-01
## rs151152524 0.000000e+00 0 1.000000e+00
## rs190265904 0.000000e+00 0 1.000000e+00
## rs147711290 0.000000e+00 0 1.000000e+00
## rs139001727 0.000000e+00 0 1.000000e+00
## rs144948620 0.000000e+00 0 1.000000e+00
## rs574683527 0.000000e+00 0 1.000000e+00
## rs543740889 0.000000e+00 0 1.000000e+00
## rs563555305 0.000000e+00 0 1.000000e+00
## rs576777264 0.000000e+00 0 1.000000e+00
## rs147945120 0.000000e+00 0 1.000000e+00
## rs373622147 0.000000e+00 0 1.000000e+00
## rs528205318 0.000000e+00 0 1.000000e+00
## rs141764184 0.000000e+00 0 1.000000e+00
## rs561192574 0.000000e+00 0 1.000000e+00
## rs530193039 0.000000e+00 0 1.000000e+00
## rs550222907 0.000000e+00 0 1.000000e+00
## rs570302987 0.000000e+00 0 1.000000e+00
## rs192666663 0.000000e+00 0 1.000000e+00
## rs147099383 8.816327e-02 1 7.665253e-01
## rs371020678 0.000000e+00 0 1.000000e+00
## rs535174935 0.000000e+00 0 1.000000e+00
## rs551060886 2.989617e+00 1 8.380000e-02
## rs555056776 0.000000e+00 0 1.000000e+00
## rs568267971 0.000000e+00 0 1.000000e+00
## rs537212420 0.000000e+00 0 1.000000e+00
## rs556994660 0.000000e+00 0 1.000000e+00
## rs145347702 0.000000e+00 0 1.000000e+00
## rs11911394 8.416865e-01 1 3.589149e-01
## rs552589240 0.000000e+00 0 1.000000e+00
## rs115849825 2.336398e-03 1 9.614482e-01
## rs189040987 0.000000e+00 0 1.000000e+00
## rs561412657 0.000000e+00 0 1.000000e+00
## rs192259532 0.000000e+00 0 1.000000e+00
## rs544023509 0.000000e+00 0 1.000000e+00
## rs928871 2.391696e-01 1 6.248066e-01
## rs532882572 0.000000e+00 0 1.000000e+00
## rs116865960 2.336398e-03 1 9.614482e-01
## rs566085188 0.000000e+00 0 1.000000e+00
## rs528452414 0.000000e+00 0 1.000000e+00
## rs548713344 0.000000e+00 0 1.000000e+00
## rs115596471 1.211511e-01 1 7.277894e-01
## rs149173609 0.000000e+00 0 1.000000e+00
## rs184365262 0.000000e+00 0 1.000000e+00
## rs570544092 0.000000e+00 0 1.000000e+00
## rs539767738 0.000000e+00 0 1.000000e+00
## rs6517669 1.541892e+00 1 2.143364e-01
## rs187208295 0.000000e+00 0 1.000000e+00
## rs541396507 0.000000e+00 0 1.000000e+00
## rs117656646 0.000000e+00 0 1.000000e+00
## rs577944614 0.000000e+00 0 1.000000e+00
## rs143270468 0.000000e+00 0 1.000000e+00
## rs148304071 0.000000e+00 0 1.000000e+00
## rs417443 0.000000e+00 0 1.000000e+00
## rs546531769 0.000000e+00 0 1.000000e+00
## rs140458174 0.000000e+00 0 1.000000e+00
## rs145570856 1.211511e-01 1 7.277894e-01
## rs35899679 5.156998e-02 1 8.203536e-01
## rs562045100 0.000000e+00 0 1.000000e+00
## rs417888 1.564176e-01 1 6.924764e-01
## rs35041537 1.564176e-01 1 6.924764e-01
## rs570629130 2.336398e-03 1 9.614482e-01
## rs539218228 0.000000e+00 0 1.000000e+00
## rs546512831 0.000000e+00 0 1.000000e+00
## rs566457884 0.000000e+00 0 1.000000e+00
## rs535348902 0.000000e+00 0 1.000000e+00
## rs554986236 0.000000e+00 0 1.000000e+00
## rs115429336 2.545009e-01 1 6.139235e-01
## rs191763249 0.000000e+00 0 1.000000e+00
## rs557531433 0.000000e+00 0 1.000000e+00
## rs577381045 0.000000e+00 0 1.000000e+00
## rs546388820 0.000000e+00 0 1.000000e+00
## rs138161288 0.000000e+00 0 1.000000e+00
## rs573613070 2.336398e-03 1 9.614482e-01
## rs542530281 0.000000e+00 0 1.000000e+00
## rs562157895 0.000000e+00 0 1.000000e+00
## rs10154090 4.320000e+00 1 3.766692e-02
## rs183851466 0.000000e+00 0 1.000000e+00
## rs73357642 1.608864e-01 1 6.883417e-01
## rs533445403 0.000000e+00 0 1.000000e+00
## rs391099 3.701149e-01 1 5.429416e-01
## rs566323147 0.000000e+00 0 1.000000e+00
## rs73357644 1.608864e-01 1 6.883417e-01
## rs142659685 0.000000e+00 0 1.000000e+00
## rs568846685 0.000000e+00 0 1.000000e+00
## rs537412369 0.000000e+00 0 1.000000e+00
## rs557853844 0.000000e+00 0 1.000000e+00
## rs577689706 0.000000e+00 0 1.000000e+00
## rs9305745 5.061791e+00 1 2.445894e-02
## rs189298425 0.000000e+00 0 1.000000e+00
## rs180792570 0.000000e+00 0 1.000000e+00
## rs34769294 1.887624e-01 1 6.639487e-01
## rs392370 8.533594e-02 1 7.701925e-01
## rs562203987 0.000000e+00 0 1.000000e+00
## rs79617378 0.000000e+00 0 1.000000e+00
## rs544308006 0.000000e+00 0 1.000000e+00
## rs555667242 0.000000e+00 0 1.000000e+00
## rs117688613 0.000000e+00 0 1.000000e+00
## rs374261644 9.433138e-03 1 9.226276e-01
## rs185726773 2.336398e-03 1 9.614482e-01
## rs149275684 0.000000e+00 0 1.000000e+00
## rs144481222 2.336398e-03 1 9.614482e-01
## rs569026066 0.000000e+00 0 1.000000e+00
## rs531609155 2.336398e-03 1 9.614482e-01
## rs551122576 0.000000e+00 0 1.000000e+00
## rs571211518 0.000000e+00 0 1.000000e+00
## rs2838040 8.509473e-01 1 3.562845e-01
## rs139458208 2.336398e-03 1 9.614482e-01
## rs430915 5.156998e-02 1 8.203536e-01
## rs536023174 0.000000e+00 0 1.000000e+00
## rs401371 3.737024e-01 1 5.409930e-01
## rs144988776 0.000000e+00 0 1.000000e+00
## rs544880781 0.000000e+00 0 1.000000e+00
## rs190685013 0.000000e+00 0 1.000000e+00
## rs577570888 0.000000e+00 0 1.000000e+00
## rs181778800 0.000000e+00 0 1.000000e+00
## rs147138431 0.000000e+00 0 1.000000e+00
## rs386638 0.000000e+00 0 1.000000e+00
## rs76314085 0.000000e+00 0 1.000000e+00
## rs402303 1.727439e+00 1 1.887382e-01
## rs531545827 0.000000e+00 0 1.000000e+00
## rs375827195 0.000000e+00 0 1.000000e+00
## rs112517486 2.336398e-03 1 9.614482e-01
## rs527242422 0.000000e+00 0 1.000000e+00
## rs7275220 2.587902e-01 1 6.109526e-01
## rs567340892 2.336398e-03 1 9.614482e-01
## rs555725136 0.000000e+00 0 1.000000e+00
## rs569432639 2.336398e-03 1 9.614482e-01
## rs538161583 0.000000e+00 0 1.000000e+00
## rs141027872 0.000000e+00 0 1.000000e+00
## rs577632959 0.000000e+00 0 1.000000e+00
## rs540209049 0.000000e+00 0 1.000000e+00
## rs149798653 0.000000e+00 0 1.000000e+00
## rs573828363 0.000000e+00 0 1.000000e+00
## rs145738510 0.000000e+00 0 1.000000e+00
## rs562587205 0.000000e+00 0 1.000000e+00
## rs140605461 0.000000e+00 0 1.000000e+00
## rs544983081 0.000000e+00 0 1.000000e+00
## rs562000328 0.000000e+00 0 1.000000e+00
## rs527305819 0.000000e+00 0 1.000000e+00
## rs145662876 0.000000e+00 0 1.000000e+00
## rs566903241 0.000000e+00 0 1.000000e+00
## rs371497524 0.000000e+00 0 1.000000e+00
## rs549801931 0.000000e+00 0 1.000000e+00
## rs567208488 0.000000e+00 0 1.000000e+00
## rs569362310 0.000000e+00 0 1.000000e+00
## rs113773731 2.336398e-03 1 9.614482e-01
## rs551681725 0.000000e+00 0 1.000000e+00
## rs2838041 6.808736e-01 1 4.092861e-01
## rs28369457 0.000000e+00 0 1.000000e+00
## rs553787982 0.000000e+00 0 1.000000e+00
## rs573692135 0.000000e+00 0 1.000000e+00
## rs186425543 0.000000e+00 0 1.000000e+00
## rs143167369 0.000000e+00 0 1.000000e+00
## rs190970011 0.000000e+00 0 1.000000e+00
## rs376235035 0.000000e+00 0 1.000000e+00
## rs2838042 1.312079e-02 1 9.088050e-01
## rs572410855 0.000000e+00 0 1.000000e+00
## rs373196115 0.000000e+00 0 1.000000e+00
## rs3787950 7.554645e-01 1 3.847516e-01
## rs61735793 0.000000e+00 0 1.000000e+00
## rs114363287 2.336398e-03 1 9.614482e-01
## rs200169208 0.000000e+00 0 1.000000e+00
## rs532085624 0.000000e+00 0 1.000000e+00
## rs61735792 0.000000e+00 0 1.000000e+00
## rs149527323 2.336398e-03 1 9.614482e-01
## rs201679623 0.000000e+00 0 1.000000e+00
## rs547639377 0.000000e+00 0 1.000000e+00
## rs376143876 0.000000e+00 0 1.000000e+00
## rs138651919 0.000000e+00 0 1.000000e+00
## rs187460831 0.000000e+00 0 1.000000e+00
## rs141232947 0.000000e+00 0 1.000000e+00
## rs199824558 0.000000e+00 0 1.000000e+00
## rs201093031 0.000000e+00 0 1.000000e+00
## rs61735791 0.000000e+00 0 1.000000e+00
## rs141685390 0.000000e+00 0 1.000000e+00
## rs554868303 0.000000e+00 0 1.000000e+00
## rs61735790 2.336398e-03 1 9.614482e-01
## rs199636550 0.000000e+00 0 1.000000e+00
## rs138056981 0.000000e+00 0 1.000000e+00
## rs564806692 0.000000e+00 0 1.000000e+00
## rs532276666 0.000000e+00 0 1.000000e+00
## rs143523726 0.000000e+00 0 1.000000e+00
## rs565142477 0.000000e+00 0 1.000000e+00
## rs147147950 0.000000e+00 0 1.000000e+00
## rs547701911 2.336398e-03 1 9.614482e-01
## rs374203194 0.000000e+00 0 1.000000e+00
## rs530306757 0.000000e+00 0 1.000000e+00
## rs140259165 0.000000e+00 0 1.000000e+00
## rs569947342 0.000000e+00 0 1.000000e+00
## rs150066796 0.000000e+00 0 1.000000e+00
## rs144364265 0.000000e+00 0 1.000000e+00
## rs565896303 0.000000e+00 0 1.000000e+00
## rs190899605 0.000000e+00 0 1.000000e+00
## rs554443322 0.000000e+00 0 1.000000e+00
## rs574708698 0.000000e+00 0 1.000000e+00
## rs543339705 0.000000e+00 0 1.000000e+00
## rs556966925 2.336398e-03 1 9.614482e-01
## rs9975782 0.000000e+00 0 1.000000e+00
## rs545504807 0.000000e+00 0 1.000000e+00
## rs558991814 0.000000e+00 0 1.000000e+00
## rs183141812 0.000000e+00 0 1.000000e+00
## rs561709263 0.000000e+00 0 1.000000e+00
## rs561308071 0.000000e+00 0 1.000000e+00
## rs530471976 0.000000e+00 0 1.000000e+00
## rs114911304 6.064554e-02 1 8.054785e-01
## rs570041169 0.000000e+00 0 1.000000e+00
## rs532513568 0.000000e+00 0 1.000000e+00
## rs552725417 0.000000e+00 0 1.000000e+00
## rs566203795 0.000000e+00 0 1.000000e+00
## rs535103772 0.000000e+00 0 1.000000e+00
## rs548402221 0.000000e+00 0 1.000000e+00
## rs568162250 2.336398e-03 1 9.614482e-01
## rs537357445 0.000000e+00 0 1.000000e+00
## rs556937466 0.000000e+00 0 1.000000e+00
## rs576447507 0.000000e+00 0 1.000000e+00
## rs539019775 0.000000e+00 0 1.000000e+00
## rs58146697 9.367695e-02 1 7.595537e-01
## rs187426170 0.000000e+00 0 1.000000e+00
## rs541550816 0.000000e+00 0 1.000000e+00
## rs192019778 0.000000e+00 0 1.000000e+00
## rs574894406 0.000000e+00 0 1.000000e+00
## rs544110139 0.000000e+00 0 1.000000e+00
## rs563800355 0.000000e+00 0 1.000000e+00
## rs142425263 0.000000e+00 0 1.000000e+00
## rs552635816 0.000000e+00 0 1.000000e+00
## rs57901138 1.307071e+00 1 2.529256e-01
## rs559771004 0.000000e+00 0 1.000000e+00
## rs528640390 2.336398e-03 1 9.614482e-01
## rs548349124 0.000000e+00 0 1.000000e+00
## rs568256706 0.000000e+00 0 1.000000e+00
## rs184876485 0.000000e+00 0 1.000000e+00
## rs550874821 0.000000e+00 0 1.000000e+00
## rs146845793 0.000000e+00 0 1.000000e+00
## rs539333191 0.000000e+00 0 1.000000e+00
## rs552591104 0.000000e+00 0 1.000000e+00
## rs55760462 8.510729e-02 1 7.704919e-01
## rs535545388 0.000000e+00 0 1.000000e+00
## rs140625413 0.000000e+00 0 1.000000e+00
## rs188241223 0.000000e+00 0 1.000000e+00
## rs543681861 0.000000e+00 0 1.000000e+00
## rs563853944 0.000000e+00 0 1.000000e+00
## rs577650141 0.000000e+00 0 1.000000e+00
## rs564210410 2.336398e-03 1 9.614482e-01
## rs559790300 0.000000e+00 0 1.000000e+00
## rs144276163 0.000000e+00 0 1.000000e+00
## rs548680244 0.000000e+00 0 1.000000e+00
## rs562014334 0.000000e+00 0 1.000000e+00
## rs531013422 0.000000e+00 0 1.000000e+00
## rs550537161 0.000000e+00 0 1.000000e+00
## rs914184 0.000000e+00 0 1.000000e+00
## rs78217567 0.000000e+00 0 1.000000e+00
## rs546447704 0.000000e+00 0 1.000000e+00
## rs566146941 0.000000e+00 0 1.000000e+00
## rs184500277 0.000000e+00 0 1.000000e+00
## rs189067157 0.000000e+00 0 1.000000e+00
## rs181972300 2.336398e-03 1 9.614482e-01
## rs537935645 0.000000e+00 0 1.000000e+00
## rs557632131 0.000000e+00 0 1.000000e+00
## rs34256269 0.000000e+00 0 1.000000e+00
## rs455045 2.597386e+00 1 1.070401e-01
## rs553672299 0.000000e+00 0 1.000000e+00
## rs573213706 0.000000e+00 0 1.000000e+00
## rs541887371 0.000000e+00 0 1.000000e+00
## rs34983238 2.545009e-01 1 6.139235e-01
## rs369342013 0.000000e+00 0 1.000000e+00
## rs544297320 0.000000e+00 0 1.000000e+00
## rs563948676 0.000000e+00 0 1.000000e+00
## rs415731 1.893508e+00 1 1.688068e-01
## rs546510868 0.000000e+00 0 1.000000e+00
## rs1003030 7.881619e-01 1 3.746565e-01
## rs184859933 0.000000e+00 0 1.000000e+00
## rs548742822 0.000000e+00 0 1.000000e+00
## rs562332833 0.000000e+00 0 1.000000e+00
## rs189181802 9.433138e-03 1 9.226276e-01
## rs56695953 2.041588e-01 1 6.513847e-01
## rs571074719 2.336398e-03 1 9.614482e-01
## rs148901354 0.000000e+00 0 1.000000e+00
## rs575335056 0.000000e+00 0 1.000000e+00
## rs367885466 0.000000e+00 0 1.000000e+00
## rs557606231 2.336398e-03 1 9.614482e-01
## rs573493693 0.000000e+00 0 1.000000e+00
## rs542500773 0.000000e+00 0 1.000000e+00
## rs555525436 0.000000e+00 0 1.000000e+00
## rs181592444 0.000000e+00 0 1.000000e+00
## rs544315388 0.000000e+00 0 1.000000e+00
## rs115968373 2.336398e-03 1 9.614482e-01
## rs533229030 0.000000e+00 0 1.000000e+00
## rs540070564 0.000000e+00 0 1.000000e+00
## rs143576576 0.000000e+00 0 1.000000e+00
## rs529218967 0.000000e+00 0 1.000000e+00
## rs138094318 0.000000e+00 0 1.000000e+00
## rs186510586 0.000000e+00 0 1.000000e+00
## rs531300154 0.000000e+00 0 1.000000e+00
## rs11701576 7.881619e-01 1 3.746565e-01
## rs571396442 0.000000e+00 0 1.000000e+00
## rs143562195 2.336398e-03 1 9.614482e-01
## rs547257017 0.000000e+00 0 1.000000e+00
## rs145841350 0.000000e+00 0 1.000000e+00
## rs535778149 0.000000e+00 0 1.000000e+00
## rs138661151 2.142432e-02 1 8.836289e-01
## rs141128014 2.336398e-03 1 9.614482e-01
## rs150715708 0.000000e+00 0 1.000000e+00
## rs557714443 0.000000e+00 0 1.000000e+00
## rs3761373 7.881619e-01 1 3.746565e-01
## rs531410451 0.000000e+00 0 1.000000e+00
## rs140037718 0.000000e+00 0 1.000000e+00
## rs573737830 0.000000e+00 0 1.000000e+00
## rs542575245 0.000000e+00 0 1.000000e+00
## rs562776963 0.000000e+00 0 1.000000e+00
## rs531615786 0.000000e+00 0 1.000000e+00
## rs551432766 0.000000e+00 0 1.000000e+00
## rs191457025 0.000000e+00 0 1.000000e+00
## rs422761 1.030573e+00 1 3.100241e-01
## rs181091055 0.000000e+00 0 1.000000e+00
## rs186168224 0.000000e+00 0 1.000000e+00
## rs189570078 0.000000e+00 0 1.000000e+00
## rs549289994 0.000000e+00 0 1.000000e+00
## rs389001 1.030573e+00 1 3.100241e-01
## rs538139969 0.000000e+00 0 1.000000e+00
## rs557775656 0.000000e+00 0 1.000000e+00
## rs183398094 0.000000e+00 0 1.000000e+00
## rs146445857 0.000000e+00 0 1.000000e+00
## rs534034788 0.000000e+00 0 1.000000e+00
## rs2838043 2.041588e-01 1 6.513847e-01
## rs187648498 0.000000e+00 0 1.000000e+00
## rs149424945 0.000000e+00 0 1.000000e+00
## rs3761374 0.000000e+00 0 1.000000e+00
## rs576447019 0.000000e+00 0 1.000000e+00
## rs190673599 0.000000e+00 0 1.000000e+00
## rs34205539 4.429130e-01 1 5.057198e-01
## rs533552237 0.000000e+00 0 1.000000e+00
## rs77014365 0.000000e+00 0 1.000000e+00
## rs183722985 9.433138e-03 1 9.226276e-01
## rs146374360 0.000000e+00 0 1.000000e+00
## rs188198121 0.000000e+00 0 1.000000e+00
## rs529737187 0.000000e+00 0 1.000000e+00
## rs549593213 0.000000e+00 0 1.000000e+00
## rs571275590 0.000000e+00 0 1.000000e+00
## rs531884512 0.000000e+00 0 1.000000e+00
## rs78528160 6.064554e-02 1 8.054785e-01
## rs429524 1.597633e-01 1 6.893745e-01
## rs139661079 2.336398e-03 1 9.614482e-01
## rs553861191 0.000000e+00 0 1.000000e+00
## rs567783337 0.000000e+00 0 1.000000e+00
## rs144318842 2.142432e-02 1 8.836289e-01
## rs556371948 0.000000e+00 0 1.000000e+00
## rs576259678 0.000000e+00 0 1.000000e+00
## rs545333396 0.000000e+00 0 1.000000e+00
## rs558858943 0.000000e+00 0 1.000000e+00
## rs572348377 0.000000e+00 0 1.000000e+00
## rs368646465 0.000000e+00 0 1.000000e+00
## rs148719900 0.000000e+00 0 1.000000e+00
## rs142296178 0.000000e+00 0 1.000000e+00
## rs527416684 0.000000e+00 0 1.000000e+00
## rs562750287 0.000000e+00 0 1.000000e+00
## rs531754656 0.000000e+00 0 1.000000e+00
## rs551677463 0.000000e+00 0 1.000000e+00
## rs151248079 0.000000e+00 0 1.000000e+00
## rs527964326 0.000000e+00 0 1.000000e+00
## rs140408843 0.000000e+00 0 1.000000e+00
## rs567552202 0.000000e+00 0 1.000000e+00
## rs115266855 1.597633e-01 1 6.893745e-01
## rs556654990 0.000000e+00 0 1.000000e+00
## rs114027875 4.732632e+01 1 6.009970e-12
## rs538863731 0.000000e+00 0 1.000000e+00
## rs2187238 4.429130e-01 1 5.057198e-01
## rs192959767 0.000000e+00 0 1.000000e+00
## rs540964406 0.000000e+00 0 1.000000e+00
## rs554550449 0.000000e+00 0 1.000000e+00
## rs435877 2.857363e-02 1 8.657673e-01
## rs543518445 0.000000e+00 0 1.000000e+00
## rs73903404 1.211511e-01 1 7.277894e-01
## rs576510757 0.000000e+00 0 1.000000e+00
## rs545232590 0.000000e+00 0 1.000000e+00
## rs182838808 0.000000e+00 0 1.000000e+00
## rs186929947 0.000000e+00 0 1.000000e+00
## rs191345780 0.000000e+00 0 1.000000e+00
## rs368878191 0.000000e+00 0 1.000000e+00
## rs530009764 0.000000e+00 0 1.000000e+00
## rs371046741 0.000000e+00 0 1.000000e+00
## rs147620204 2.142432e-02 1 8.836289e-01
## rs539123545 0.000000e+00 0 1.000000e+00
## rs552447029 2.142432e-02 1 8.836289e-01
## rs565690211 0.000000e+00 0 1.000000e+00
## rs534924497 0.000000e+00 0 1.000000e+00
## rs535472251 6.629730e-01 1 4.155123e-01
## rs554530868 0.000000e+00 0 1.000000e+00
## rs9980693 2.336398e-03 1 9.614482e-01
## rs185103560 0.000000e+00 0 1.000000e+00
## rs73903405 3.109578e-01 1 5.770931e-01
## rs142119028 9.433138e-03 1 9.226276e-01
## rs545494331 0.000000e+00 0 1.000000e+00
## rs565142599 0.000000e+00 0 1.000000e+00
## rs572507505 0.000000e+00 0 1.000000e+00
## rs541726793 0.000000e+00 0 1.000000e+00
## rs189431131 0.000000e+00 0 1.000000e+00
## rs530402877 0.000000e+00 0 1.000000e+00
## rs568959921 0.000000e+00 0 1.000000e+00
## rs550238535 2.336398e-03 1 9.614482e-01
## rs563670115 0.000000e+00 0 1.000000e+00
## rs532655542 0.000000e+00 0 1.000000e+00
## rs145877432 2.142432e-02 1 8.836289e-01
## rs565879862 0.000000e+00 0 1.000000e+00
## rs138763189 2.336398e-03 1 9.614482e-01
## rs548341299 0.000000e+00 0 1.000000e+00
## rs377403765 0.000000e+00 0 1.000000e+00
## rs193019598 0.000000e+00 0 1.000000e+00
## rs149270377 3.844785e-02 1 8.445466e-01
## rs369948528 0.000000e+00 0 1.000000e+00
## rs143460343 3.844785e-02 1 8.445466e-01
## rs147977027 0.000000e+00 0 1.000000e+00
## rs141583878 0.000000e+00 0 1.000000e+00
## rs567302726 0.000000e+00 0 1.000000e+00
## rs555134126 0.000000e+00 0 1.000000e+00
## rs536147878 0.000000e+00 0 1.000000e+00
## rs544083065 0.000000e+00 0 1.000000e+00
## rs146385718 3.844785e-02 1 8.445466e-01
## rs532511763 0.000000e+00 0 1.000000e+00
## rs546337899 0.000000e+00 0 1.000000e+00
## rs559617171 0.000000e+00 0 1.000000e+00
## rs528568210 0.000000e+00 0 1.000000e+00
## rs139700775 0.000000e+00 0 1.000000e+00
## rs568244368 0.000000e+00 0 1.000000e+00
## rs530721239 0.000000e+00 0 1.000000e+00
## rs9984012 0.000000e+00 0 1.000000e+00
## rs142914234 3.844785e-02 1 8.445466e-01
## rs539019494 0.000000e+00 0 1.000000e+00
## rs552912416 0.000000e+00 0 1.000000e+00
## rs566485026 0.000000e+00 0 1.000000e+00
## rs535531364 9.433138e-03 1 9.226276e-01
## rs77599791 2.336398e-03 1 9.614482e-01
## rs373480225 0.000000e+00 0 1.000000e+00
## rs544141276 0.000000e+00 0 1.000000e+00
## rs557656729 0.000000e+00 0 1.000000e+00
## rs185296080 0.000000e+00 0 1.000000e+00
## rs545910743 0.000000e+00 0 1.000000e+00
## rs559691670 0.000000e+00 0 1.000000e+00
## rs420737 6.095175e-01 1 4.349695e-01
## rs542265380 0.000000e+00 0 1.000000e+00
## rs561557938 0.000000e+00 0 1.000000e+00
## rs530584676 0.000000e+00 0 1.000000e+00
## rs373611430 0.000000e+00 0 1.000000e+00
## rs570441750 0.000000e+00 0 1.000000e+00
## rs533043037 0.000000e+00 0 1.000000e+00
## rs146120690 0.000000e+00 0 1.000000e+00
## rs74749793 6.912000e-01 1 4.057565e-01
## rs181446489 0.000000e+00 0 1.000000e+00
## rs555405939 0.000000e+00 0 1.000000e+00
## rs568689321 0.000000e+00 0 1.000000e+00
## rs373052874 0.000000e+00 0 1.000000e+00
## rs376189701 0.000000e+00 0 1.000000e+00
## rs577650336 0.000000e+00 0 1.000000e+00
## rs186280517 0.000000e+00 0 1.000000e+00
## rs553097317 0.000000e+00 0 1.000000e+00
## rs573152591 2.336398e-03 1 9.614482e-01
## rs541977145 0.000000e+00 0 1.000000e+00
## rs561819256 0.000000e+00 0 1.000000e+00
## rs575313753 0.000000e+00 0 1.000000e+00
## rs544037190 0.000000e+00 0 1.000000e+00
## rs563930557 0.000000e+00 0 1.000000e+00
## rs116020930 2.336398e-03 1 9.614482e-01
## rs188257146 9.433138e-03 1 9.226276e-01
## rs7277080 1.185306e+00 1 2.762784e-01
## rs395584 6.095175e-01 1 4.349695e-01
## rs548678297 0.000000e+00 0 1.000000e+00
## rs8129713 2.041588e-01 1 6.513847e-01
## rs114848359 2.336398e-03 1 9.614482e-01
## rs551244742 0.000000e+00 0 1.000000e+00
## rs139829932 9.433138e-03 1 9.226276e-01
## rs374886738 0.000000e+00 0 1.000000e+00
## rs368411333 0.000000e+00 0 1.000000e+00
## rs573178525 0.000000e+00 0 1.000000e+00
## rs73357663 8.816327e-02 1 7.665253e-01
## rs555327583 0.000000e+00 0 1.000000e+00
## rs557494549 2.142432e-02 1 8.836289e-01
## rs145728087 0.000000e+00 0 1.000000e+00
## rs148499598 0.000000e+00 0 1.000000e+00
## rs563946938 0.000000e+00 0 1.000000e+00
## rs577731955 0.000000e+00 0 1.000000e+00
## rs180818774 2.336398e-03 1 9.614482e-01
## rs186418926 0.000000e+00 0 1.000000e+00
## rs74605993 8.816327e-02 1 7.665253e-01
## rs145297649 2.336398e-03 1 9.614482e-01
## rs75373173 0.000000e+00 0 1.000000e+00
## rs531439352 0.000000e+00 0 1.000000e+00
## rs147630170 0.000000e+00 0 1.000000e+00
## rs73357664 2.805211e+00 1 9.395850e-02
## rs548025657 0.000000e+00 0 1.000000e+00
## rs569429288 0.000000e+00 0 1.000000e+00
## rs567188556 0.000000e+00 0 1.000000e+00
## rs378616 2.272587e+00 1 1.316799e-01
## rs555390803 0.000000e+00 0 1.000000e+00
## rs183758144 0.000000e+00 0 1.000000e+00
## rs187037274 2.336398e-03 1 9.614482e-01
## rs372733383 0.000000e+00 0 1.000000e+00
## rs578015897 2.336398e-03 1 9.614482e-01
## rs115129572 8.816327e-02 1 7.665253e-01
## rs553673225 0.000000e+00 0 1.000000e+00
## rs573923452 0.000000e+00 0 1.000000e+00
## rs542898542 0.000000e+00 0 1.000000e+00
## rs562479511 0.000000e+00 0 1.000000e+00
## rs191526916 0.000000e+00 0 1.000000e+00
## rs545149163 0.000000e+00 0 1.000000e+00
## rs183385542 0.000000e+00 0 1.000000e+00
## rs527380714 0.000000e+00 0 1.000000e+00
## rs187143520 2.336398e-03 1 9.614482e-01
## rs371180501 2.142432e-02 1 8.836289e-01
## rs142261174 9.433138e-03 1 9.226276e-01
## rs529550811 2.336398e-03 1 9.614482e-01
## rs549230369 0.000000e+00 0 1.000000e+00
## rs28360562 4.429130e-01 1 5.057198e-01
## rs537876986 0.000000e+00 0 1.000000e+00
## rs557712450 0.000000e+00 0 1.000000e+00
## rs571750399 0.000000e+00 0 1.000000e+00
## rs537584838 0.000000e+00 0 1.000000e+00
## rs553934496 0.000000e+00 0 1.000000e+00
## rs573715927 0.000000e+00 0 1.000000e+00
## rs542961770 0.000000e+00 0 1.000000e+00
## rs144458055 1.211511e-01 3 9.891837e-01
## rs545056304 0.000000e+00 0 1.000000e+00
## rs151189718 0.000000e+00 0 1.000000e+00
## rs114641598 3.844785e-02 1 8.445466e-01
## rs191587330 0.000000e+00 0 1.000000e+00
## rs182330267 2.336398e-03 1 9.614482e-01
## rs140230703 0.000000e+00 0 1.000000e+00
## rs55704664 0.000000e+00 0 1.000000e+00
## rs75200570 7.881619e-01 1 3.746565e-01
## rs538465291 9.433138e-03 1 9.226276e-01
## rs551797163 0.000000e+00 0 1.000000e+00
## rs150314077 9.433138e-03 1 9.226276e-01
## rs534312135 0.000000e+00 0 1.000000e+00
## rs538226704 0.000000e+00 0 1.000000e+00
## rs186605876 0.000000e+00 0 1.000000e+00
## rs137929897 0.000000e+00 0 1.000000e+00
## rs556381293 2.336398e-03 1 9.614482e-01
## rs149021153 9.433138e-03 1 9.226276e-01
## rs538933388 0.000000e+00 0 1.000000e+00
## rs145900878 1.597633e-01 1 6.893745e-01
## rs138498737 0.000000e+00 0 1.000000e+00
## rs541139957 0.000000e+00 0 1.000000e+00
## rs560845787 0.000000e+00 0 1.000000e+00
## rs574176888 0.000000e+00 0 1.000000e+00
## rs57161767 8.816327e-02 1 7.665253e-01
## rs112467088 1.307071e+00 1 2.529256e-01
## rs531884709 0.000000e+00 0 1.000000e+00
## rs578020496 0.000000e+00 0 1.000000e+00
## rs565299100 0.000000e+00 0 1.000000e+00
## rs527658940 0.000000e+00 0 1.000000e+00
## rs185419671 0.000000e+00 0 1.000000e+00
## rs567794698 0.000000e+00 0 1.000000e+00
## rs73230088 2.336398e-03 1 9.614482e-01
## rs550049184 9.433138e-03 1 9.226276e-01
## rs553875292 2.336398e-03 1 9.614482e-01
## rs573736906 0.000000e+00 0 1.000000e+00
## rs558543141 0.000000e+00 0 1.000000e+00
## rs572279833 0.000000e+00 0 1.000000e+00
## rs534353584 0.000000e+00 0 1.000000e+00
## rs8126497 0.000000e+00 0 1.000000e+00
## rs574238899 0.000000e+00 0 1.000000e+00
## rs543066660 0.000000e+00 0 1.000000e+00
## rs562751769 0.000000e+00 0 1.000000e+00
## rs75603675 8.606343e-01 1 3.535613e-01
## rs200291871 0.000000e+00 0 1.000000e+00
## rs565468881 0.000000e+00 0 1.000000e+00
## rs543548317 9.433138e-03 1 9.226276e-01
## rs527970498 0.000000e+00 0 1.000000e+00
## rs541109823 0.000000e+00 0 1.000000e+00
## rs561266954 2.336398e-03 1 9.614482e-01
## rs542471574 0.000000e+00 0 1.000000e+00
## rs146132480 2.336398e-03 1 9.614482e-01
## rs569886879 0.000000e+00 0 1.000000e+00
## rs532232636 0.000000e+00 0 1.000000e+00
## rs4303794 8.606343e-01 1 3.535613e-01
## rs11088551 8.606343e-01 1 3.535613e-01
## rs66492316 8.606343e-01 1 3.535613e-01
## rs534913861 2.336398e-03 1 9.614482e-01
## rs4303795 8.606343e-01 1 3.535613e-01
## rs567676569 0.000000e+00 0 1.000000e+00
## rs5844077 9.917355e-02 1 7.528237e-01
## rs536763630 0.000000e+00 0 1.000000e+00
## rs76833541 0.000000e+00 0 1.000000e+00
## rs4283504 1.597633e-01 1 6.893745e-01
## rs545307151 0.000000e+00 0 1.000000e+00
## rs559098267 0.000000e+00 0 1.000000e+00
## rs572699613 0.000000e+00 0 1.000000e+00
## rs141478137 0.000000e+00 0 1.000000e+00
## rs561279116 0.000000e+00 0 1.000000e+00
## rs530084541 0.000000e+00 0 1.000000e+00
## rs543815381 0.000000e+00 0 1.000000e+00
## rs12481984 8.533594e-02 1 7.701925e-01
## rs532627816 0.000000e+00 0 1.000000e+00
## rs192854367 0.000000e+00 0 1.000000e+00
## rs565820920 0.000000e+00 0 1.000000e+00
## rs184146774 0.000000e+00 0 1.000000e+00
## rs188500480 0.000000e+00 0 1.000000e+00
## rs28707508 1.482555e-01 1 7.002080e-01
## rs552257429 6.268394e-01 1 4.285170e-01
## rs536677167 0.000000e+00 0 1.000000e+00
## rs553617111 0.000000e+00 0 1.000000e+00
## rs12626358 5.634047e-02 1 8.123764e-01
## rs570494222 0.000000e+00 0 1.000000e+00
## rs539520259 0.000000e+00 0 1.000000e+00
## rs202094412 0.000000e+00 0 1.000000e+00
## rs8128074 8.816327e-02 1 7.665253e-01
## rs555328523 0.000000e+00 0 1.000000e+00
## rs575168295 0.000000e+00 0 1.000000e+00
## rs150838246 0.000000e+00 0 1.000000e+00
## rs563693995 0.000000e+00 0 1.000000e+00
## rs139258152 2.336398e-03 1 9.614482e-01
## rs563980567 0.000000e+00 0 1.000000e+00
## rs181044447 0.000000e+00 0 1.000000e+00
## rs528259213 0.000000e+00 0 1.000000e+00
## rs547961708 0.000000e+00 0 1.000000e+00
## rs561682181 0.000000e+00 0 1.000000e+00
## rs56218846 4.340811e-02 1 8.349587e-01
## rs11281229 4.340811e-02 1 8.349587e-01
## rs550208798 0.000000e+00 0 1.000000e+00
## rs8127674 4.340811e-02 1 8.349587e-01
## rs115506893 2.142432e-02 1 8.836289e-01
## rs546643133 0.000000e+00 0 1.000000e+00
## rs566166020 0.000000e+00 0 1.000000e+00
## rs8132655 8.925620e-01 1 3.447840e-01
## rs8132786 2.024095e-01 1 6.527830e-01
## rs73230090 2.336398e-03 1 9.614482e-01
## rs9978557 5.187727e-01 1 4.713657e-01
## rs557087870 0.000000e+00 0 1.000000e+00
## rs577134436 0.000000e+00 0 1.000000e+00
## rs546001096 0.000000e+00 0 1.000000e+00
## rs9978803 1.354658e+00 1 2.444654e-01
## rs572987167 0.000000e+00 0 1.000000e+00
## rs541663911 0.000000e+00 0 1.000000e+00
## rs561542595 2.336398e-03 1 9.614482e-01
## rs55953905 9.433138e-03 1 9.226276e-01
## rs9981099 1.354658e+00 1 2.444654e-01
## rs563943938 0.000000e+00 0 1.000000e+00
## rs9981177 8.925620e-01 1 3.447840e-01
## rs147244020 0.000000e+00 0 1.000000e+00
## rs566600769 0.000000e+00 0 1.000000e+00
## rs12627186 5.869029e-01 1 4.436193e-01
## rs548801187 0.000000e+00 0 1.000000e+00
## rs568628088 0.000000e+00 0 1.000000e+00
## rs537756467 0.000000e+00 0 1.000000e+00
## rs182546346 0.000000e+00 0 1.000000e+00
## rs577444816 0.000000e+00 0 1.000000e+00
## rs9981463 1.354658e+00 1 2.444654e-01
## rs140638214 0.000000e+00 0 1.000000e+00
## rs573048129 0.000000e+00 0 1.000000e+00
## rs9981563 1.354658e+00 1 2.444654e-01
## rs9981570 1.354658e+00 1 2.444654e-01
## rs7277527 1.929709e-01 1 6.604554e-01
## rs187111053 0.000000e+00 0 1.000000e+00
## rs147144661 6.064554e-02 1 8.054785e-01
## rs73357672 9.433138e-03 1 9.226276e-01
## rs192075489 0.000000e+00 0 1.000000e+00
## rs183401776 0.000000e+00 0 1.000000e+00
## rs540895190 0.000000e+00 0 1.000000e+00
## rs9981973 1.354658e+00 1 2.444654e-01
## rs529005816 0.000000e+00 0 1.000000e+00
## rs549064372 0.000000e+00 0 1.000000e+00
## rs529913299 0.000000e+00 0 1.000000e+00
## rs568691585 0.000000e+00 0 1.000000e+00
## rs115774689 1.211511e-01 1 7.277894e-01
## rs550879629 0.000000e+00 0 1.000000e+00
## rs571165092 0.000000e+00 0 1.000000e+00
## rs539741673 0.000000e+00 0 1.000000e+00
## rs553335412 0.000000e+00 0 1.000000e+00
## rs144500278 0.000000e+00 0 1.000000e+00
## rs185939056 0.000000e+00 0 1.000000e+00
## rs145572528 0.000000e+00 0 1.000000e+00
## rs575315760 0.000000e+00 0 1.000000e+00
## rs74849454 6.912000e-01 1 4.057565e-01
## rs8129192 8.925620e-01 1 3.447840e-01
## rs190253467 0.000000e+00 0 1.000000e+00
## rs9984523 5.187727e-01 1 4.713657e-01
## rs560205728 0.000000e+00 0 1.000000e+00
## rs528866229 0.000000e+00 0 1.000000e+00
## rs542622474 0.000000e+00 0 1.000000e+00
## rs562604402 0.000000e+00 0 1.000000e+00
## rs8133719 1.354658e+00 1 2.444654e-01
## rs551204780 0.000000e+00 0 1.000000e+00
## rs182567672 0.000000e+00 0 1.000000e+00
## rs539186653 0.000000e+00 0 1.000000e+00
## rs554512076 0.000000e+00 0 1.000000e+00
## rs566551507 0.000000e+00 0 1.000000e+00
## rs142565367 3.844785e-02 1 8.445466e-01
## rs555327867 0.000000e+00 0 1.000000e+00
## rs117520130 0.000000e+00 0 1.000000e+00
## rs537904115 0.000000e+00 0 1.000000e+00
## rs116880445 0.000000e+00 0 1.000000e+00
## rs116368509 1.211511e-01 1 7.277894e-01
## rs534116368 0.000000e+00 0 1.000000e+00
## rs554136103 0.000000e+00 0 1.000000e+00
## rs115356706 1.211511e-01 1 7.277894e-01
## rs13046567 2.022382e+00 1 1.549958e-01
## rs562465933 0.000000e+00 0 1.000000e+00
## rs11909349 2.966147e-01 1 5.860125e-01
## rs186851415 9.433138e-03 1 9.226276e-01
## rs558652539 0.000000e+00 0 1.000000e+00
## rs150947188 0.000000e+00 0 1.000000e+00
## rs546770118 0.000000e+00 0 1.000000e+00
## rs374815822 0.000000e+00 0 1.000000e+00
## rs570586955 0.000000e+00 0 1.000000e+00
## rs144630231 2.336398e-03 1 9.614482e-01
## rs138518282 1.211511e-01 1 7.277894e-01
## rs538350112 0.000000e+00 0 1.000000e+00
## rs575933876 9.433138e-03 1 9.226276e-01
## rs13048893 1.354658e+00 1 2.444654e-01
## rs145145381 9.433138e-03 1 9.226276e-01
## rs553999473 0.000000e+00 0 1.000000e+00
## rs576001583 0.000000e+00 0 1.000000e+00
## rs112719709 2.336398e-03 1 9.614482e-01
## rs556204609 0.000000e+00 0 1.000000e+00
## rs138904458 9.433138e-03 1 9.226276e-01
## rs545113252 0.000000e+00 0 1.000000e+00
## rs564861368 0.000000e+00 0 1.000000e+00
## rs572093544 0.000000e+00 0 1.000000e+00
## rs184156179 0.000000e+00 0 1.000000e+00
## rs560264501 0.000000e+00 0 1.000000e+00
## rs529550773 0.000000e+00 0 1.000000e+00
## rs549228416 0.000000e+00 0 1.000000e+00
## rs562920187 0.000000e+00 0 1.000000e+00
## rs555250676 2.336398e-03 1 9.614482e-01
## rs551736035 0.000000e+00 0 1.000000e+00
## rs13050241 1.354658e+00 1 2.444654e-01
## rs559113271 0.000000e+00 0 1.000000e+00
## rs144780284 0.000000e+00 0 1.000000e+00
## rs373939423 0.000000e+00 0 1.000000e+00
## rs536567653 0.000000e+00 0 1.000000e+00
## rs556158219 0.000000e+00 0 1.000000e+00
## rs569913599 0.000000e+00 0 1.000000e+00
## rs538292644 0.000000e+00 0 1.000000e+00
## rs558518683 2.336398e-03 1 9.614482e-01
## rs532503621 0.000000e+00 0 1.000000e+00
## rs540682335 0.000000e+00 0 1.000000e+00
## rs188751108 0.000000e+00 0 1.000000e+00
## rs573865390 0.000000e+00 0 1.000000e+00
## rs8129232 1.354658e+00 1 2.444654e-01
## rs148125307 0.000000e+00 0 1.000000e+00
## rs4818243 0.000000e+00 0 1.000000e+00
## rs545118901 0.000000e+00 0 1.000000e+00
## rs143900931 0.000000e+00 0 1.000000e+00
## rs147254430 0.000000e+00 0 1.000000e+00
## rs140678733 1.211511e-01 1 7.277894e-01
## rs4339034 1.354658e+00 1 2.444654e-01
## rs529945723 0.000000e+00 0 1.000000e+00
## rs6517670 8.925620e-01 1 3.447840e-01
## rs193193107 0.000000e+00 0 1.000000e+00
## rs538923979 0.000000e+00 0 1.000000e+00
## rs558275427 0.000000e+00 0 1.000000e+00
## rs187282381 0.000000e+00 0 1.000000e+00
## rs4437046 1.354658e+00 1 2.444654e-01
## rs554269982 0.000000e+00 0 1.000000e+00
## rs111419457 6.064554e-02 1 8.054785e-01
## rs142914934 2.336398e-03 1 9.614482e-01
## rs556713596 0.000000e+00 0 1.000000e+00
## rs4635604 2.325824e+00 1 1.272428e-01
## rs545642867 0.000000e+00 0 1.000000e+00
## rs565153477 0.000000e+00 0 1.000000e+00
## rs527659955 0.000000e+00 0 1.000000e+00
## rs28751065 0.000000e+00 0 1.000000e+00
## rs4441561 1.354658e+00 1 2.444654e-01
## rs530212386 0.000000e+00 0 1.000000e+00
## rs374197073 0.000000e+00 0 1.000000e+00
## rs57859815 2.041588e-01 1 6.513847e-01
## rs532317051 0.000000e+00 0 1.000000e+00
## rs552417594 0.000000e+00 0 1.000000e+00
## rs191752543 0.000000e+00 0 1.000000e+00
## rs534348818 0.000000e+00 0 1.000000e+00
## rs139152704 2.142432e-02 1 8.836289e-01
## rs9979885 0.000000e+00 0 1.000000e+00
## rs6517671 7.500000e-01 1 3.864762e-01
## rs8134056 7.500000e-01 1 3.864762e-01
## rs8133813 7.500000e-01 1 3.864762e-01
## rs9789878 5.869029e-01 1 4.436193e-01
## rs8133685 7.500000e-01 1 3.864762e-01
## rs8134861 7.500000e-01 1 3.864762e-01
## rs541479740 0.000000e+00 0 1.000000e+00
## rs149941881 0.000000e+00 0 1.000000e+00
## rs574991676 0.000000e+00 0 1.000000e+00
## rs543862636 0.000000e+00 0 1.000000e+00
## rs563357000 0.000000e+00 0 1.000000e+00
## rs542029547 0.000000e+00 0 1.000000e+00
## rs552069008 0.000000e+00 0 1.000000e+00
## rs180846756 0.000000e+00 0 1.000000e+00
## rs528365377 0.000000e+00 0 1.000000e+00
## rs547786160 0.000000e+00 0 1.000000e+00
## rs9982595 2.545009e-01 1 6.139235e-01
## rs114933510 9.433138e-03 1 9.226276e-01
## rs570432763 0.000000e+00 0 1.000000e+00
## rs9982757 1.787771e+00 1 1.811980e-01
## rs117700839 0.000000e+00 0 1.000000e+00
## rs185854865 0.000000e+00 0 1.000000e+00
## rs9982761 1.354658e+00 1 2.444654e-01
## rs190429798 0.000000e+00 0 1.000000e+00
## rs574745113 0.000000e+00 0 1.000000e+00
## rs370266410 0.000000e+00 0 1.000000e+00
## rs181769907 0.000000e+00 0 1.000000e+00
## rs577255863 0.000000e+00 0 1.000000e+00
## rs186175806 0.000000e+00 0 1.000000e+00
## rs190501444 0.000000e+00 0 1.000000e+00
## rs2210282 1.354658e+00 1 2.444654e-01
## rs548047761 0.000000e+00 0 1.000000e+00
## rs144999434 0.000000e+00 0 1.000000e+00
## rs139881855 2.336398e-03 1 9.614482e-01
## rs550291786 0.000000e+00 0 1.000000e+00
## rs570569709 0.000000e+00 0 1.000000e+00
## rs143256411 0.000000e+00 0 1.000000e+00
## rs181383657 0.000000e+00 0 1.000000e+00
## rs9983001 0.000000e+00 0 1.000000e+00
## rs535359170 0.000000e+00 0 1.000000e+00
## rs531373761 0.000000e+00 0 1.000000e+00
## rs146670856 0.000000e+00 0 1.000000e+00
## rs185608459 0.000000e+00 0 1.000000e+00
## rs557073855 0.000000e+00 0 1.000000e+00
## rs2097001 1.354658e+00 1 2.444654e-01
## rs576787597 0.000000e+00 0 1.000000e+00
## rs552764896 0.000000e+00 0 1.000000e+00
## rs2097002 1.354658e+00 1 2.444654e-01
## rs76041556 5.869029e-01 1 4.436193e-01
## rs562050387 0.000000e+00 0 1.000000e+00
## rs191201574 0.000000e+00 0 1.000000e+00
## rs2097003 1.354658e+00 1 2.444654e-01
## rs563883165 0.000000e+00 0 1.000000e+00
## rs2097004 1.058968e+00 1 3.034506e-01
## rs563489305 0.000000e+00 0 1.000000e+00
## rs143397179 9.433138e-03 1 9.226276e-01
## rs77947664 6.064554e-02 1 8.054785e-01
## rs528728595 0.000000e+00 0 1.000000e+00
## rs548789098 0.000000e+00 0 1.000000e+00
## rs568781224 0.000000e+00 0 1.000000e+00
## rs537298231 0.000000e+00 0 1.000000e+00
## rs557333375 0.000000e+00 0 1.000000e+00
## rs8130911 1.354658e+00 1 2.444654e-01
## rs539567413 0.000000e+00 0 1.000000e+00
## rs553080083 0.000000e+00 0 1.000000e+00
## rs572987230 0.000000e+00 0 1.000000e+00
## rs541664199 0.000000e+00 0 1.000000e+00
## rs11414078 1.354658e+00 1 2.444654e-01
## rs555270521 0.000000e+00 0 1.000000e+00
## rs575499123 0.000000e+00 0 1.000000e+00
## rs140600701 0.000000e+00 0 1.000000e+00
## rs182532730 0.000000e+00 0 1.000000e+00
## rs527393608 0.000000e+00 0 1.000000e+00
## rs540215789 0.000000e+00 0 1.000000e+00
## rs6517672 3.685077e-01 1 5.438188e-01
## rs187705808 2.336398e-03 1 9.614482e-01
## rs62219272 3.109578e-01 1 5.770931e-01
## rs138465570 0.000000e+00 0 1.000000e+00
## rs531250546 0.000000e+00 0 1.000000e+00
## rs192040797 0.000000e+00 0 1.000000e+00
## rs114970614 6.064554e-02 1 8.054785e-01
## rs7276907 3.507074e-01 1 5.537130e-01
## rs546580569 0.000000e+00 0 1.000000e+00
## rs566672751 0.000000e+00 0 1.000000e+00
## rs535512482 0.000000e+00 0 1.000000e+00
## rs555550283 0.000000e+00 0 1.000000e+00
## rs575314488 0.000000e+00 0 1.000000e+00
## rs538019747 0.000000e+00 0 1.000000e+00
## rs66953820 1.208040e-03 1 9.722736e-01
## rs577990422 0.000000e+00 0 1.000000e+00
## rs146495790 0.000000e+00 0 1.000000e+00
## rs76043846 1.268920e-01 1 7.216765e-01
## rs114020457 8.816327e-02 1 7.665253e-01
## rs542323159 2.336398e-03 1 9.614482e-01
## rs73230101 7.057946e-02 1 7.904950e-01
## rs531040705 0.000000e+00 0 1.000000e+00
## rs551138646 2.142432e-02 1 8.836289e-01
## rs564481444 0.000000e+00 0 1.000000e+00
## rs533443329 0.000000e+00 0 1.000000e+00
## rs139037537 3.737024e-01 1 5.409930e-01
## rs377096823 0.000000e+00 0 1.000000e+00
## rs112714648 1.011790e-02 1 9.198777e-01
## rs10211909 5.266526e-02 1 8.184888e-01
## rs569316663 0.000000e+00 0 1.000000e+00
## rs537683306 0.000000e+00 0 1.000000e+00
## rs557926869 0.000000e+00 0 1.000000e+00
## rs187495341 0.000000e+00 0 1.000000e+00
## rs192013336 0.000000e+00 0 1.000000e+00
## rs553817708 0.000000e+00 0 1.000000e+00
## rs572596202 0.000000e+00 0 1.000000e+00
## rs10211872 1.597633e-01 1 6.893745e-01
## rs542586667 0.000000e+00 0 1.000000e+00
## rs10211876 1.597633e-01 1 6.893745e-01
## rs576078506 0.000000e+00 0 1.000000e+00
## rs544527140 0.000000e+00 0 1.000000e+00
## rs200265359 2.442196e-02 1 8.758161e-01
## rs7282202 8.175472e+00 1 4.246063e-03
## rs533457313 0.000000e+00 0 1.000000e+00
## rs57197209 1.138472e+01 1 7.405075e-04
## rs7281137 1.597633e-01 1 6.893745e-01
## rs529302729 0.000000e+00 0 1.000000e+00
## rs35304539 1.597633e-01 1 6.893745e-01
## rs549346140 0.000000e+00 0 1.000000e+00
## rs569297723 0.000000e+00 0 1.000000e+00
## rs567232784 0.000000e+00 0 1.000000e+00
## rs7281696 1.597633e-01 1 6.893745e-01
## rs11700751 0.000000e+00 0 1.000000e+00
## rs534102071 0.000000e+00 0 1.000000e+00
## rs187943964 2.336398e-03 1 9.614482e-01
## rs567551843 0.000000e+00 0 1.000000e+00
## rs143295085 2.142432e-02 1 8.836289e-01
## rs146681963 0.000000e+00 0 1.000000e+00
## rs139187042 2.336398e-03 1 9.614482e-01
## rs145677425 9.189000e-01 1 3.377640e-01
## rs3037876 5.985027e-01 1 4.391499e-01
## rs2410431 3.431863e-01 1 5.579958e-01
## rs71285081 9.835236e-01 1 3.213304e-01
## rs149890385 3.844785e-02 1 8.445466e-01
## rs192835740 0.000000e+00 0 1.000000e+00
## rs2210283 1.597633e-01 1 6.893745e-01
## rs529586385 0.000000e+00 0 1.000000e+00
## rs184368284 0.000000e+00 0 1.000000e+00
## rs562858847 0.000000e+00 0 1.000000e+00
## rs531748551 2.336398e-03 1 9.614482e-01
## rs11701348 1.268920e-01 1 7.216765e-01
## rs11701382 1.268920e-01 1 7.216765e-01
## rs527497781 0.000000e+00 0 1.000000e+00
## rs374098316 0.000000e+00 0 1.000000e+00
## rs549672486 0.000000e+00 0 1.000000e+00
## rs536542935 0.000000e+00 0 1.000000e+00
## rs144880274 0.000000e+00 0 1.000000e+00
## rs569616563 0.000000e+00 0 1.000000e+00
## rs538371853 0.000000e+00 0 1.000000e+00
## rs549245809 0.000000e+00 0 1.000000e+00
## rs189341647 0.000000e+00 0 1.000000e+00
## rs73903407 1.211511e-01 1 7.277894e-01
## rs149015327 0.000000e+00 0 1.000000e+00
## rs181965615 0.000000e+00 0 1.000000e+00
## rs574295783 0.000000e+00 0 1.000000e+00
## rs543291662 0.000000e+00 0 1.000000e+00
## rs184851438 0.000000e+00 0 1.000000e+00
## rs189179693 0.000000e+00 0 1.000000e+00
## rs545386863 0.000000e+00 0 1.000000e+00
## rs565387759 0.000000e+00 0 1.000000e+00
## rs56134964 1.211511e-01 1 7.277894e-01
## rs12482673 1.597633e-01 1 6.893745e-01
## rs561125651 0.000000e+00 0 1.000000e+00
## rs529976954 0.000000e+00 0 1.000000e+00
## rs550025541 2.336398e-03 1 9.614482e-01
## rs181582940 8.816327e-02 1 7.665253e-01
## rs538259425 0.000000e+00 0 1.000000e+00
## rs7281569 8.483839e-02 1 7.708446e-01
## rs565487725 0.000000e+00 0 1.000000e+00
## rs538562119 0.000000e+00 0 1.000000e+00
## rs554107686 0.000000e+00 0 1.000000e+00
## rs186717374 0.000000e+00 0 1.000000e+00
## rs143799618 0.000000e+00 0 1.000000e+00
## rs7281987 3.685077e-01 1 5.438188e-01
## rs7282236 5.537449e-01 1 4.567916e-01
## rs545422998 0.000000e+00 0 1.000000e+00
## rs538165269 0.000000e+00 0 1.000000e+00
## rs80321213 0.000000e+00 0 1.000000e+00
## rs572611425 0.000000e+00 0 1.000000e+00
## rs541520587 0.000000e+00 0 1.000000e+00
## rs555164860;rs34666868 1.417171e+00 1 2.338697e-01
## rs73231904 0.000000e+00 0 1.000000e+00
## rs191214512 0.000000e+00 0 1.000000e+00
## rs549673724 0.000000e+00 0 1.000000e+00
## rs563606222 0.000000e+00 0 1.000000e+00
## rs532273790 0.000000e+00 0 1.000000e+00
## rs574884040 0.000000e+00 0 1.000000e+00
## rs73903408 1.211511e-01 1 7.277894e-01
## rs534441083 0.000000e+00 0 1.000000e+00
## rs548288699 0.000000e+00 0 1.000000e+00
## rs139183850 0.000000e+00 0 1.000000e+00
## rs6517673 9.711553e-02 1 7.553190e-01
## rs556434796 0.000000e+00 0 1.000000e+00
## rs570408392 0.000000e+00 0 1.000000e+00
## rs141467809 3.844785e-02 1 8.445466e-01
## rs559025205 0.000000e+00 0 1.000000e+00
## rs572447702 0.000000e+00 0 1.000000e+00
## rs77325194 0.000000e+00 0 1.000000e+00
## rs555058445 0.000000e+00 0 1.000000e+00
## rs73231906 4.937457e-02 1 8.241552e-01
## rs543652613 0.000000e+00 0 1.000000e+00
## rs73903409 1.211511e-01 1 7.277894e-01
## rs181093220 0.000000e+00 0 1.000000e+00
## rs138139129 0.000000e+00 0 1.000000e+00
## rs559195470 0.000000e+00 0 1.000000e+00
## rs571854174 0.000000e+00 0 1.000000e+00
## rs547836968 0.000000e+00 0 1.000000e+00
## rs189662163 0.000000e+00 0 1.000000e+00
## rs530517698 0.000000e+00 0 1.000000e+00
## rs550390247 0.000000e+00 0 1.000000e+00
## rs570247596 0.000000e+00 0 1.000000e+00
## rs149541145 0.000000e+00 0 1.000000e+00
## rs35247320 2.336398e-03 1 9.614482e-01
## rs146023440 2.336398e-03 1 9.614482e-01
## rs535024632 0.000000e+00 0 1.000000e+00
## rs554617100 0.000000e+00 0 1.000000e+00
## rs114264164 6.064554e-02 1 8.054785e-01
## rs543537596 0.000000e+00 0 1.000000e+00
## rs556768034 0.000000e+00 0 1.000000e+00
## rs576650666 0.000000e+00 0 1.000000e+00
## rs11908828 7.974953e-02 1 7.776371e-01
## rs141787670 0.000000e+00 0 1.000000e+00
## rs564388907 0.000000e+00 0 1.000000e+00
## rs146241351 9.433138e-03 1 9.226276e-01
## rs144274178 3.737024e-01 1 5.409930e-01
## rs561517039 3.737024e-01 1 5.409930e-01
## rs530663805 0.000000e+00 0 1.000000e+00
## rs550381527 0.000000e+00 0 1.000000e+00
## rs570532512 0.000000e+00 0 1.000000e+00
## rs139229175 2.336398e-03 1 9.614482e-01
## rs142553463 2.336398e-03 1 9.614482e-01
## rs566376925 0.000000e+00 0 1.000000e+00
## rs73357702 3.658352e-01 1 5.452832e-01
## rs548600814 0.000000e+00 0 1.000000e+00
## rs568344205 0.000000e+00 0 1.000000e+00
## rs537426013 0.000000e+00 0 1.000000e+00
## rs183394810 0.000000e+00 0 1.000000e+00
## rs186254487 0.000000e+00 0 1.000000e+00
## rs539252974 0.000000e+00 0 1.000000e+00
## rs190661548 0.000000e+00 0 1.000000e+00
## rs150513522 0.000000e+00 0 1.000000e+00
## rs541742281 0.000000e+00 0 1.000000e+00
## rs561750829 0.000000e+00 0 1.000000e+00
## rs139447364 0.000000e+00 0 1.000000e+00
## rs544190588 0.000000e+00 0 1.000000e+00
## rs182915405 6.064554e-02 1 8.054785e-01
## rs150062514 0.000000e+00 0 1.000000e+00
## rs534024498 0.000000e+00 0 1.000000e+00
## rs559973220 0.000000e+00 0 1.000000e+00
## rs371744816 0.000000e+00 0 1.000000e+00
## rs548625038 0.000000e+00 0 1.000000e+00
## rs114644785 2.142432e-02 1 8.836289e-01
## rs8134378 1.147395e-01 1 7.348115e-01
## rs550514293 0.000000e+00 0 1.000000e+00
## rs56095453 1.273219e-01 1 7.212252e-01
## rs188479461 0.000000e+00 0 1.000000e+00
## rs75115548 1.597633e-01 1 6.893745e-01
## rs192746461 0.000000e+00 0 1.000000e+00
## rs8134657 1.147395e-01 1 7.348115e-01
## rs555388366 0.000000e+00 0 1.000000e+00
## rs575453095 0.000000e+00 0 1.000000e+00
## rs183127646 0.000000e+00 0 1.000000e+00
## rs187713047 3.844785e-02 1 8.445466e-01
## rs577849576 0.000000e+00 0 1.000000e+00
## rs62219275 0.000000e+00 0 1.000000e+00
## rs191337624 0.000000e+00 0 1.000000e+00
## rs528696415 0.000000e+00 0 1.000000e+00
## rs147734689 2.336398e-03 1 9.614482e-01
## rs562247792 0.000000e+00 0 1.000000e+00
## rs531106884 2.336398e-03 1 9.614482e-01
## rs550602035 0.000000e+00 0 1.000000e+00
## rs570496645 0.000000e+00 0 1.000000e+00
## rs533278510 0.000000e+00 0 1.000000e+00
## rs546709342 0.000000e+00 0 1.000000e+00
## rs142630406 0.000000e+00 0 1.000000e+00
## rs185091034 0.000000e+00 0 1.000000e+00
## rs555489343 0.000000e+00 0 1.000000e+00
## rs13050216 0.000000e+00 0 1.000000e+00
## rs35616234 5.537449e-01 1 4.567916e-01
## rs538126737 2.336398e-03 1 9.614482e-01
## rs35221358 0.000000e+00 0 1.000000e+00
## rs147924692 2.336398e-03 1 9.614482e-01
## rs148964290 6.064554e-02 1 8.054785e-01
## rs4818244 0.000000e+00 0 1.000000e+00
## rs143692782 0.000000e+00 0 1.000000e+00
## rs573635382 0.000000e+00 0 1.000000e+00
## rs189418914 0.000000e+00 0 1.000000e+00
## rs11088552 4.573750e-02 1 8.306537e-01
## rs531142275 0.000000e+00 0 1.000000e+00
## rs140912147 0.000000e+00 0 1.000000e+00
## rs192376344 0.000000e+00 0 1.000000e+00
## rs533169706 0.000000e+00 0 1.000000e+00
## rs546930405 0.000000e+00 0 1.000000e+00
## rs566853226 0.000000e+00 0 1.000000e+00
## rs144895242 0.000000e+00 0 1.000000e+00
## rs11088553 3.658352e-01 1 5.452832e-01
## rs569015070 0.000000e+00 0 1.000000e+00
## rs373900580 0.000000e+00 0 1.000000e+00
## rs558100984 0.000000e+00 0 1.000000e+00
## rs185173399 0.000000e+00 0 1.000000e+00
## rs190114929 0.000000e+00 0 1.000000e+00
## rs553727290 0.000000e+00 0 1.000000e+00
## rs573695555 0.000000e+00 0 1.000000e+00
## rs542693019 0.000000e+00 0 1.000000e+00
## rs555710972 0.000000e+00 0 1.000000e+00
## rs575455297 0.000000e+00 0 1.000000e+00
## rs577318123 0.000000e+00 0 1.000000e+00
## rs11910588 8.130363e-01 1 3.672241e-01
## rs543681865 0.000000e+00 0 1.000000e+00
## rs562387304 0.000000e+00 0 1.000000e+00
## rs573651044 0.000000e+00 0 1.000000e+00
## rs529400107 0.000000e+00 0 1.000000e+00
## rs74504161 3.737024e-01 1 5.409930e-01
## rs181420525 2.336398e-03 1 9.614482e-01
## rs371636071 0.000000e+00 0 1.000000e+00
## rs58282962 1.510079e-01 1 6.975740e-01
## rs186433432 0.000000e+00 0 1.000000e+00
## rs534085241 2.336398e-03 1 9.614482e-01
## rs141824966 0.000000e+00 0 1.000000e+00
## rs567231355 0.000000e+00 0 1.000000e+00
## rs535844612 0.000000e+00 0 1.000000e+00
## rs188347907 0.000000e+00 0 1.000000e+00
## rs575689705 0.000000e+00 0 1.000000e+00
## rs538049916 0.000000e+00 0 1.000000e+00
## rs557925898 0.000000e+00 0 1.000000e+00
## rs577955731 0.000000e+00 0 1.000000e+00
## rs540902787 0.000000e+00 0 1.000000e+00
## rs560597275 0.000000e+00 0 1.000000e+00
## rs114451099 2.142432e-02 1 8.836289e-01
## rs180808723 0.000000e+00 0 1.000000e+00
## rs115361224 2.336398e-03 1 9.614482e-01
## rs531839145 0.000000e+00 0 1.000000e+00
## rs551642864 0.000000e+00 0 1.000000e+00
## rs565059986 0.000000e+00 0 1.000000e+00
## rs527619469 0.000000e+00 0 1.000000e+00
## rs146252482 0.000000e+00 0 1.000000e+00
## rs567413133 0.000000e+00 0 1.000000e+00
## rs138015397 0.000000e+00 0 1.000000e+00
## rs375750497 0.000000e+00 0 1.000000e+00
## rs569456729 0.000000e+00 0 1.000000e+00
## rs4818245 2.404855e-01 1 6.238557e-01
## rs557962349 0.000000e+00 0 1.000000e+00
## rs149489375 5.187727e-01 1 4.713657e-01
## rs143911765 6.014703e-01 1 4.380176e-01
## rs554458606 0.000000e+00 0 1.000000e+00
## rs376752665 0.000000e+00 0 1.000000e+00
## rs543057150 0.000000e+00 0 1.000000e+00
## rs562861772 0.000000e+00 0 1.000000e+00
## rs113192256 0.000000e+00 0 1.000000e+00
## rs11911639 0.000000e+00 0 1.000000e+00
## rs190908789 0.000000e+00 0 1.000000e+00
## rs11911661 1.165795e+00 1 2.802669e-01
## rs62219277 6.014703e-01 1 4.380176e-01
## rs561239256 2.336398e-03 1 9.614482e-01
## rs529939909 0.000000e+00 0 1.000000e+00
## rs2838045 5.187727e-01 1 4.713657e-01
## rs569293185 0.000000e+00 0 1.000000e+00
## rs532102372 0.000000e+00 0 1.000000e+00
## rs551872615 0.000000e+00 0 1.000000e+00
## rs565599424 0.000000e+00 0 1.000000e+00
## rs534411308 0.000000e+00 0 1.000000e+00
## rs554071332 0.000000e+00 0 1.000000e+00
## rs75579634 3.844785e-02 1 8.445466e-01
## rs183944200 2.336398e-03 1 9.614482e-01
## rs11702754 3.658352e-01 1 5.452832e-01
## rs576528439 0.000000e+00 0 1.000000e+00
## rs140537307 3.109578e-01 1 5.770931e-01
## rs533054515 0.000000e+00 0 1.000000e+00
## rs73231912 6.033647e-01 1 4.372971e-01
## rs73231913 6.033647e-01 1 4.372971e-01
## rs560913568 0.000000e+00 0 1.000000e+00
## rs11700636 1.130450e-01 1 7.367037e-01
## rs186810667 0.000000e+00 0 1.000000e+00
## rs150011210 0.000000e+00 0 1.000000e+00
## rs531991508 0.000000e+00 0 1.000000e+00
## rs145227429 0.000000e+00 0 1.000000e+00
## rs12151964 2.336398e-03 1 9.614482e-01
## rs528166702 0.000000e+00 0 1.000000e+00
## rs548004905 0.000000e+00 0 1.000000e+00
## rs567738745 0.000000e+00 0 1.000000e+00
## rs371450260 0.000000e+00 0 1.000000e+00
## rs556910361 0.000000e+00 0 1.000000e+00
## rs570321028 0.000000e+00 0 1.000000e+00
## rs539037582 0.000000e+00 0 1.000000e+00
## rs149163386 9.433138e-03 1 9.226276e-01
## rs145013758 0.000000e+00 0 1.000000e+00
## rs191770613 0.000000e+00 0 1.000000e+00
## rs138803889 0.000000e+00 0 1.000000e+00
## rs183116984 0.000000e+00 0 1.000000e+00
## rs140419297 0.000000e+00 0 1.000000e+00
## rs188249333 2.336398e-03 1 9.614482e-01
## rs548879701 8.816327e-02 1 7.665253e-01
## rs576726249 0.000000e+00 0 1.000000e+00
## rs191576885 0.000000e+00 0 1.000000e+00
## rs568364386 0.000000e+00 0 1.000000e+00
## rs528217706 2.336398e-03 1 9.614482e-01
## rs548004327 0.000000e+00 0 1.000000e+00
## rs58987663 5.325570e-02 1 8.174920e-01
## rs561789686 0.000000e+00 0 1.000000e+00
## rs182320551 0.000000e+00 0 1.000000e+00
## rs145529392 0.000000e+00 0 1.000000e+00
## rs570408544 0.000000e+00 0 1.000000e+00
## rs186603866 0.000000e+00 0 1.000000e+00
## rs115668254 1.597633e-01 1 6.893745e-01
## rs565911015 0.000000e+00 0 1.000000e+00
## rs192580203 0.000000e+00 0 1.000000e+00
## rs185634591 0.000000e+00 0 1.000000e+00
## rs574547899 0.000000e+00 0 1.000000e+00
## rs148034673 0.000000e+00 0 1.000000e+00
## rs7280868 3.113611e-01 1 5.768461e-01
## rs576763786 0.000000e+00 0 1.000000e+00
## rs545757847 2.336398e-03 1 9.614482e-01
## rs559569997 0.000000e+00 0 1.000000e+00
## rs190106715 0.000000e+00 0 1.000000e+00
## rs541767767 0.000000e+00 0 1.000000e+00
## rs2838046 3.113611e-01 1 5.768461e-01
## rs530618388 0.000000e+00 0 1.000000e+00
## rs550436579 0.000000e+00 0 1.000000e+00
## rs563853849 0.000000e+00 0 1.000000e+00
## rs532881056 0.000000e+00 0 1.000000e+00
## rs552614787 0.000000e+00 0 1.000000e+00
## rs77013654 0.000000e+00 0 1.000000e+00
## rs534713766 0.000000e+00 0 1.000000e+00
## rs548561732 0.000000e+00 0 1.000000e+00
## rs568258055 0.000000e+00 0 1.000000e+00
## rs537135183 0.000000e+00 0 1.000000e+00
## rs556766125 0.000000e+00 0 1.000000e+00
## rs576650492 0.000000e+00 0 1.000000e+00
## rs370678258 0.000000e+00 0 1.000000e+00
## rs116635332 2.336398e-03 1 9.614482e-01
## rs141713948 0.000000e+00 0 1.000000e+00
## rs114127218 2.336398e-03 1 9.614482e-01
## rs193154554 0.000000e+00 0 1.000000e+00
## rs575445630 0.000000e+00 0 1.000000e+00
## rs544300558 0.000000e+00 0 1.000000e+00
## rs150137927 0.000000e+00 0 1.000000e+00
## rs9982091 3.113611e-01 1 5.768461e-01
## rs9982095 1.728000e-01 1 6.776353e-01
## rs9981925 3.113611e-01 1 5.768461e-01
## rs9982014 3.113611e-01 1 5.768461e-01
## rs548330476 0.000000e+00 0 1.000000e+00
## rs9982037 3.113611e-01 1 5.768461e-01
## rs76044243 0.000000e+00 0 1.000000e+00
## rs550678350 0.000000e+00 0 1.000000e+00
## rs138556619 0.000000e+00 0 1.000000e+00
## rs539259554 0.000000e+00 0 1.000000e+00
## rs553128764 0.000000e+00 0 1.000000e+00
## rs541665279 0.000000e+00 0 1.000000e+00
## rs535724590 0.000000e+00 0 1.000000e+00
## rs555426974 0.000000e+00 0 1.000000e+00
## rs575331705 2.336398e-03 1 9.614482e-01
## rs544190631 0.000000e+00 0 1.000000e+00
## rs557893989 0.000000e+00 0 1.000000e+00
## rs577666096 0.000000e+00 0 1.000000e+00
## rs114668878 3.109578e-01 1 5.770931e-01
## rs79244027 3.109578e-01 1 5.770931e-01
## rs528814730 9.433138e-03 1 9.226276e-01
## rs184947944 0.000000e+00 0 1.000000e+00
## rs549014969 0.000000e+00 0 1.000000e+00
## rs368561007 0.000000e+00 0 1.000000e+00
## rs78509661 0.000000e+00 0 1.000000e+00
## rs367901010 0.000000e+00 0 1.000000e+00
## rs533455791 0.000000e+00 0 1.000000e+00
## rs73359717 3.658352e-01 1 5.452832e-01
## rs566637509 0.000000e+00 0 1.000000e+00
## rs374686397 0.000000e+00 0 1.000000e+00
## rs78117715 0.000000e+00 0 1.000000e+00
## rs569131078 0.000000e+00 0 1.000000e+00
## rs147361655 0.000000e+00 0 1.000000e+00
## rs557779439 0.000000e+00 0 1.000000e+00
## rs538688932 0.000000e+00 0 1.000000e+00
## rs556005692 0.000000e+00 0 1.000000e+00
## rs533866185 0.000000e+00 0 1.000000e+00
## rs188585610 0.000000e+00 0 1.000000e+00
## rs575910538 0.000000e+00 0 1.000000e+00
## rs73231921 0.000000e+00 0 1.000000e+00
## rs372655953 0.000000e+00 0 1.000000e+00
## rs575531458 0.000000e+00 0 1.000000e+00
## rs544239350 2.336398e-03 1 9.614482e-01
## rs564489865 0.000000e+00 0 1.000000e+00
## rs571196173 0.000000e+00 0 1.000000e+00
## rs547090650 0.000000e+00 0 1.000000e+00
## rs144815871 5.187727e-01 1 4.713657e-01
## rs369745916 0.000000e+00 0 1.000000e+00
## rs566009716 0.000000e+00 0 1.000000e+00
## rs549150397 0.000000e+00 0 1.000000e+00
## rs148563236 2.142432e-02 1 8.836289e-01
## rs568717616 0.000000e+00 0 1.000000e+00
## rs554701196 0.000000e+00 0 1.000000e+00
## rs551426680 0.000000e+00 0 1.000000e+00
## rs181040610 0.000000e+00 0 1.000000e+00
## rs533521056 0.000000e+00 0 1.000000e+00
## rs553841651 0.000000e+00 0 1.000000e+00
## rs142946068 0.000000e+00 0 1.000000e+00
## rs151083361 0.000000e+00 0 1.000000e+00
## rs139808880 0.000000e+00 0 1.000000e+00
## rs114296721 2.142432e-02 1 8.836289e-01
## rs884166 1.738793e+00 1 1.872918e-01
## rs59484614 2.024095e-01 1 6.527830e-01
## rs564596145 0.000000e+00 0 1.000000e+00
## rs143657279 0.000000e+00 0 1.000000e+00
## rs185165758 0.000000e+00 0 1.000000e+00
## rs555672759 2.336398e-03 1 9.614482e-01
## rs146793180 2.336398e-03 1 9.614482e-01
## rs535215353 0.000000e+00 0 1.000000e+00
## rs562652985 0.000000e+00 0 1.000000e+00
## rs113504625 0.000000e+00 0 1.000000e+00
## rs551672207 0.000000e+00 0 1.000000e+00
## rs188011118 2.336398e-03 1 9.614482e-01
## rs533786198 0.000000e+00 0 1.000000e+00
## rs547075987 0.000000e+00 0 1.000000e+00
## rs117076712 0.000000e+00 0 1.000000e+00
## rs4327307 1.481481e-01 1 7.003114e-01
## rs555567995 0.000000e+00 0 1.000000e+00
## rs569130560 0.000000e+00 0 1.000000e+00
## rs538012773 0.000000e+00 0 1.000000e+00
## rs558434142 0.000000e+00 0 1.000000e+00
## rs182251131 0.000000e+00 0 1.000000e+00
## rs540668517 0.000000e+00 0 1.000000e+00
## rs374623858 0.000000e+00 0 1.000000e+00
## rs145568947 2.336398e-03 1 9.614482e-01
## rs187279637 0.000000e+00 0 1.000000e+00
## rs145134740 0.000000e+00 0 1.000000e+00
## rs1041449 1.075746e-01 1 7.429225e-01
## rs373907767 0.000000e+00 0 1.000000e+00
## rs377062920 0.000000e+00 0 1.000000e+00
## rs149110927 3.844785e-02 1 8.445466e-01
## rs547402085 0.000000e+00 0 1.000000e+00
## rs191773870 0.000000e+00 0 1.000000e+00
## rs373706792 0.000000e+00 0 1.000000e+00
## rs116326720 9.433138e-03 1 9.226276e-01
## rs569370190 0.000000e+00 0 1.000000e+00
## rs538435670 0.000000e+00 0 1.000000e+00
## rs558064990 2.336398e-03 1 9.614482e-01
## rs560247965 0.000000e+00 0 1.000000e+00
## rs572217056 0.000000e+00 0 1.000000e+00
## rs554479236 0.000000e+00 0 1.000000e+00
## rs143146719 3.844785e-02 1 8.445466e-01
## rs543002921 0.000000e+00 0 1.000000e+00
## rs148242543 0.000000e+00 0 1.000000e+00
## rs112245925 0.000000e+00 0 1.000000e+00
## rs181783732 0.000000e+00 0 1.000000e+00
## rs564918415 0.000000e+00 0 1.000000e+00
## rs527609581 0.000000e+00 0 1.000000e+00
## rs71318571 0.000000e+00 0 1.000000e+00
## rs560813074 2.336398e-03 1 9.614482e-01
## rs146950418 0.000000e+00 0 1.000000e+00
## rs137913263 0.000000e+00 0 1.000000e+00
## rs569458711 2.336398e-03 1 9.614482e-01
## rs143380592 2.336398e-03 1 9.614482e-01
## rs73903411 8.065844e-01 1 3.691320e-01
## rs374097887 0.000000e+00 0 1.000000e+00
## rs185923975 0.000000e+00 0 1.000000e+00
## rs77667933 0.000000e+00 0 1.000000e+00
## rs567914620 0.000000e+00 0 1.000000e+00
## rs536694562 0.000000e+00 0 1.000000e+00
## rs556589974 0.000000e+00 0 1.000000e+00
## rs576517652 0.000000e+00 0 1.000000e+00
## rs150507723 0.000000e+00 0 1.000000e+00
## rs558980572 0.000000e+00 0 1.000000e+00
## rs75579202 2.545009e-01 1 6.139235e-01
## rs374915335 0.000000e+00 0 1.000000e+00
## rs541176607 0.000000e+00 0 1.000000e+00
## rs4334565 2.161690e+00 1 1.414890e-01
## rs574387604 0.000000e+00 0 1.000000e+00
## rs572013450 2.336398e-03 1 9.614482e-01
## rs182567491 2.142432e-02 1 8.836289e-01
## rs713399 6.912000e-01 1 4.057565e-01
## rs76728412 1.312079e-02 1 9.088050e-01
## rs190956901 0.000000e+00 0 1.000000e+00
## rs184096964 0.000000e+00 0 1.000000e+00
## rs139423940 0.000000e+00 0 1.000000e+00
## rs193020956 1.211511e-01 1 7.277894e-01
## rs149635596 0.000000e+00 0 1.000000e+00
## rs371526198 0.000000e+00 0 1.000000e+00
## rs713400 0.000000e+00 0 1.000000e+00
## rs185299923 9.433138e-03 1 9.226276e-01
You can save a text file version of this table using the following code:
write.table(HWE,file="TMPRSS2_YRI.HWE.txt")
Now, you’ll have a saved copy of the results of your HWE analysis of all the variants in ACE2 or TMPRSS2 for your population in your workspace!
If you’ve got the text file saved, you can close R Studio and complete the rest of this exercise using the analysis results in the text file. You can look at it in the SCC using either gedit:
gedit TMPRSS2_YRI.HWE.txt &
Or in the terminal space itself using the less command (don’t forget, to leave the less output, press ‘q’):
less TMPRSS2_YRI.HWE.txt
Now, on the right of your table is a list of P-values. Look over your P-values to see if any of them are significant (p < 0.05). If any of your P-values are significant, this means that that SNP is not in Hardy-Weinberg Equilibrium!
Here’s some easy R code to filter out JUST those SNPs that are significant:
library(tidyverse)
HWE.sig <-
HWE %>%
as_tibble(HWE) %>%
mutate(snp = rownames(HWE)) %>%
rename(chi2 = 'chi^2', pval = 'Pr(chi^2 >)') %>%
select(snp,chi2,df,pval) %>%
filter(pval <= 0.05)
HWE.sig
## # A tibble: 13 x 4
## snp chi2 df pval
## <chr> <dbl> <dbl> <dbl>
## 1 rs7278739 4.30 1 3.82e- 2
## 2 rs2070790 4.05 1 4.42e- 2
## 3 rs2070787 4.05 1 4.42e- 2
## 4 rs2070786 4.23 1 3.98e- 2
## 5 rs28524972 4.54 1 3.32e- 2
## 6 rs144192191 5.17 1 2.30e- 2
## 7 rs141400887 4.63 1 3.14e- 2
## 8 rs9983252 4.12 1 4.23e- 2
## 9 rs10154090 4.32 1 3.77e- 2
## 10 rs9305745 5.06 1 2.45e- 2
## 11 rs114027875 47.3 1 6.01e-12
## 12 rs7282202 8.18 1 4.25e- 3
## 13 rs57197209 11.4 1 7.41e- 4
In this case the YRI population has 13 SNPs out of HWE in TMPRSS2. Record the SNP ID numbers that have significant P-values for your population, because they will be important for what we do next.
Hopefully, some of you have found that there are SNPs in your population that are not in Hardy-Weinberg Equilibrium. That’s exciting! But, think back to the Wigginton paper. There is a large Type 1 Error rate for Hardy-Weinberg tests performed using a Chi-Squared test. To account for this, we will test our significant SNPs for Type 1 Error using the Shiny App calculator Becca De Camp created for us. Here’s how we will do that:
Here, you will be interacting with a user interface on this webpage that was created to calculate “true” Hardy-Weinberg Equilibrium P-values for one SNP at a time. The input for this UI will be the genotype frequencies from your population (i.e. the number of people from the population with the heterozygous genotype, the number of people with one homozygous genotype, etc.). You will find this information on Ensembl by doing the following:
Make a shortlist of FOUR SNP ID numbers of the SNPs that you found not to be in Hardy-Weinberg Equilibrium. When making this list, check the degrees of freedom, column. This UI will only accept SNPs with only two alleles, or one degree of freedom. Therefore, even though there are SNPs with more than two alleles that are not in Hardy-Weinberg Equilibrium, please ignore them for now. If you do not have four SNPs that meet these criteria, fill in your list with SNPs at random. Ultimately, EVERYONE should be looking at FOUR SNPs each.
Go to the ACE2 or TMPRSS2 gene page on Ensembl and click on the Variant table tab so you get this familiar page:
In the search bar at the top right of the variant table, type in the SNP ID number for one of your SNPs that’s not in Hardy-Weinberg Equilibrium and search it. I will use the ACE2 SNP rs2074192 to demonstrate this.
When you have found your SNP in the table, you should first make note of the consequence type of the SNP. For example, rs2074192 is an intronic variant. How do you think the consequence type will affect the phenotypic outcome of this SNP?
Next, click on the SNP ID to get to the SNP page. From there, click on the Population Genetics tab to get the population genetics table.
Find the Genotype Frequency column in this table:
Find your population in the table. Record the genotype frequency count numbers (heterozygous, homozygous 1 and homozygous 2). These will be your input to the UI below. Using this method, find these genotype frequencies for each SNP you want to test.
NOTE: before you’ve entered any data into the UI, it will show a result error (“Error: missing value where TRUE/FALSE needed”). This just means there’s nothing entered yet. Once you’ve entered your genotype counts, that error message will be replaced by the calculated “True” Hardy-Weinberg Equilibrium P-value for your SNP.
Record your results, and make sure to record the SNP ID of each SNP you are testing. The last thing we’ll do is think about what the calculations you did today mean.
Think about the results you produced today in the context of your population. Here are some guiding questions to help you:
If your population had SNPs that were not in Hardy-Weinberg Equilibrium, what kind of consequence did the SNP(s) have? How would that consequence type affect the phenotype?
Does a SNP not being in Hardy-Weinberg Equilibrium mean that evolutionary forces are acting on it? How might evolutionary forces be acting on ACE2 if a SNP is not in Hardy-Weinberg Equilibrium?
You can actually run the same analysis we just did using a single line of code in vcftools, thanks to the hardy option.
If you’d like to try, load the vcftools module in the SCC:
module load vcftools
Now, enter the following code. This is essentially telling vcftools to run the “true” HWE test on each variant site, and to give the output as a file called “TMPRSS2_YRI.hwe”:
vcftools --gzvcf TMPRSS2_YRI.vcf.gz --hardy --out TMPRSS2_YRI
Once it’s run, you can read the output file using the less command:
less TMPRSS2_YRI.hwe
Notice that the output is giving us both more and less information than our analysis in R Studio. Rather than the variant name, we’re given just the variant position (‘POS’) on which chromosome (‘CHR’). We’re also, however, given the observed (‘OBS(HOM1/HET/HOM2)’) and expected (‘E(HOM1/HET/HOM2)’) counts of individuals under “true” HWE with each genotype, along with the “true” HWE Chi-Squared (‘ChiSq_HWE’) and P-values (‘P_HWE’). It’s also giving us P-values regarding whether there is a significant deficit (‘P_HET_DEFICIT’) or excess (‘P_HET_EXCESS’) in heterozygotes, which can tell us a bit more about why the population is out of Hardy-Weinberg equilibrium (we’ll learn more about this in class).
To see only the significant SNPs for this output, we just need to import the file into R:
vcftools.HWE<-read.table("TMPRSS2_YRI.hwe",header=TRUE)
head(vcftools.HWE)
## CHR POS OBS.HOM1.HET.HOM2. E.HOM1.HET.HOM2. ChiSq_HWE P_HWE
## 1 21 41464332 108/0/0 108.00/0.00/0.00 NaN 1
## 2 21 41464372 107/1/0 107.00/1.00/0.00 0.002336398 1
## 3 21 41464375 108/0/0 108.00/0.00/0.00 NaN 1
## 4 21 41464414 106/2/0 106.01/1.98/0.01 0.009433138 1
## 5 21 41464462 108/0/0 108.00/0.00/0.00 NaN 1
## 6 21 41464497 108/0/0 108.00/0.00/0.00 NaN 1
## P_HET_DEFICIT P_HET_EXCESS
## 1 1 1.0000000
## 2 1 1.0000000
## 3 1 1.0000000
## 4 1 0.9953488
## 5 1 1.0000000
## 6 1 1.0000000
And then filter by P_HWE
, which represents the HWE p-value:
library(tidyverse)
vcftools.HWE.sig <-
vcftools.HWE %>%
as_tibble(vcftools.HWE) %>%
select(CHR,POS,ChiSq_HWE,P_HWE) %>%
filter(P_HWE <= 0.05)
vcftools.HWE.sig
## # A tibble: 10 x 4
## CHR POS ChiSq_HWE P_HWE
## <int> <int> <dbl> <dbl>
## 1 21 41469311 4.30 0.0429
## 2 21 41470276 4.23 0.0476
## 3 21 41470844 4.54 0.0417
## 4 21 41473914 5.17 0.0397
## 5 21 41475153 4.63 0.0395
## 6 21 41491796 4.32 0.0469
## 7 21 41492147 5.06 0.0285
## 8 21 41500731 47.3 0.0140
## 9 21 41518025 8.18 0.00526
## 10 21 41518039 11.4 0.00107