Question: How to find rs number of SNP when UCSC data base position known in reverse strand sequence
0
gravatar for Nawaj Shaikh
17 months ago by
National AIDS Research Institute, Pune, India
Nawaj Shaikh0 wrote:

I have screened a human gene to find out its association with infection. For the same i have used UCSC data base as reference sequence. I have noted some intronic and exon base SNPs. To identify its rs number i used Predict SNP online tool but it gives following massage (The user's given refererence alleles of more than 50% input variants do not match the reference alleles of selected genome assembly (GRCh37/hg19 and GRCh38/hg38 supported). Please make sure that you selected correct genome assembly.) When i am choosing GRCh38/hg38 this option for genome assembly, same massage coming. Please anybody tell me How to find rs number of SNP when UCSC data base position known in reverse strand sequence? Is there another online tool use to identify rs number when reference sequence is reveres strand? or any other .... Please....

snp • 704 views
ADD COMMENTlink modified 17 months ago by ATpoint14k • written 17 months ago by Nawaj Shaikh0

Hello,

how does your input looks like? How should your expected output look like?

fin swimmer

ADD REPLYlink written 17 months ago by finswimmer11k

I just need rs number in out put.

ADD REPLYlink written 17 months ago by Nawaj Shaikh0

And your input? Please post an exact example.

ADD REPLYlink written 17 months ago by finswimmer11k

For example, I have screened a gene TLR9 for polymorphism and i found variation (G to T) at chromosome position chr3: 52222681 (UCSC data base) (rev), Hence i have to search what is rs number allotted for the same position. which online tools i can use to identify the same? Kindly guide me on same. Thanking you,

ADD REPLYlink written 17 months ago by Nawaj Shaikh0
0
gravatar for ATpoint
17 months ago by
ATpoint14k
Germany
ATpoint14k wrote:

The user's given refererence alleles of more than 50% input variants do not match the reference alleles of selected genome assembly.

This sounds like a typical 1-based/0-based coordinate error. Check if the online tool expects 0-based (BED-like) coordinates or 1-based (genome browser) coordinates. Which tool did you use for variant calling?

If you called SNVs they are probably in VCF format? If so, you can use the Ensemble Variant Effect Predictor with the --check-existing option. It will give you dbSNP IDs in case of overlaps with your data and optionally a plethora of other annotations, such as allele frequencies in the 1000 Genomes project, polyphen scores etc.

ADD COMMENTlink modified 17 months ago • written 17 months ago by ATpoint14k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1131 users visited in the last hour