How to annotate only closest SNP using SnpEff?
0
0
Entering edit mode
18 months ago
rimgubaev ▴ 220

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 • 437 views
ADD COMMENT

Login before adding your answer.

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