how to extract LOF variants from vcf file?
0
0
Entering edit mode
6.4 years ago
reza ▴ 300

hi all

i have a vcf file resulted from whole genome sequencing project and now i want to extract heterozygote and homozygote loss function (LOF) non synonymous SNPs. how can i do it?

thanks in advance

variants vcf Loss function variants • 2.3k views
ADD COMMENT
1
Entering edit mode

search this site for SNPEff and / or VEP

ADD REPLY
1
Entering edit mode

See the Variant Effect Predictor from Ensemble. An awesome tool, but a bit bulky and tricky at times. Did cost me a lot of nerves but I think well worth it. Here are the docs, come back if you have questions.

ADD REPLY
0
Entering edit mode

i annotated vcf file using snpeff but i do not now how find LOFs

ADD REPLY
0
Entering edit mode

i found that i must add -lof option to my command to vcf annotation as following:

java -jar snpEff.jar eff -v -lof BB_12 WC247_final_snps.vcf > WC247_final_snps_lof.vcf

but there is no LOF tags in INFO fields. what part of my work is fault? it is worth to mention that building of database for annotation done by myself in snpEff.

ADD REPLY

Login before adding your answer.

Traffic: 2015 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