This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Identification of SNPs located on genes

Hi guys,

I have a VCF file with SNPs located on the whole genome and I want to identify those that are located on genes. I have converted VCF file to Genotype matrix and wrote a python script that goes through both GM file and GTF one to find SNPs located in genes. But this takes a huge amount of time. Does somebody know if there exists a standard tool that solves my problem in an optimal time-efficient way?

Thanks!

snp vcf
• 2,162 views • link
updated by WouterDeCoster WouterDeCoster • written by valerie valerie

1 answer

If you have a VCF file then you can easily use SnpEff, Annovar or VEP for annotation.

• 0 views • link
written by WouterDeCoster WouterDeCoster

Thank you for your advice! Actually my goal is to study genetic differences between two groups of mice (3 mice in each group). My VCF file contains information about all the samples. After converting to genotype matrix I selected only those SNPs that differ two groups. And now I want to understand what genes are involved in this differences. Is it possible to use the tools you mentioned to solve this problem?

• 0 views • link
written by valerie valerie

Most convenient would probably be to filter your vcf now with what you found when selecting variants differing between the two groups.

• 0 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this