How to annotate only closest SNP using SnpEff?
0
0
Entering edit mode
4.8 years ago
rimgubaev ▴ 330

Hello everyone! I am trying to annotate only the closest SNP using SnpEff software with the following command:

java -Xmx16G -jar snpEff.jar closest organismDB my.vcf > my.annot.closest.vcf

However when I then try to extract the fields into a table using such command:

cat my.annot.closest.vcf | ./scripts/vcfEffOnePerLine.pl | java -jar  SnpSift.jar extractFields - CHROM POS ID AF "CLOSEST[*]GENEID" > my.annot.closest.vcf.txt

I am getting a lot of columns instead of the five: CHROM, POS, ID, AF, GENEID. If anyone faced such problem please help!

SnpEff annotation SNP • 1.0k views
ADD COMMENT

Login before adding your answer.

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