Skip to content

Instantly share code, notes, and snippets.

@msuiche
Created January 30, 2017 15:29
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save msuiche/b4b35c7cbb5baa1e4ba7a3d115c2987d to your computer and use it in GitHub Desktop.
Save msuiche/b4b35c7cbb5baa1e4ba7a3d115c2987d to your computer and use it in GitHub Desktop.
Powershell script to parse your SNP
# http://slatestarcodex.com/2014/11/12/how-to-use-23andme-irresponsibly/
$indexSNPName = 0
$indexSNPValue = 3
$fileName = "genome_Peter_Parker_Full_20170130051759.txt"
Write-Host " --- "
$rs909525 = gci $fileName | Select-String 'rs909525\t'
$rs909525 = $rs909525.Line.Split("`t")
Write-Host "[+]"$rs909525[$indexSNPName] "(Warrior Gene): ("$rs909525[$indexSNPValue]")"
switch ($rs909525[$indexSNPValue]) {
"T" { " Not Present. Probably MAOA 4 or 5 repeats: not Warrior Gene."}
"C" { " Probably one Warrior Gene and one non-Warrior Gene."}
"CC" { " Perhaps MAOA 3 repeats: Warrior Gene?"}
default { " [?] Can't identify: '" + $rs909525[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs53576 = gci $fileName | Select-String 'rs53576\t'
$rs53576 = $rs53576.Line.Split("`t")
Write-Host "[+]"$rs53576[$indexSNPName] "(Cuddle/Trust Hormone - emotional warmth - oxytocin receptor) : ("$rs53576[$indexSNPValue]")"
switch ($rs53576[$indexSNPValue]) {
"AA" { " Lack of empathy? empathy"}
"AG" { " Lack of empathy? empathy"}
"GG" { " More empathy and sociability. Optimistic and empathetic; handle stress well"}
default { " [?] Can't identify: '" + $rs53576[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs4680 = gci $fileName | Select-String 'rs4680\t'
$rs4680 = $rs4680.Line.Split("`t")
Write-Host "[+]"$rs4680[$indexSNPName] "(Worrier gene): ("$rs4680[$indexSNPValue]")"
switch ($rs4680[$indexSNPValue]) {
"AA" { " Smart and Anxious. advantage in memory and attention tasks"}
"AG" { " Inbetween (Smart Anxious) & (Dumb and Calm)"}
"GG" { " Dumb and Calm. (G;G) carriers deprived of sleep respond quite well to 2x 100mg modafinil in terms of improved vigor and well-being"}
default { " [?] Can't identify: '" + $rs4680[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs7632287 = gci $fileName | Select-String 'rs7632287\t'
$rs7632287 = $rs7632287.Line.Split("`t")
Write-Host "[+]"$rs7632287[$indexSNPName] "(Divorce gene (female) oxytocin receptor): ("$rs7632287[$indexSNPValue]")"
switch ($rs7632287[$indexSNPValue]) {
"AA" { " Highest risk"}
"AG" { " Mid risk"}
"GG" { " Lower risk. "}
default { " [?] Can't identify: '" + $rs7632287[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs11174811 = gci $fileName | Select-String 'rs11174811\t'
$rs11174811 = $rs11174811.Line.Split("`t")
Write-Host "[+]"$rs11174811[$indexSNPName] "(spousal satisfaction (male) oxytocin receptor): ("$rs11174811[$indexSNPValue]")"
switch ($rs11174811[$indexSNPValue]) {
"AA" { " Highest spousal satisfaction "}
"AC" { " Mid spousal satisfaction "}
"CC" { " Lowest spousal satisfaction "}
default { " [?] Can't identify: '" + $rs11174811[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs25531 = gci $fileName | Select-String 'rs25531\t'
$rs25531 = $rs25531.Line.Split("`t")
Write-Host "[+]"$rs25531[$indexSNPName] "(Orchid gene (kids) serotonin transporter.): ("$rs25531[$indexSNPValue]")"
switch ($rs25531[$indexSNPValue]) {
"TT" { " Orchid (high potential). long form of 5-HTTLPR. less sensitive to pain"}
"CC" { " Normal. short form of 5-HTTLPR. lower levels of serotonin, slightly less happy, benefits from more support "}
default { " [?] Can't identify: '" + $rs25531[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs1800955 = gci $fileName | Select-String 'rs1800955\t'
$rs1800955 = $rs1800955.Line.Split("`t")
Write-Host "[+]"$rs1800955[$indexSNPName] "(Adventure gene - dopamine receptor): ("$rs1800955[$indexSNPValue]")"
switch ($rs1800955[$indexSNPValue]) {
"CC" { " Novelty seeker. increased susceptibility to novelty seeking"}
"CT" { " Novelty seeker. increased susceptibility to novelty seeking"}
"TT" { " Normal"}
default { " [?] Can't identify: '" + $rs1800955[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs2760118 = gci $fileName | Select-String 'rs2760118\t'
$rs2760118 = $rs2760118.Line.Split("`t")
Write-Host "[+]"$rs2760118[$indexSNPName] "(Smart gene - succinate semialdehyde dehydrogenase): ("$rs2760118[$indexSNPValue]")"
switch ($rs2760118[$indexSNPValue]) {
"CC" { " Smarter and can expect to live longer (strong)"}
"CT" { " Smarter and can expect to live longer (mid)"}
"TT" { " Normal"}
default { " [?] Can't identify: '" + $rs2760118[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs6311 = gci $fileName | Select-String 'rs6311\t'
$rs6311 = $rs6311.Line.Split("`t")
Write-Host "[+]"$rs6311[$indexSNPName] "(Social Success of Young Males - Serotonin receptor): ("$rs6311[$indexSNPValue]")"
switch ($rs6311[$indexSNPValue]) {
"CC" { " Popular. 3.6x increased risk of sexual dysfunction when taking SSRI Antidepressants"}
"CT" { " Popular. Normal risk of sexual dysfunction when taking SSRI Antidepressants."}
"TT" { " Normal (lower) risk of sexual dysfunction when taking SSRI Antidepressants."}
default { " [?] Can't identify: '" + $rs6311[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs6265 = gci $fileName | Select-String 'rs6265\t'
$rs6265 = $rs6265.Line.Split("`t")
Write-Host "[+]"$rs6265[$indexSNPName] "(depression-linked BDNF system - resistence to depression): ("$rs6265[$indexSNPValue]")"
switch ($rs6265[$indexSNPValue]) {
"CC" { " introversion; depression resistant; impaired motor skills learning"}
"CT" { " Slightly increased risk for ADHD or depression"}
"TT" { " Normal"}
default { " [?] Can't identify: '" + $rs6265[$indexSNPValue] + "'" }
}
Write-Host " --- "
$rs41310927 = gci $fileName | Select-String 'rs41310927\t'
$rs41310927 = $rs41310927.Line.Split("`t")
Write-Host "[+]"$rs41310927[$indexSNPName] "(Tonal language): ("$rs41310927[$indexSNPValue]")"
switch ($rs41310927[$indexSNPValue]) {
"CC" { " Great at tonal languages (Asians)"}
"CT" { " inbetween ability to differentiate tone" }
"TT" { " lowest ability to differentiate tones"}
default { " [?] Can't identify: '" + $rs41310927[$indexSNPValue] + "'" }
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment