vcftools QC process
I am going to have QC process of every CHR1 that I merge. I want to do HWE and call rate. Thant's waht I did with HEW process ./vcftools --vcf life_chr1.vcf --hwe 5e-8 --recode --out output_HWE
Is this right? the output_HWE is the final data??
How can I do with call rate process??(missing ratio)
Thank you very much!
• 2,220 views
•
link
written
by
chengyi31000 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
VCFTools MAF Filter needs Efficiency Inquiry
written by S •
-
Apply hard-call-threshold to pgen files using PLINK2
written by xuan.zhou2017 •
-
Merging multiple .ped/.map files for PLINK analisys purposes
written by giusdalt95 •
-
calculate ti/tv ratio in WGS data
written by chengyi31000 •
-
Call missing variants in VCF as reference allele
written by olavur
-
sequences data QC
written by chengyi31000 •
-
how to cut the gz. file by chromosome
written by chengyi31000 •
-
vcftools flags --mac and --max-mac giving an empty output file
written by vasilipankratov •
-
How to handle missing values when combining multiple VCF files?
written by Ruyang •
-
vcftools subset data missing all the information in FORMAT column
written by MAPK
