I have a vcf file of SNPs and another vcf file for my Indels. During SNP calling step by GATK, I created them separately instead of using -glm (calling together SNPs and Indels). Now to get the consensus sequence of my mapped genome, I want to put them together in the same vcf file. Should I merge these two vcf files or concatenate them to get the proper variant vcf file, so I can consider SNPs and Indels for my consensus sequence?
•
written
by
bioinfo
2 answers
You should merge them using the vcf-merge utility that is part of the vcftools package.
Concatenation would be appropriate if you had separate files for each chromosome, and simply wanted to join them 'end-to-end' into a single file. In this case, your SNPs and indels need to be inter-woven ie. merged.
If you haven't used vcftools before, you can find it here:
http://vcftools.sourceforge.net/
Specifically, you can read about vcf-merge here:
http://vcftools.sourceforge.net/docs.html#merge
You will need to compress the vcf files with bgzip and index with tabix before you can run any vcftools functions on them.
What should I do to merge VCF files column wise?
The format is as follows
CHROM POSID REF ALT QUAL FILTER INFO FORMAT CAST_EiJ
I want to add a column extra in this, so that my format will be
CHROM POSID REF ALT QUAL FILTER INFO FORMAT CAST_EiJ C57BL6J
in which all the columns except last will be common.
Thanks.
You want to merge, you might also want to take a look at bcftools merge. bcftools is usually pretty quick. I am not sure if this is a requirement but I always make sure the files are sorted, bgzip compressed and indexed, ideally you do not want any overlapping sample IDs either.
Depending on how large your files are, you might want to take the first few rows from each one and use them to test the various options to make sure you are happy with the result, before going onto a time consuming big merge.
bcftools merge expects different samples in VCF files when merge is used, else one needs to use --force-samples
Log in to answer this question.
-
merge vcf files with overlapping positions and create consensus genotype
written by hans.recknagel •
-
Use of GenotypeGVCFs in population genetic studies
written by kk.mahsa
-
Remapping genomic coordinates to account for indels
written by jergosh
-
Merge SNP & INDEL vcf files
written by mplace •
-
SNP calling from a DeNovo genome assembly
written by harishk0201
-
Which contigs formed the consensus sequence
written by low.weiquan •
-
Artemis Reads The Snps But Not The Indel Vcf File
written by bioinfo
-
Checking Frameshifts And Stop Codons In Gene
written by bioinfo
-
Gatk Variant Discovery And Filtering Process For 500K Target Regions
written by Nino •
-
Gatk Multi-Sample Calling
written by Bioscientist
I have two folder snp and indel , both folder have vcf file chromosome wise 1...22 , I want to make one vcf file of all vcf file of both folder ....so
How to use vcftool perl script...
When I run vcf-conct command (vcftools) then occur error "Empty fields in the header line, the column 106 is empty, removing." line ....