Hi Pierre, thanks for your reply. For the first step, which should be the input file? I have checked the 1000 genomes release files, they seems only contain the sequencing information for variant sites (eg SNPs, INDELs, SVS) but no invariant sites or monomorphic sites.
How To Get Non-Snp Sites Or Invariant Sites Based On 1000 Genomes Project
Hi everyone, I need to find the human genomic sites which are invariant among all the individuals (identified by sequencing in all individuals and also without SNP). I know the sites which are not included in the 38 million SNPs identified by the 1000 genome project are possible candidates. However, not all of them are necessarily sequenced or covered in all individuals. One possible way to get such invariant sites is to first find the sites which are covered in all individuals and then subtract the SNP sites. Does anybody know how to get such information or where to get it (exclude extracting such information from the raw mapping data)? Thanks very much.
• 6,082 views
•
link
updated
by
esha.sharma •
•
written
by
Gangcai
•
written
by
Gangcai
2 answers
use GATK SelectVariants to keep the variant where all the samples have been covered; something like (not tested):
-select 'vc.getGenotype("SAMPLE1").isCalled() && vc.getGenotype("SAMPLE2").isCalled() && vc.getGenotype("SAMPLE3").isCalled() .... etc...'
see also GATK multi-sample VCF VariantFiltration
use vcftools to remove the SNP from the VCF of the 1K genomes projects http://vcftools.sourceforge.net/docs.html#isec
• 151 views
•
link
updated
by
Ram
•
written
by
Pierre Lindenbaum
•
written
by
Pierre Lindenbaum
hello Gangcai, did u get any way to get non-SNP or invariant sites based on 1000 genome project ??? I am trying Unified Genotyper of GATK-2.7-2 for the same, will it give any direction for finding invariants in exomes ???
• 0 views
•
link
written
by
esha.sharma •
Log in to answer this question.
More posts like this
-
Bcftools consensus unexpected output
written by selplat21 •
-
Manually add invariant sites from a FASTA to VCF?
written by tal.shalev •
-
How to extract genomic positions of SNPs using snp-sites?
written by pltbiotech_tkarthi
-
Calculation of heterozygosity
written by shinken123
-
Determining efficient solvents for extracting metabolites
written by Natasha •
-
IBS in Plink gives only one cluster ??
written by manon.devisser •
-
Genome sites from BPM file
written by Jautis
-
Viral haplotyping using pair end information
written by lovett05 •
-
Why some SNP databases don't contain the strand information
written by camelbbs
-
High Variation Within Biological Replicates For Significant Genes Detected By Deseq/Edger
written by Gangcai
The vcf file by definition can record monomorphic sites (sites without alternate alleles), however the vcf files from the 1000 genomes release (ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/) seems only contain SNPs, INDELS, SVS but no monomorphic or invariant sites.