Asslamu Alikum
I have vcf file look like
chr2 104 . T <non_ref> . . . T:AD:DP:GQ:PL 0/0:2,0:2:3:0,3,45
chr2 105 . T <non_ref> . . . GT:AD:DP:GQ:PL 0/0:1,1:2:0:0,0,0
chr2 104 . T <non_ref> . . . T:AD:DP:GQ:PL 0/0:2,0:2:3:0,3,45
chr2 105 . T <non_ref>. . . GT:AD:DP:GQ:PL 0/0:1,1:2:0:0,0,0
Is it possible to add nucleotide in ref column?
•
written
by
misbahabas •
1 answer
It is, but it would require providing a fasta file of the reference when initially calling the VCF from a bam file in GATK.
I'm not sure what was used to make your VCF file, but I'm going to hazard a guess that it was HaplotypeCaller, in which case the relevant documentation can be found here (with the reference fasta file fed to the -R argument).
Log in to answer this question.
-
VCFtools doesn't keep any variants using GATK output
written by fcarolinebe •
-
sed (or similar) command to remove part of the FORMAT from a vcf file?
written by Jautis
-
Why invariant blocks in GATK consistently have very low quality scores (but not variant sites)
written by yixianshuiesuan •
-
Vcf produced after ApplyVSQR step
written by BAGeno
-
HOW EDIT VCF FILE COLUMNS
written by misbahabas •
-
Convert specific regions in bam file to vcf and annotating the vcf
written by komal.rathi
-
Calculate Allele balance by sample
written by kirannbishwa01
-
QD value from GATK in joint VCF
written by deepmehra.atgc •
-
problems with MAF for MutSigCV (vcf2maf)
written by A. Domingues
-
Id And Ref Are Empty
written by Zhshqzyc •
