split vcf, choice
Hi,
I have to split a multisample VCF (the output of GATK4) into single vcfs (per sample). What is the best code? Then I have to annotate them with AnnovaR.
vcf-subset -c samplename -e inputvcf > sample.vcf
bcftools view -Oz -s samplename inputvcf > sample.vcf.gz
Is it correct to add the -e option in the vcf-subset code? or I'm going to lose data?
Best
• 1,614 views
•
link
updated
by
likame8586 •
•
written
by
emma.a •
•
written
by
emma.a •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
The allele with index 2 is not defined in the REF/ALT columns" error after AF filtering with gnomAD…
written by erikatatianacs •
-
How to extract and analyze variant sequences from gnomAD v4 exome VCF
written by Rohan •
-
BCFtools consensus
written by Peerzada •
-
Concensus from 1000 genome project
written by Peerzada •
-
for loop in a bash script
written by Erika •
-
How to filter SNP and INDEL separately for VCF files
written by YL •
-
split vcf file into multiple files
written by lz1903 •
-
Add tag to FORMAT tab of vcf
written by cmiliotis •
-
Problem with annovar in GATK4 output
written by ahmad mousavi
-
Split a VCF file into individual sample files
written by win
If that particular sample doesn't contain a variant i.e either not called or "reference homozygous", it will not be printed, which makes sense.