Answers.Body: <p>A csv file with RSID, chromosome, position and result is enough for what you want to do and these are the core columns of a VCF (which is just a TSV with some headers describing how it was made).</p>
  <p>Given you have ~700k rows I suspect that yes there will be genetic health-related SNVs (single nucleotide variant). Disclaimer; these data could include info pertaining to your health and your family's health, I strongly advise that you speak to a genetic counsellor to understand these things.</p>
  <p>Gnomad is a population database so generally speaking (plenty of exceptions of course) if there is a SNV in there more than a few times its probably not causing disease. <a href="https://www.ncbi.nlm.nih.gov/clinvar/" rel="nofollow noreferrer">https://www.ncbi.nlm.nih.gov/clinvar/</a> and <a href="https://omim.org/" rel="nofollow noreferrer">https://omim.org/</a> are example of disease databases.</p>
  <p>Here's a pandas example with clinvar (<a href="https://ftp.ncbi.nlm.nih.gov/pub/clinvar/vcf_GRCh38/" rel="nofollow noreferrer">https://ftp.ncbi.nlm.nih.gov/pub/clinvar/vcf_GRCh38/</a> you'll need to match your version).</p>
  <pre><code>import pandas as pd
  df = pd.read_csv(
      'clinvar.vcf.gz',
      sep='\t',
      comment = '#',
      header=None, 
      names = ['CHROM','POS','ID','REF','ALT','QUAL','FILTER','INFO']
  )
  df.head()
  
  CHROM   POS ID  REF ALT QUAL    FILTER  INFO
  0   1   930188  846933  G   A   .   .   ALLELEID=824438;CLNDISDB=MedGen:CN517202;CLNDN...
  1   1   930203  972363  C   T   .   .   ALLELEID=959431;CLNDISDB=MedGen:CN517202;CLNDN...
  2   1   930248  789256  G   A   .   .   AF_ESP=0.00347;AF_EXAC=0.00622;AF_TGP=0.00280;...
  3   1   930275  969662  T   G   .   .   ALLELEID=959432;CLNDISDB=MedGen:CN517202;CLNDN...
  4   1   930336  843786  G   A   .   .   ALLELEID=824439;CLNDISDB=MedGen:CN517202;CLNDN...
  </code></pre>
  <p>Then lets say one of your 700k SNVs is at chrom 1, position 930336 and you have an A</p>
  <pre><code>print(*df.query('(CHROM == 1) &amp; (POS == 930336) &amp; (ALT == &quot;A&quot;)')['INFO'].str.split(';'))
  </code></pre>
  <p>Gives a list from which you could pull out CLNSIG which here is Uncertain_significance</p>
  <pre><code>['ALLELEID=824439', 'CLNDISDB=MedGen:CN517202', 'CLNDN=not_provided', 'CLNHGVS=NC_000001.11:g.930336G&gt;A', 'CLNREVSTAT=criteria_provided,_single_submitter', 'CLNSIG=Uncertain_significance', 'CLNVC=single_nucleotide_variant', 'CLNVCSO=SO:0001483', 'GENEINFO=SAMD11:148398', 'MC=SO:0001583|missense_variant', 'ORIGIN=1']
  </code></pre>
  <p>So you could easily parse your 700k SNVs and get their respective clinical info this way. Please think about this carefully before you do it though, there are many implications including life insurance and mental health! I'm of the opinion that this sort of data is useful but MANY people disagree. Further, this advice is general in nature only, I accept no liabilty for what you or anyone else chooses to do with this publicly available data.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2021-02-11T21:57:23.887Z
Answers.LastActivityDate: 2021-02-11T21:57:23.887Z
Answers.UserDisplayName: Liam McIntyre
Body: <p>I had my DNA tested by Myheritage and they sent me a csv file with RSID, Chromosome, position and result (which base) with about 700,000 rows.
  I understand most analyses of DNA use VCF files but is there anything i can do with this csv file i.e. check for genetic health-related genes?
  I am not a bioinformaticist but I am a scientist and I can use R and python.  Ive heard of the Gnomad database but not sure if i can match things to my csv file.</p>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2021-02-11T15:19:03.910Z
Comments.Text: I think people often use promethease for this kind of thing.
Comments.UserDisplayName: user438383
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2021-02-13T16:06:40.963Z
Comments.Text: Check out data.table and merge with gnomAD or better, CLINVAR. That might give you some pointers.
Comments.UserDisplayName: Ram RS
Tags: <vcf><phylogenetics>
Title: How can I use my Myheritage DNA results file for further analysis?
ThreadText: # How can I use my Myheritage DNA results file for further analysis?
  
  (Asked by: username_0 on 2021-02-11)
  
  I had my DNA tested by Myheritage and they sent me a csv file with RSID, Chromosome, position and result (which base) with about 700,000 rows.
  I understand most analyses of DNA use VCF files but is there anything i can do with this csv file i.e. check for genetic health-related genes?
  I am not a bioinformaticist but I am a scientist and I can use R and python. Ive heard of the Gnomad database but not sure if i can match things to my csv file.
  
  **Tags:** vcf, phylogenetics
  
  > **Comments:**
  > **username_1**: I think people often use promethease for this kind of thing.  
  > **username_2**: Check out data.table and merge with gnomAD or better, CLINVAR. That might give you some pointers.  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_3 on 2021-02-11. Score: 1
  
  A csv file with RSID, chromosome, position and result is enough for what you want to do and these are the core columns of a VCF (which is just a TSV with some headers describing how it was made).
  
  Given you have ~700k rows I suspect that yes there will be genetic health-related SNVs (single nucleotide variant). Disclaimer; these data could include info pertaining to your health and your family's health, I strongly advise that you speak to a genetic counsellor to understand these things.
  
  Gnomad is a population database so generally speaking (plenty of exceptions of course) if there is a SNV in there more than a few times its probably not causing disease. <https://www.ncbi.nlm.nih.gov/clinvar/> and <https://omim.org/> are example of disease databases.
  
  Here's a pandas example with clinvar (<https://ftp.ncbi.nlm.nih.gov/pub/clinvar/vcf_GRCh38/> you'll need to match your version).
  
  ```
  import pandas as pd
  df = pd.read_csv(
      'clinvar.vcf.gz',
      sep='\t',
      comment = '#',
      header=None, 
      names = ['CHROM','POS','ID','REF','ALT','QUAL','FILTER','INFO']
  )
  df.head()
  
  CHROM   POS ID  REF ALT QUAL    FILTER  INFO
  0   1   930188  846933  G   A   .   .   ALLELEID=824438;CLNDISDB=MedGen:CN517202;CLNDN...
  1   1   930203  972363  C   T   .   .   ALLELEID=959431;CLNDISDB=MedGen:CN517202;CLNDN...
  2   1   930248  789256  G   A   .   .   AF_ESP=0.00347;AF_EXAC=0.00622;AF_TGP=0.00280;...
  3   1   930275  969662  T   G   .   .   ALLELEID=959432;CLNDISDB=MedGen:CN517202;CLNDN...
  4   1   930336  843786  G   A   .   .   ALLELEID=824439;CLNDISDB=MedGen:CN517202;CLNDN...
  
  ```
  
  Then lets say one of your 700k SNVs is at chrom 1, position 930336 and you have an A
  
  ```
  print(*df.query('(CHROM == 1) & (POS == 930336) & (ALT == "A")')['INFO'].str.split(';'))
  
  ```
  
  Gives a list from which you could pull out CLNSIG which here is Uncertain\_significance
  
  ```
  ['ALLELEID=824439', 'CLNDISDB=MedGen:CN517202', 'CLNDN=not_provided', 'CLNHGVS=NC_000001.11:g.930336G>A', 'CLNREVSTAT=criteria_provided,_single_submitter', 'CLNSIG=Uncertain_significance', 'CLNVC=single_nucleotide_variant', 'CLNVCSO=SO:0001483', 'GENEINFO=SAMD11:148398', 'MC=SO:0001583|missense_variant', 'ORIGIN=1']
  
  ```
  
  So you could easily parse your 700k SNVs and get their respective clinical info this way. Please think about this carefully before you do it though, there are many implications including life insurance and mental health! I'm of the opinion that this sort of data is useful but MANY people disagree. Further, this advice is general in nature only, I accept no liabilty for what you or anyone else chooses to do with this publicly available data.
  
  
  ---
  
  