Counting of SNPs in different genomic regions

[shis]

Hi,
After annotating the SNP.vcf files using SnpEff v 4.2, I would like to count how many SNPs are found in different genomic regions such as upstream, gene, exon, intron, CDS, 5' UTR, 3' UTR, downstream and intergenic. To perform the SNPs count, I am looking for a program or perl- or python script for few days in different forums. I have annotated SNP.VCF and gene.txt file obtained from SNpEff. Can anyone help me - how can I count the number of SNPs found in different genomic regions? I would highly appreciate your help. Thanks

[gringer]

I've forgotten how SnpEff outputs information, but in general using 'grep -c' is a quick (but mostly manual) way.

If you want to count lots of things at once, isolate out the categories (e.g. using 'cut'), then run 'sort' and 'uniq -c'. Alternatively, import the data into R and use 'table'.

[shis]

Thanks a lot Gringer. I used SnpSift Extract Fields to get the effect of SNPs in different regions. Then the number of SNPs in different regions are counted in Excel. Finally I got SNPs number in different genomic regions.