Answers.Body: <p>From the latest <a href="https://samtools.github.io/hts-specs/VCFv4.3.pdf" rel="noreferrer">VCF specs</a> (page 8):</p>
  
  <ol start="4">
  <li>REF - reference base(s): 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;</li>
  </ol>
  
  <p>So I think you're right</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-09-06T14:07:27.647Z
Answers.LastActivityDate: 2017-09-06T14:07:27.647Z
Answers.UserDisplayName: Emily_Ensembl
Body: <p>I am writing a small script to write a VCF file from short multiple sequence alignments, and stumbled into an interesting corner case.</p>
  
  <p>This is <a href="https://www.ncbi.nlm.nih.gov/variation/tools/reporter/docs/examples#vcf-section-3" rel="noreferrer">the way a deletion should be represented in VCF format</a> (I assume version 4+):</p>
  
  <pre><code>NC_000016.9  2138199 rs1799758   GTGAG   G  .   .  RS=1799758;...
  </code></pre>
  
  <p>But what happens if the deletion is at the first position? Immagine the following two aligned sequences:</p>
  
  <pre><code>&gt;reference
  AGTC
  &gt;sample
  -GTC
  </code></pre>
  
  <p>How should it be encoded in VCF format? According to <a href="http://seqanswers.com/forums/showthread.php?t=30884" rel="noreferrer">this thread on SEQAnswers</a>, it might be necessary to include the first position after the end of the deletion, like this:</p>
  
  <pre><code>chromosome  1  AG G  .  .  [...]
  </code></pre>
  
  <p>Is this a VCF 4+ correct way to encode such variants?</p>
  
Tags: <vcf><variation>
Title: How to represent a deletion at position 1 in a VCF file?
ThreadText: # How to represent a deletion at position 1 in a VCF file?
  
  (Asked by: username_0 on 2017-09-06)
  
  I am writing a small script to write a VCF file from short multiple sequence alignments, and stumbled into an interesting corner case.
  
  This is [the way a deletion should be represented in VCF format](https://www.ncbi.nlm.nih.gov/variation/tools/reporter/docs/examples#vcf-section-3) (I assume version 4+):
  
  ```
  NC_000016.9  2138199 rs1799758   GTGAG   G  .   .  RS=1799758;...
  
  ```
  
  But what happens if the deletion is at the first position? Immagine the following two aligned sequences:
  
  ```
  >reference
  AGTC
  >sample
  -GTC
  
  ```
  
  How should it be encoded in VCF format? According to [this thread on SEQAnswers](http://seqanswers.com/forums/showthread.php?t=30884), it might be necessary to include the first position after the end of the deletion, like this:
  
  ```
  chromosome  1  AG G  .  .  [...]
  
  ```
  
  Is this a VCF 4+ correct way to encode such variants?
  
  **Tags:** vcf, variation
  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_1 on 2017-09-06. Score: 7 | Accepted answer
  
  From the latest [VCF specs](https://samtools.github.io/hts-specs/VCFv4.3.pdf) (page 8):
  
  4. REF - reference base(s): 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;
  
  So I think you're right
  
  
  ---
  
  