Skip to content

Python class which reads and interprets 23andMe raw genotype data.

License

Notifications You must be signed in to change notification settings

naturale0/GenoMe

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

7 Commits
 
 
 
 
 
 
 
 

Repository files navigation

GenoMe: interpret your 23andMe data

GenoMe is a Python module which reads and interprets 23andMe raw genotype data. It Support simple wellness report from 23andMe, and searching your genotype with related traits at SNPedia.

Quick-look

>>> from GenoMe import Geno
>>> g = Geno("/path/to/your/raw/data.txt")
  • Geno.report_wellness() prints out wellness report provided by 23andMe API.
>>> g.report_wellness()
=============================== WELLNESS REPORT =============================== 

< Alcohol Flush Reaction >
--------------------------------------------------------------------------------
The marker we tested comes in two different forms, the G variant and the A varia
nt. The A variant results in an enzyme that is less efficient at breaking down a
cetaldehyde. The A variant is also known as c.1510G>A, Glu487Lys, and Glu504Lys.
 This marker has been studied the most in people of East Asian descent.

  rsID: rs671
  Yours: G;G - homozygous
  Detail: 2 | Alcohol Flush: Normal, doesn't flush. Normal hangovers. Normal ris
k of Alcoholism. Normal risk of Esophageal Cancer. Disulfiram is effective for a
lcoholism. | 60.2% of JPT

...
  • Geno.type_id(rsid) prints out specific SNP of input rsid and its related trait written in SNPedia.
>>> g.type_id("rs28897696")
rs28897696 | G;G | 0 | normal | 0.0% of JPT
rsid chr pos genotype
492260 rs28897696 17 38469446 GG
  • Geno.type_trait(trait) prints out all SNPs related to input trait, written in SNPedia.
>>> g.type_trait("baldness")
rs6152 | G;G | 0.5 | able to go bald | 100.0% of JPT
rs2223841 | T;T | 1.2 | more likely to go bald before age 40 | 100.0% of JPT
rs6152 | G;G | 0.5 | able to go bald | 100.0% of JPT
rs2223841 | T;T | 1.2 | more likely to go bald before age 40 | 100.0% of JPT
rs2180439 | C;T | 2.5 | Increased risk of Male Pattern Baldness. | 44.2% of JPT
rs11683401 | C;T |  | --% of JPT
rs2073963 | G;G | 2.5 | increased risk of baldness | 18.6% of JPT
rs1511061 | T;T |  | --% of JPT
 * more details at https://www.snpedia.com/index.php/baldness

Reference Population

You can also change the population used to report genotype ratio, other than 'JPT'. report_wellness, type_id, and type_trait supports ref_population option to input population code.

>>> g.report_wellness(ref_population="CEU")
>>> g.type_id("rs28897696", ref_population="CHB")
>>> g.type_trait("baldness", ref_population="ASW")

For more information about populations and their code, see this page.

About

Python class which reads and interprets 23andMe raw genotype data.

Topics

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages