text: user1574: I'm trying to play around with my family's raw genetic data from Ancestry.com so that I can do some genetic analyses, and I noticed in the AncestryDNA raw data files that some regions of the genome have more than one variant associated with them and I would like to know why this is. The end goal is get this into a regularly formatted VCF, but I feel that these double variants will probably mess things up. So for example one might see something similar to this in the tab-delimited file:
  rs6010164      22    50515270    T    T
  rs267607236    22    50515270    I    I
  
  (Header line: rsid    chromosome    position    allele1    allele2)
  I understand that when rsid numbers were created, sometimes multiple rsid numbers were accidentally assigned to the same region of the genome, but  I do not understand how two separate variants can be assigned to one single region in the genome. I could be wrong but as I understand it, I should only see one line here (either a "T    T" or an "I    I") because genotypes are reported in comparison to the human reference genome and so each position should uniquely refer to a single nucleotide in the said human reference genome.
  But in this example case one line says there is a "T    T" genotype at this position, and the other line says there is an insertion(see "How do I read my raw DNA data?" here).  One possible explanation I thought of is if this would just be saying that a T was inserted, but if this is the case, then as I understand it the reference allele present in the individual should still reported, just in conjunction with the allele at the index. (So assuming an extra "A" was inserted after the T, which could be the reference or alternative allele, something like this should be recorded: rs6010164  22  50515270    TA  TA).
  Can anyone help me understand why there are two separate variants associated with chromosomal position? Thanks!
  
  user581: It's pretty common for SNPs and indels that map to the same position in the genome to be assigned different rsIDs, particularly when the indel spans multiple base pairs in the reference. You can get more info about these particular variants here and here. Or you can see how both of these variants relate to the reference sequence here.
  Without info about the reference allele it's pretty confusing (particularly because they appear to be reporting the opposite strand from the database (G vs. C/T), but I would interpret this to mean that both alleles match the reference for rs267607236 (they have the insertion), but do not match the reference for rs6010164 (G-C pair rather than T-A). This is also the most common genotype at this position, being found in almost 80% of people in most populations. 
  
  user298: First of all, take these data with a very large pinch of salt. This sort of targeted analysis is not designed to produce high quality genetic data but, to give an idea of a sample's ancestry. Given a file like the one you describe with no quality information, you have no way of knowing how reliable any of the variants called actually is. From personal experience with another, similar company, I expect that they aren't reliable at all. 
  I have tested four individuals using a commercial tool like the one you describe and while I don't doubt the validity of their main results (those displayed on their website), their raw data showed combinations of variants that were extremely unlikely to be real given that the individuals were still breathing.
  This isn't to say that such companies don't do their job well, only that their job is limited to what they actually display. When you look into their raw data, you're on your own. Presumably, they apply various filters before displaying the main results, but these are most likely not applied to the raw data. So if you want to do your own analyses, this is unlikely to be a reliable dataset. Also, I would strongly urge you to seek the help of an expert when interpreting any data you do have. 
  
  That said, the rsIDs are simply identifiers of a specific genetic variation, not of a location. Any genomic position can have many (essentially infinite) possible variations. For example, a C residue at position 100 of chromosome 1 could become a G, a T or an A. Each of these will have its own, separate rsID. Now, if you also consider insertions and deletions, you can see why you can have essentially limitless possible variations at that specific position. So no, having more than one rsID describing the same position is not strange at all. 
  The way these tests work is they look at a list of possible variants, each of which has its own rsID, and they report what they find for that rsID in the sample. Then, they will report what they find at the rsID's position for each allele. This is explained in the FAQ article you linked to:
  
  The possible observations are A for adenine, C for cytosine, G for guanine, T for thymine, I for insertion’ and ‘D for deletion, or 0 for missing data. Column one provides the identifier (including the #rsID where possible). Columns two and three contain the chromosome and basepair position of the variant using human reference 37.1 coordinates. Columns four and five contain the two alleles observed at this variant (genotype)
  
  Let's have a look at the two rsIds you mention.
  
  rs6010164: This is a single nucleotide variant (SNV) which describes a C residue where the reference genome has a T. 
  rs267607236: this is a deletion at position 50076841, so a very different variant to the above. In this one, the sequence TGAG is deleted. 
  
  Your results are:
  rs6010164      22    50515270    T    T
  rs267607236    22    50515270    I    I
  
  So, your sample doesn't have rs6010164 since it has the same sequence at that position as the reference. However, it does have the insertion as indicated by the two Is.
  The Is and Ds of the insertions and deletions (indels) are going to be the easy bit here. The only way you can parse this sort of file and convert it to VCF is to get the list of target SNV rsIDs, check what they describe as the reference and as the variant and then, for every one, check which of the two is shown in your file.
  If you're not an expert in this field, you will probably be better off hiring someone to do this for you. There are actually companies that offer this sort of service and can take these tab-delimited files as input. I even work for one! 
  
original_text: 1574: <p>I'm trying to play around with my family's raw genetic data from Ancestry.com so that I can do some genetic analyses, and I noticed in the AncestryDNA raw data files that some regions of the genome have more than one variant associated with them and I would like to know why this is. The end goal is get this into a regularly formatted VCF, but I feel that these double variants will probably mess things up. So for example one might see something similar to this in the tab-delimited file:</p>
  
  <pre><code>rs6010164      22    50515270    T    T
  rs267607236    22    50515270    I    I
  </code></pre>
  
  <p>(<em>Header line</em>: rsid    chromosome    position    allele1    allele2)</p>
  
  <p>I understand that when rsid numbers were created, sometimes multiple rsid numbers were accidentally assigned to the same region of the genome, but  I do not understand how two separate variants can be assigned to one single region in the genome. I could be wrong but as I understand it, I should only see one line here (either a "T    T" or an "I    I") because genotypes are reported in comparison to the human reference genome and so each position should uniquely refer to a single nucleotide in the said human reference genome.</p>
  
  <p>But in this example case one line says there is a "T    T" genotype at this position, and the other line says there is an insertion(see "How do I read my raw DNA data?" <a href="https://www.ancestry.com/dna/en/legal/us/faq#raw-6" rel="nofollow noreferrer">here</a>).  One possible explanation I thought of is if this would just be saying that a T was inserted, but if this is the case, then as I understand it the reference allele present in the individual should still reported, just in conjunction with the allele at the index. (So assuming an extra "A" was inserted after the T, which could be the reference or alternative allele, something like this should be recorded: rs6010164  22  50515270    TA  TA).</p>
  
  <p>Can anyone help me understand why there are two separate variants associated with chromosomal position? Thanks!</p>
   581: <p>It's pretty common for SNPs and indels that map to the same position in the genome to be assigned different rsIDs, particularly when the indel spans multiple base pairs in the reference. You can get more info about these particular variants <a href="https://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=6010164" rel="nofollow noreferrer">here</a> and <a href="https://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=267607236" rel="nofollow noreferrer">here</a>. Or you can see how both of these variants relate to the reference sequence <a href="https://www.ncbi.nlm.nih.gov/variation/view/?q=rs267607236&amp;filters=source%3Adbsnp&amp;assm=GCF_000001405.25" rel="nofollow noreferrer">here</a>.</p>
  
  <p>Without info about the reference allele it's pretty confusing (particularly because they appear to be reporting the opposite strand from the database (G vs. C/T), but I would interpret this to mean that both alleles match the reference for rs267607236 (they have the insertion), but do not match the reference for rs6010164 (G-C pair rather than T-A). This is also the most common genotype at this position, being found in almost 80% of people in most populations. </p>
   298: <p>First of all, <strong>take these data with a very large pinch of salt</strong>. This sort of targeted analysis is not designed to produce high quality genetic data but, to give an idea of a sample's ancestry. Given a file like the one you describe with no quality information, you have no way of knowing how reliable any of the variants called actually is. From personal experience with another, similar company, I expect that they aren't reliable at all. </p>
  
  <p>I have tested four individuals using a commercial tool like the one you describe and while I don't doubt the validity of their main results (those displayed on their website), their raw data showed combinations of variants that were extremely unlikely to be real given that the individuals were still breathing.</p>
  
  <p>This isn't to say that such companies don't do their job well, only that their job is limited to what they actually display. When you look into their raw data, you're on your own. Presumably, they apply various filters before displaying the main results, but these are most likely not applied to the raw data. So if you want to do your own analyses, this is unlikely to be a reliable dataset. Also, I would strongly urge you to seek the help of an expert when interpreting any data you do have. </p>
  
  <hr>
  
  <p>That said, the rsIDs are simply identifiers of a specific genetic variation, not of a location. Any genomic position can have many (essentially infinite) possible variations. For example, a C residue at position 100 of chromosome 1 could become a G, a T or an A. Each of these will have its own, separate rsID. Now, if you also consider insertions and deletions, you can see why you can have essentially limitless possible variations at that specific position. So no, having more than one rsID describing the same position is not strange at all. </p>
  
  <p>The way these tests work is they look at a list of possible variants, each of which has its own rsID, and they report what they find for that rsID in the sample. Then, they will report what they find at the rsID's position for each allele. This is explained in <a href="https://www.ancestry.com/dna/en/legal/us/faq#raw-6" rel="nofollow noreferrer">the FAQ article</a> you linked to:</p>
  
  <blockquote>
    <p>The possible observations are A for adenine, C for cytosine, G for guanine, T for thymine, I for insertion’ and ‘D for deletion, or 0 for missing data. Column one provides the identifier (including the #rsID where possible). Columns two and three contain the chromosome and basepair position of the variant using human reference 37.1 coordinates. Columns four and five contain the two alleles observed at this variant (genotype)</p>
  </blockquote>
  
  <p>Let's have a look at the two rsIds you mention.</p>
  
  <ol>
  <li><p><a href="https://varsome.com/variant/hg38/rs6010164" rel="nofollow noreferrer">rs6010164</a>: This is a single nucleotide variant (SNV) which describes a C residue where the reference genome has a T. </p></li>
  <li><p><a href="https://varsome.com/variant/hg38/rs267607236" rel="nofollow noreferrer">rs267607236</a>: this is a deletion at position 50076841, so a very different variant to the above. In this one, the sequence TGAG is deleted. </p></li>
  </ol>
  
  <p>Your results are:</p>
  
  <pre><code>rs6010164      22    50515270    T    T
  rs267607236    22    50515270    I    I
  </code></pre>
  
  <p>So, your sample doesn't have rs6010164 since it has the same sequence at that position as the reference. However, it does have the insertion as indicated by the two Is.</p>
  
  <p>The Is and Ds of the insertions and deletions (indels) are going to be the easy bit here. The only way you can parse this sort of file and convert it to VCF is to get the list of target SNV rsIDs, check what they describe as the reference and as the variant and then, for every one, check which of the two is shown in your file.</p>
  
  <p>If you're not an expert in this field, you will probably be better off hiring someone to do this for you. There are actually companies that offer this sort of service and can take these tab-delimited files as input. I even work for one! </p>
  