Answers.Body: <p>The alternative deletion representation that you are describing, i.e. using <code>.</code> to represent alternate alleles that are not present, is allowed within the definition of the VCF format (in other words, this representation can be found in a standard format). An example of this is provided on the <a href="https://en.wikipedia.org/wiki/Variant_Call_Format" rel="nofollow noreferrer">Wikipedia page for VCF</a>:</p>
  <pre><code>#CHROM POS      ID         REF   ALT    QUAL  FILTER
  20     1230237  .          T     .      47    PASS
  </code></pre>
  <p>While missing values are allowed for <em>ALT</em>, they are not typically used because they lead to ambiguous variant definitions.</p>
  <p>For missing <em>REF</em> values, the <a href="https://samtools.github.io/hts-specs/VCFv4.2.pdf" rel="nofollow noreferrer">VCF format specification</a> has the following explanation:</p>
  <blockquote>
  <p>Each base must be one of A,C,G,T,N (case insensitive). Multiple bases are permitted. The value in the POS field refers to the position of the first base in the String. For simple insertions and deletions in which either the REF or one of the ALT alleles would otherwise be null/empty, the REF and ALT Strings must include the base before the event (which must be reflected in the POS field), unless the event
  occurs at position 1 on the contig in which case it must include the base after the event....</p>
  </blockquote>
  <p>To work around this ambiguity, VCF files are commonly left-normalised so that the specified sequences and variants can be precisely placed within the reference genome. Left normalisation can be carried out using <a href="https://samtools.github.io/bcftools/bcftools.html#norm" rel="nofollow noreferrer"><code>bcftools norm</code></a>.</p>
  <p>More details on left normalisation can be found here:</p>
  <p><a href="https://genome.sph.umich.edu/wiki/Variant_Normalization#Left_alignment" rel="nofollow noreferrer">https://genome.sph.umich.edu/wiki/Variant_Normalization#Left_alignment</a></p>
  <blockquote>
  <ul>
  <li>A variant is normalized if and only if it is parsimonious and left aligned.</li>
  <li>A variant is parsimonious if and only if it is represented in as few nucleotides as possible
  without an allele of length 0.</li>
  <li>A variant is left aligned if and only if it is no longer possible to shift its position to the left while keeping the length of all its alleles constant.</li>
  </ul>
  </blockquote>
  <p>In other words, these are <em><strong>non-parsimonious variants</strong></em>.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2022-10-09T11:35:50.587Z
Answers.LastActivityDate: 2022-10-10T11:50:43.680Z
Answers.LastEditDate: 2022-10-10T11:50:43.680Z
Answers.UserDisplayName: gringer
Body: <p>I have recently been working with variants which are represented in a non-standard format. This alternative format uses a period instead of prefixing variants with a reference base. What is this format called?</p>
  <p>There is no specific output format associated with it. It's part of some calling code I've inherited. I figured that this existed elsewhere. But, I could be wrong.</p>
  <p><strong>Expected Format:</strong></p>
  <p>Insertion:</p>
  <pre><code>pos=4
  alt=ATTT
  ref=A
  </code></pre>
  <p>Deletion:</p>
  <pre><code>pos=4
  alt=A
  ref=ATTT
  </code></pre>
  <p><strong>Alternative Format:</strong></p>
  <p>Insertion:</p>
  <pre><code>pos=5
  alt=TTT
  ref=.
  </code></pre>
  <p>Deletion:</p>
  <pre><code>pos=5
  alt=.
  ref=TTT
  </code></pre>
  
Tags: <variant-calling><variants>
Title: Name this alternative INDEL format
ThreadText: # Name this alternative INDEL format
  
  (Asked by: username_0 on 2022-10-07)
  
  I have recently been working with variants which are represented in a non-standard format. This alternative format uses a period instead of prefixing variants with a reference base. What is this format called?
  
  There is no specific output format associated with it. It's part of some calling code I've inherited. I figured that this existed elsewhere. But, I could be wrong.
  
  **Expected Format:**
  
  Insertion:
  
  ```
  pos=4
  alt=ATTT
  ref=A
  
  ```
  
  Deletion:
  
  ```
  pos=4
  alt=A
  ref=ATTT
  
  ```
  
  **Alternative Format:**
  
  Insertion:
  
  ```
  pos=5
  alt=TTT
  ref=.
  
  ```
  
  Deletion:
  
  ```
  pos=5
  alt=.
  ref=TTT
  
  ```
  
  **Tags:** variant-calling, variants
  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_1 on 2022-10-09. Score: 1 | Accepted answer
  
  The alternative deletion representation that you are describing, i.e. using `.` to represent alternate alleles that are not present, is allowed within the definition of the VCF format (in other words, this representation can be found in a standard format). An example of this is provided on the [Wikipedia page for VCF](https://en.wikipedia.org/wiki/Variant_Call_Format):
  
  ```
  #CHROM POS      ID         REF   ALT    QUAL  FILTER
  20     1230237  .          T     .      47    PASS
  
  ```
  
  While missing values are allowed for *ALT*, they are not typically used because they lead to ambiguous variant definitions.
  
  For missing *REF* values, the [VCF format specification](https://samtools.github.io/hts-specs/VCFv4.2.pdf) has the following explanation:
  
  > Each base must be one of A,C,G,T,N (case insensitive). Multiple bases are permitted. The value in the POS field refers to the position of the first base in the String. For simple insertions and deletions in which either the REF or one of the ALT alleles would otherwise be null/empty, the REF and ALT Strings must include the base before the event (which must be reflected in the POS field), unless the event
  > occurs at position 1 on the contig in which case it must include the base after the event....
  
  To work around this ambiguity, VCF files are commonly left-normalised so that the specified sequences and variants can be precisely placed within the reference genome. Left normalisation can be carried out using [`bcftools norm`](https://samtools.github.io/bcftools/bcftools.html#norm).
  
  More details on left normalisation can be found here:
  
  <https://genome.sph.umich.edu/wiki/Variant_Normalization#Left_alignment>
  
  > * A variant is normalized if and only if it is parsimonious and left aligned.
  > * A variant is parsimonious if and only if it is represented in as few nucleotides as possible
  >   without an allele of length 0.
  > * A variant is left aligned if and only if it is no longer possible to shift its position to the left while keeping the length of all its alleles constant.
  
  In other words, these are ***non-parsimonious variants***.
  
  
  ---
  
  