Hello
When merging a list of bed files in plink I get an error
Error: 1785 variants with 3+ alleles present.
The files were originally in vcf format and were converted to bed files using vcftools (with the --plink flag). Is there a simple way to scan through the bed files with plink or another tool and remove them?
•
written
by
jan.breitling •
Not sure what you mean exactly. But have you tried
--biallelic-only strict
in plink 1.90_beta_3o? One can find further info on page https://www.cog-genomics.org/plink2/input
•
written
by
fengchongwang •
--biallelic-only is retired. Use e.g. --max-alleles 2 instead
•
written
by
Shicheng Guo
is there no way to do it in the vcf sry ped file..otherwise is has to be reconverted before...
If you want to keep some or all of those variants, and then treat least common alternate allele calls as missing, you should perform the merge with another tool, and then use plink --vcf to import from the merged VCF.
However, if you just want to get rid of all the triallelic variants, refer to the last batch of sample commands under https://www.cog-genomics.org/plink2/data#merge3. The .missnp file generated during the failed merge is designed to be used with --exclude.
•
written
by
chrchang523
I am merging a whole list of files and the problem is that I don't know from which subfile the variants in the mssnp file you mention come from. Is there a way to use this field to exclude the snps still?
•
written
by
jan.breitling •
Yes, it's safe to --exclude [prefix].missnp on every single fileset. Nothing bad happens if a variant named in the .missnp file is not in the current fileset.
•
written
by
chrchang523
Log in to answer this question.
