I have used snpEff to annotate 1000genomes phase 3 vcf files. However, I seem to be getting a large number of warnings that read: 'WARNING_REF_DOES_NOT_MATCH_GENOME'. Too many of these warnings indicate a lack of match up between vcf file and the database snpEff is using.
I have been using the hg38 database, and get this warning in almost 25% of my variants. This surprises me as I am using data direct from 1000 genomes and the latest hg38 database and was hoping someone may shed some light on why this is.
N.b. This is more a matter of interest than anything else, as I have downloaded the gencode annotated vcf files that 1000genomes website provide and can use them instead.
0 answers
No answers yet.
Log in to answer this question.
-
Snpeff run with downloaded genome database
written by ttom
-
Selscan error: Variant physical position must be strictly increasing.
written by spiral01
-
Discrepencies between ANNOVAR and Gencode annotations of 1000 genomes
written by spiral01
-
Annotating Altai Neanderthal and Denisovan variants
written by spiral01
-
Show warnings and errors in snpeff summary
written by misbahabas •
-
Estimating site information from 1000 genomes vcf data
written by spiral01
-
PyVCF error: KeyError: 'ANN' when appending annotation fields to a list
written by spiral01
-
Ensembl 86 annotations with SnpEff for GRCh37
written by raghus606 •
-
Rooting phylogenetic trees in R using Ape's root function: Trees don't seem to root.
written by spiral01
-
Chromosome X variants in 1000genomes Phase3
written by Ram
Check to make sure that
chromosome namesmatch in your VCF files and the genome you are using.