I have a multisample vcf file from a study species with a pretty fragmented reference genome. I already applied the basic filtering steps to my vcf using vcftools. 1- Sample filtering: filter out low coverage samples 2- Site and genotype filtering with the following vcftools parameters: --minDP 3 --max-missing 0.8 --maf 0.05 I am now trying to run imputation of missing genotypes using BEAGLE 5.1 and it is running smoothly until it reaches a small scaffold that only has 1 position in the vcf. So clearly I need to erase scaffolds with too few positions from my vcf file, but here is my question: how many positions per scaffold is too few? I tried to look into the BEAGLE documentation but I couldn't find an answer, hopefully someone here can help?!
1 answer
Solved it! I removed all scaffolds harbouring less than 3 SNP in my vcf file and it worked smoothly! Good luck to anyone else having the same problem!
Log in to answer this question.
-
Problem after filtering vcf files
written by yoser4 •
-
Which VCF file to select for generating Multiple sequence alignment of a gene?
written by anikcropscience •
-
Impute genotype probabilities with BEAGLE
written by dario.galanti92 •
-
Hardy-Weinberg Equilibrium and MAF filtering post imputation
written by andreas131298 •
-
Beagle 5.2 Imputation Issue, understanding output
written by June •
-
unable to impute with beagle after filtering VCF table using vcftools
written by ziv_attia •
-
--max-missing filter in vcftools removes most of my SNP data
written by truebeliever24 •
-
Filter max-missing multiple time removes more SNPs
written by QPaps04 •
-
Imputation of vcf file
written by mab658
-
vcf SNP file minor allele frequency filter?
written by sun.nation
