I'm trying to create a VCF formatted file from a pairwise alignment of two genomes. Downstream I'm going to use this to calibrate some deep-sequencing data where I know the input sequence (and thus the variants).
I'm just running into an odd edge-case in the VCF format. Following the specifications described in http://samtools.github.io/hts-specs/VCFv4.2.pdf I'm running into the edge-case where there is a "gap" in the first base of the query sequence with respect to the reference:
>Ref
ATTTGT
>Query
-TTTGT
Based on my understanding of the the VCF format I need to put the preceding & current base in the "REF" column and the preceding base in the ALT column. The POS should also be 0?
Any suggestions?
0 answers
No answers yet.
Log in to answer this question.
-
Where to find info on VCF file format?
written by _quantum_girl_ •
-
What does the position column in the VCF file from PostprocessGermlineCNVCalls refer to?
written by Frieda •
-
structural variation viewer
written by qqq2395 •
-
About Foreign Fields in VCF (4.3)
written by Ibrahim Tanyalcin
-
ALT <*> Interpretation
written by Ada •
-
How do I embed licence information in a VCF file?
written by hyanwong
-
Valid or invalid vcf? REF nucleotide is not matching reference in vcf lines where ALT equals <CNV>
written by dschika
-
Why do indel variants need a leading base according to the VCF spec?
written by William
-
In VCF files, what character encodings are used / legal?
written by Andrew Thomas •
-
Interpreting Gaps at Pos 0 in Terms of VCF
written by pld
duplicate of Interpreting Gaps at Pos 0 in Terms of VCF
Yup, had the wrong search terms.