Your command should be
plink --allow-no-sex --23file jovialisrawdnafile.txt --make-bed --out jovialis
After that you need to rename your sample which will be called ID001 in the .fam file (use a code editor) then filter.
23ame V3 format has the highest coverage, you should use that or the all_SNPs option, the latter has some glitches I've noticed so v3 is the best (imo).
WGSExtract can produce a raw dna file with all SNPs from hg38.
But, if you want to use the vcf you can use DNA kit studio that will convert it to 23afile txt that plink can read.
If you want some feedback, you are wasting your time.
The highest SNP count possible you'll get with the Reich data is 1.24mil SNPs. You can hit that number easy if you extract a raw file from your BAM file.
Then you can convert with plink.
Example, you can't make Turkey_BA work with Yamnaya (outgroup from Southern Arc):
target left weight se z
<chr> <chr> <dbl> <dbl> <dbl>
1 Turkey_EBA Turkey_N 0.425 0.0402 10.6
2 Turkey_EBA CHG...
This site uses cookies to help personalise content, tailor your experience and to keep you logged in if you register.
By continuing to use this site, you are consenting to our use of cookies.