Question: Using SnpSift to filter VCF for protein_coding regions not working
gravatar for goatsrunfaster
11 months ago by
goatsrunfaster20 wrote:

I have annotated a VCF file using snpEff which has worked wonderfully. Now I am trying to filter that file based on protein-coding regions but just get an outfile file with nothing in it, and the error below:

java -Xmx4g -jar /network/rit/home/sb939359/turnerlab/bin/snpEff/SnpSift.jar filter "ANN[*].EFFECT has 'protein_coding'" all_no_outgroups.recode.ann.vcf > codons.vcf

Exception in thread "main" java.lang.RuntimeException: Error reading file 'all_no_outgroups.recode.ann.vcf'. Line ignored: Line (0): 'null' at org.snpeff.fileIterator.VcfFileIterator.readNext( at org.snpeff.fileIterator.VcfFileIterator.readNext( at org.snpeff.fileIterator.FileIterator.hasNext( at at at at org.snpsift.SnpSift.main(

Any help is appreciated.

snp next-gen genome • 208 views
ADD COMMENTlink written 11 months ago by goatsrunfaster20

is it really a vcf file ? can you read it with bcftools view all_no_outgroups.recode.ann.vcf

ADD REPLYlink written 11 months ago by Pierre Lindenbaum131k

Apparently it is not. When I view the file before processing with snpEff using bcftools there is no problem. When I open the new annotated file (which I ran through snpEff) in bcftools I get the message: input/output error.

Note that viewing the .vcf in a text editor they look identical with the exception of the added ANN field.

What is going on here?

ADD REPLYlink written 11 months ago by goatsrunfaster20

This was my fault, the file became corrupt or something. I re-did the original filtering and then SnpSift worked. thanks!


ADD REPLYlink written 11 months ago by goatsrunfaster20
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 853 users visited in the last hour