This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get rsid during the variant annotation using VEP (offline mode)

Dear all,

I've got a large list of variant as VCF format per chromosome (human), which annotated using VEP (release 98) as offline mode. But, I didn't get the corresponding rsid at the output. Could you please help me out how I can obtain the "rsid" during the annotation?

Thanks

vep ensembl annotation
• 5,339 views • link
updated by Emily Emily • written by seta seta

1 answer

--check_existing=1

Checks for the existence of known variants that are co-located with your input. By default the alleles are compared and variants on an allele-specific basis - to compare only coordinates, use --no_check_alleles.

Some databases may contain variants with unknown (null) alleles and these are included by default; to exclude them use --exclude_null_alleles.

See this page for more details.

Not used by default

• 0 views • link
updated by zx8754 zx8754 • written by Emily Emily

Thanks Emily; I annotated VCF file of chromosome 22 and the "Existing_variation" column of output is empty, so what's happened in your view?

• 0 views • link
written by seta seta

Please share exact code you used.

• 0 views • link
written by ATpoint ATpoint

Then maybe you have novel variants.

• 0 views • link
written by Emily Emily

it's my mistake; actually, I thought the "check_existing" is considered by VEP as default. Now, with adding this item, the issue was solved.

• 0 views • link
written by seta seta

If issue is resolved, please accept as an answer.

• 0 views • link
written by zx8754 zx8754

Log in to answer this question.

More posts like this