Question: ANNOVAR annotate cuts short the VCF file ?
0
gravatar for BC357
3.8 years ago by
BC35770
United States
BC35770 wrote:

I uploaded a standard VCF that contains 57651 lines (or 57471 lines variant data) and ran it with ANNOVAR , under usegalaxy.org, to annotate the variants. To my surprise, the output from ANNOVAR contains only 24072 variants. When I aligned the two tables, it's clear ANNOVAR "removed" some of the variants (both novel and known SNPs, granted ANNOVAR on Galaxy only uses snp137) in the original VCF.

Is there some default setting in ANNOVAR that makes some cutoff for variant annotation ? I don't see how I can change any parameter when running ANNOVAR on usegalaxy. Any explanation for the short list of ourput or any suggestion what I could do to re-annotate the VCF? Great many thanks.

annovar • 2.1k views
ADD COMMENTlink modified 3.8 years ago • written 3.8 years ago by BC35770
0
gravatar for Jennifer Hillman Jackson
3.8 years ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

Does the .vcf dataset contain data from multiple samples? If so, this could be the reason for the problem. Do you see this in the info comments of the output?

WARNING to old ANNOVAR users: this program no longer does line-to-line conversion for multi-sample VCF files. If you want to include all variants in output, use '-format vcf4old' instead.

Our team is reviewing the tool wrapper for an update. No ticket has been created yet, but I can provide more feedback next week, if this is your case (please confirm yes/no).

I am assuming that you are using the tool on the public Main Galaxy instance at http://usegalaxy.org. If not, can you reproduce the error on Main? We may ask for you to share the history from Main for troubleshooting if this is not the issue and it is reproducible there (privately, not here in Biostars).

Thanks, Jen, Galaxy team

ADD COMMENTlink written 3.8 years ago by Jennifer Hillman Jackson25k

comment cancelled and re-posted as an answer.

ADD REPLYlink modified 3.8 years ago • written 3.8 years ago by BC35770

I am having the same problem. Is there a solution to the problem yet?

ADD REPLYlink written 3.4 years ago by hadaarjan0
0
gravatar for BC357
3.8 years ago by
BC35770
United States
BC35770 wrote:

Sorry for missing your feedback. Indeed, I was using the tool on the public Main Galaxy and did notice the warning of ANNOVAR no longer does line-to-line conversion. From your response, it sounds like the team is reviewing an update to address the issue of ANNOVAR. If correct, is there any fix/solution to this issue of ANNOVAR not doing "line-by-line annotation" during the interim? Thank you.
 

ADD COMMENTlink modified 3.8 years ago • written 3.8 years ago by BC35770

Hello,

You could split the vcf file by sample, which would be the easiest thing to do. You can always merge again after. 

The other option is to update the tool wrapper itself. This would need testing, etc. The comments from the binary are for command-line execution.

I made a ticket for the upgrade request. There is no time-line yet. https://trello.com/c/9puPsRdq

Best, Jen, Galaxy team

ADD REPLYlink written 3.8 years ago by Jennifer Hillman Jackson25k

I split the vcf file by sample and used it as input for ANNOVAR. It worked fine. No warnings. How do I merge the output files for each sample from ANNOVAR?

ADD REPLYlink written 3.4 years ago by hadaarjan0

To merge VCF datasets, see the tool "VCFcombine".

To merge BAM/SAM, Bedgraph, and other formats (tabular, interval, bed, gtf), perform a tool search using the keyword "merge". 

Best, Jen, Galaxy team

ADD REPLYlink written 3.4 years ago by Jennifer Hillman Jackson25k

ANNOVAR gives me tabular data. So I have three tabular data that I need to merge. Is there a right tool on Galaxy main for merging 3 tabular data?

If not can you please point me in right direction?

Thank you for the response.

ADD REPLYlink written 3.4 years ago by hadaarjan0

I have the same problem. Can anyone offer a solution? TIA

ADD REPLYlink written 3.3 years ago by pastpeak0
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 16.09
Traffic: 175 users visited in the last hour