Minimum # markers per chromosome/scaffold for imputation with Beagle 5.1
1
0
Entering edit mode
4.0 years ago

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?!

snp imputation beagle vcf • 952 views
ADD COMMENT
0
Entering edit mode
4.0 years ago

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!

ADD COMMENT

Login before adding your answer.

Traffic: 1485 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6