Locus in my snps.vcf file
1
0
Entering edit mode
8.3 years ago
abrahamthz • 0

Hello sorry to trouble you, but I am working with Arabidopsis and In my final file result using GATK I have this information from snps in my sample:

#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  PA01                                                            
1   32210   0   T   C   13.22   LowQual AC=1    AF=1.00 AN=1    BaseQRankSum=-0.389 ClippingRankSum=0.000   DP=11   FS=7.404    MLEAC=1 MLEAF=1.00  MQ=60.00    MQRankSum=0.000 QD=1.20 ReadPosRankSum=-0.275   SOR=3.352   GT:AD:DP:GQ:PL  1:5 6:11:43:43  0

Do you know how I can get the information of gene description or localization from those snps? (For example AT5G28540)

Thanks in advance

Abraham

SNP GATK • 2.3k views
ADD COMMENT
1
Entering edit mode

Have you tried to annotate your snps using a tool like SnpEff? It annotates and predicts the effects of variants on genes (and of course, you will get the gene name).

ADD REPLY
0
Entering edit mode
8.3 years ago

Submit your VCF file to ANNOVAR, it will provide all the annotation you need:

http://wannovar.usc.edu/

ADD COMMENT
1
Entering edit mode

WebAnnovar is limited to human data and therefore unsuitable for the OP's Arabidopsis sample.

ADD REPLY
0
Entering edit mode

Ah yes that is also another thing. In that case VEP will be a better resource or command line annovar which comes with a bundle for species if I reckon correctly

ADD REPLY
0
Entering edit mode

That is true, my mistake.

ADD REPLY
0
Entering edit mode

webannovar has a restriction to the number of variants that it takes as input in a sinlge vcf file as far as I remember, so it depends on the total number of variants in the file of OP, if its less than 500 then web annovar is fine else use the command line , or VEP or SNPeff as suggested by @iarun

ADD REPLY

Login before adding your answer.

Traffic: 883 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6