text: user2519: Let's say I'm sequencing a donor, and I have 20X coverage. And let's say all the reads covering a particular base has the same allele at the position. Then, is that position considered a "variant" with only 1 allele in that position or is that not even considered a variant?
  The reason I'm asking is, I was given some read data, and it's in the form location:variant - location:variant.
  But there was often the case where for all the variants at a particular location, there was only one allele. But this data is labelled to only have variants. Hence, how is a particular location deemed a "variant."
  For example, lets say I'm told that the following is "variant reads given in the form read# location:variant - location:variant"
  read 1 6523432: A - 6710432 T
  read 2 6523432: C - 6710432 G
  read 3 6523432: A - 6710432 T
  read 4 6523432: C - 6710432 G
  
  Then, why are these two locations considered variants? It seems like it's just a heterozygous locations. Why were these locations considered to carry variants? Similarly, 
  read 1 6523432: A - 6710432 T
  read 2 6523432: A - 6710432 T
  read 3 6523432: A - 6710432 T
  read 4 6523432: A - 6710432 T
  
  In those reads, it just seems like it's homozygous. Why is it being called a variant? How was it classified as such?
  Put another way, if I'm given 100000 reads, and I said I only want to look at informative reads, and so I'm given a subset with only 100 reads, and those 100 reads look something like the above, then how were these reads selected?
  
  user298: 
  Set 1
  read 1 6523432: A - 6710432 T
  read 2 6523432: C - 6710432 G
  read 3 6523432: A - 6710432 T
  read 4 6523432: C - 6710432 G
  
  I understand this as showing the location of 2 variants, supported by 2 reads each. Now, if the reference has an A at position 6523432, your sample could be heterozygous for an A->C variant at position 6523432. If the reference has a C at that position, then your sample could be heterozygous for C->A. If the reference has something else again, for example G, then your sample could be heterozygous for G->A/C. Similarly for location 6710432.
  Set 2
  read 1 6523432: A - 6710432 T
  read 2 6523432: A - 6710432 T
  read 3 6523432: A - 6710432 T
  read 4 6523432: A - 6710432 T
  
  In this case, if the reference has an A at position 6523432, the reads do not support a variant. However, if the reference has anything else, then the reads would support a homozygous variant e.g. G->A. And the same for position 6710432.
  
  Now, in all of these cases, you cannot know whether a variant is real without taking into account the number of reads that support it, the quality of the alignment, and the quality of the reads. That's why we do this using dedicated variant callers, software that is designed to take all these factors into account and decide whether to call a variant at that position. 
  But in all cases, the variant must be defined with respect to a reference genome. You cannot call a variant just by looking at reads without any knowledge of a reference. A "variant" just means "not what is found in the reference". It is largely a matter of random chance whether something is a "variant" or "normal". If the individual(s) whose genome was taken as a reference has the same sequence, then you call it normal. If not, then you call it a variant. There are many cases where the standard hg19 or hg38 genomes have a "normal" residue which is actually far less common than the known "variant". There's nothing special about the reference genome, after all, it's just the one we chose to call a reference.
  
original_text: 2519: <p>Let's say I'm sequencing a donor, and I have 20X coverage. And let's say all the reads covering a particular base has the same allele at the position. Then, is that position considered a "variant" with only 1 allele in that position or is that not even considered a variant?</p>
  
  <p>The reason I'm asking is, I was given some read data, and it's in the form location:variant - location:variant.</p>
  
  <p>But there was often the case where for all the variants at a particular location, there was only one allele. But this data is labelled to only have variants. Hence, how is a particular location deemed a "variant."</p>
  
  <p>For example, lets say I'm told that the following is <strong>"variant reads given in the form read# location:variant - location:variant"</strong></p>
  
  <pre><code>read 1 6523432: A - 6710432 T
  read 2 6523432: C - 6710432 G
  read 3 6523432: A - 6710432 T
  read 4 6523432: C - 6710432 G
  </code></pre>
  
  <p>Then, why are these two locations considered variants? It seems like it's just a heterozygous locations. Why were these locations considered to carry variants? Similarly, </p>
  
  <pre><code>read 1 6523432: A - 6710432 T
  read 2 6523432: A - 6710432 T
  read 3 6523432: A - 6710432 T
  read 4 6523432: A - 6710432 T
  </code></pre>
  
  <p>In those reads, it just seems like it's homozygous. Why is it being called a variant? How was it classified as such?</p>
  
  <p>Put another way, if I'm given 100000 reads, and I said I only want to look at <strong>informative reads</strong>, and so I'm given a subset with only 100 reads, and those 100 reads look something like the above, then how were these reads selected?</p>
   298: <ol>
  <li><p>Set 1</p>
  
  <pre><code>read 1 6523432: A - 6710432 T
  read 2 6523432: C - 6710432 G
  read 3 6523432: A - 6710432 T
  read 4 6523432: C - 6710432 G
  </code></pre>
  
  <p>I understand this as showing the location of 2 variants, supported by 2 reads each. Now, if the reference has an A at position 6523432, your sample could be heterozygous for an A->C variant at position 6523432. If the reference has a C at that position, then your sample could be heterozygous for C->A. If the reference has something else again, for example G, then your sample could be heterozygous for G->A/C. Similarly for location 6710432.</p></li>
  <li><p>Set 2</p>
  
  <pre><code>read 1 6523432: A - 6710432 T
  read 2 6523432: A - 6710432 T
  read 3 6523432: A - 6710432 T
  read 4 6523432: A - 6710432 T
  </code></pre>
  
  <p>In this case, if the reference has an A at position 6523432, the reads do not support a variant. However, if the reference has anything else, then the reads would support a homozygous variant e.g. G->A. And the same for position 6710432.</p></li>
  </ol>
  
  <p>Now, in all of these cases, you cannot know whether a variant is real without taking into account the number of reads that support it, the quality of the alignment, and the quality of the reads. That's why we do this using dedicated variant callers, software that is designed to take all these factors into account and decide whether to call a variant at that position. </p>
  
  <p>But in all cases, the variant must be defined with respect to a reference genome. You cannot call a variant just by looking at reads without any knowledge of a reference. A "variant" just means "not what is found in the reference". It is largely a matter of random chance whether something is a "variant" or "normal". If the individual(s) whose genome was taken as a reference has the same sequence, then you call it normal. If not, then you call it a variant. There are many cases where the standard hg19 or hg38 genomes have a "normal" residue which is actually far less common than the known "variant". There's nothing special about the reference genome, after all, it's just the one we chose to call a reference.</p>
  