Hello,
I have about 4200 array genotyping VCFs (from the Illumina Infinium CoreExome-24 Kit) and I have merged them using bcftools merge. The chip has 500K exonic SNPs. These are trio data - which means 1700 of them are probands, mother and father respectively.
Following that, I converted them to Plink format using plinkv2 --make-bed function. I included the fam information later on by generating a fam file.
Finally, I ran KING --kinship but am getting a very high error rate for within family relatedness, which I know is not correct.
Are there other QC steps I should have been running at the merge step? I have removed all monomorphic SNPs and multi allelic sites, too.
•
written
by
Anupriya •
1 answer
if you already have VCF try https://github.com/brentp/somalier
Log in to answer this question.
-
Validate RNAseq reads against array genotypes?
written by S •
-
Trio Phasing using SHAPEIT2 and BEAGLE4
written by Anupriya Dalmia •
-
KING relatedness inference hanging
written by marcus.tutert •
-
KING struggle: Relatedness
written by German.M.Demidov
-
Merging VCF files into one
written by sara.aglantarek •
-
GWAS in mixed populations
written by USA_225478 •
-
Infer pedigree from Kinship coefficients
written by landscape95
-
Trio SNP data to test script (preferably from Infinium Global Screening Array)
written by kurt.cagri •
-
Merging a PLINK file with itself after duplicating sampleIDs and rsIDs
written by aritra90
-
Trio-based analysis: pipeline for phasing the data
written by newbee •
