Question: Running vcf2maf on snpeff-annotated vcfs
0
gravatar for becky3.a
19 months ago by
becky3.a0
becky3.a0 wrote:

Hi,

I have vcf files that were annotated using snpeff v4.3 and I would like to convert them to maf using vcf2maf.

My first run included the VEP component (below) but I got odd results when I run MutsigCV using the maf file I had generated using vcf2maf:

perl vcf2maf.pl --input-vcf file.vcf --output-maf file.maf --ref-fasta ~/GRCh37-lite.fa --vep-forks 4 --vep-path ~/vep --ref-fasta $VEP_DATA/homo_sapiens/86_GRCh37/Homo_sapiens.GRCh37.75.dna.primary_assembly.fa.gz --species homo_sapiens --tumor-id tid --normal-id nid

I then tried another run (below) that skips the VEP annotation step:

perl vcf2maf.pl --input-vcf file.vcf --output-maf file.maf --ref-fasta ~/GRCh37-lite.fa --species homo_sapiens --tumor-id ddd --normal-id dddnn

I got the following errors:

Use of uninitialized value in list assignment at ~/vcf2maf.pl line 649, line 2186.
Use of uninitialized value in list assignment at ~/vcf2maf.pl line 649, line 2186.
Use of uninitialized value in list assignment at ~/vcf2maf.pl line 649, line 2186.
Use of uninitialized value in list assignment at ~/vcf2maf.pl line 649, line 2186.
Use of uninitialized value in list assignment at ~/vcf2maf.pl line 649, line 2186.
....

Where could I be going wrong and how best can I go about this?

Thanks

vcf2maf snpeff • 821 views
ADD COMMENTlink modified 18 months ago • written 19 months ago by becky3.a0

Please add more relevant tags, like vcf and vcf2maf. software error conveys next to nothing of value, and better tags mean quicker, more targeted responses.

Also, please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.
code_formatting

ADD REPLYlink written 19 months ago by RamRS25k

To skip the VEP step I read somewhere (probably github) that in the same directory:

cp input.vcf input.vep.vcf

The first step of the script is to annotate using VEP, but if there is already a .vep.vcf file present it skips it.

ADD REPLYlink written 16 months ago by Barry Digby310
0
gravatar for sutturka
19 months ago by
sutturka170
USA
sutturka170 wrote:

This is a known issue and could happen because of newer annotation types added in SnpEff. Please see the GitHub issues page for vcf2maf.

ADD COMMENTlink written 19 months ago by sutturka170

Hi Sutturka, thanks for your feedback. I have been following the vcf2maf GitHub page and issues. I posted an issue with a similar question but I haven't received feedback. Thanks again.

ADD REPLYlink modified 19 months ago • written 19 months ago by becky3.a0

This issue on github page is very similar and comments on this issue may be helpful.

ADD REPLYlink modified 18 months ago • written 18 months ago by sutturka170
0
gravatar for becky3.a
18 months ago by
becky3.a0
becky3.a0 wrote:

Thanks for your response. I think the only way around this is to forget snpeff and use vep. This becomes a problem for teams that receive and work with vcf files already annotated with snpeff.

ADD COMMENTlink written 18 months ago by becky3.a0
Please log in to add an answer.

Help
Access

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