locateVariants returns an empty result
Hi! I'm having trouble using VariantAnnotation. When I try to assign variants (be it somatic variants or indels) to genes I get the following error:
all <- locateVariants(common[[1]], txdb, AllVariants())
Warning message:
In .returnEmpty() : none of seqlevels(query) match seqlevels(subject)
What does this error mean? I have about 13000 indels and a similar number of somatic variants, so it's very unlikely none are in gene regions. Using ht19 for txdb. Thanks for any tips!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Tutorial: How to align RNA and DNA 3D structures using ARTEMIS
written by Eugene F. Baulin •
-
Issue with VCF format while using Pharmcat
written by taniamahmood38 •
-
Getting more hits running locateVariant than actual number of variants
written by ali.movasati76 •
-
rDGIdb - VCF input error using locateVariant command (VariantAnnotation)
written by tarapcsakszabolcs •
-
SciClone on WES data
written by mar.ark.parr •
-
VariantAnnotation parallel usage
written by russhh
-
Mouse Genome Build And Gtf Conflict In Tophat Downloads
written by k.nirmalraman
-
Tutorial: Machine Learning For Prediction of Relapse in Cancer - Part 2 - Building A Random Forest …
written by Nicholas Spies
