Recreating QC of 1000 Genomes project - removing non overlapping SNPs
Hi everyone,
I am attempting to recreate the the quality control analysis performed in the 1000 genomes project (http://tcag.ca/documents/tools/omni25_qcReport.pdf).
I am fairly new to performing QC on a dataset, and am currently stuck on section 5.1 of the analysis. The 1000 genomes data was assembled with 2 different sets of markers, and in the analysis only the overlapping 2,150,028 SNPs in both datasets were used.
Does anyone know how I can remove any of the non-overlapping SNPs? Could this be done in PLINK, and should it be performed on the .vcf file?
• 1,114 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
nick.gold •
•
written
by
nick.gold •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Keep SNP's with call rate greater than 30 in VCF
written by Smilesky •
-
GWASQC
written by karthick •
-
Calculating number of SNPs in linkage disequilibrium at different thresholds from 1000 genomes proj…
written by octopuslegs11 •
-
how many genomes from BLAST hit should be considered for phylogeny tree analysis?
written by DanielC
-
Where to download official reference and alternative alleles from 1000 Genomes?
written by Mr Locuace
-
PLINK error: "Could not place marker for left window edge"
written by ebrudermanver
-
(beginner) Where to find a good tutorial for QC in gwas data, genotype imputation
written by Laura •
-
Significance of PLINK epistasis p-values
written by viavialin •
-
Clarification On P-Value Threshold For Gwas
written by mlscmahe
-
How To Uncompress The 1000 Genome Vcf.Gz File
written by GPR
have a look at
bcftools isec