Question: (Closed) Is it possible to run VEP on a merged VCF file derived from a 30 VCFs from the separate reference gene alignments?
0
gravatar for mathisdead
11 months ago by
mathisdead10
mathisdead10 wrote:

Is it possible to run VEP on a merged VCF file derived from a 30 VCFs from the separate reference gene alignments? And is there any rational way to run VEP on a batch of separate VCFs?

human vep annotation vcf • 475 views
ADD COMMENTlink modified 11 months ago by Emily_Ensembl17k • written 11 months ago by mathisdead10
2

Better title: Running VEP on VCF with samples aligned to separate reference genomes (.... What?)

How did you merge VCFs with samples aligned to different ref genomes? Any proper merge operation would need exactly one ref genome to perform the merge.

ADD REPLYlink written 11 months ago by RamRS20k

please validate / comment your previous question: Are there any human genome variant callers and BAM indexing and sorting tools to run on Windows 10? ;

ADD REPLYlink written 11 months ago by Pierre Lindenbaum116k

Hello mathisdead!

We believe that this post does not fit the main topic of this site.

Not a real question (bot)

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

ADD REPLYlink written 10 months ago by Michael Dondrup45k
1
gravatar for Emily_Ensembl
11 months ago by
Emily_Ensembl17k
EMBL-EBI
Emily_Ensembl17k wrote:

VCF format should not allow you to have different reference genomes. You should retain your 30 different files if they use different references.

Running the VEP is simply a case of running a command on the command line. If you want to run it 30 times and you don't feel like writing our the command over and over, you can just use your favourite scripting language to write a short script that runs through your list and plugs the appropriate values into the standard VEP command.

Edit: added highlighting to the text to make it clear what the OP needs to to next.

ADD COMMENTlink modified 11 months ago • written 11 months ago by Emily_Ensembl17k

thx, show a Python or Ruby code of such a script, pls

ADD REPLYlink written 11 months ago by mathisdead10

Did you try writing a script? Give us an idea of the effort you put in and we'll help you get better at scripting.

ADD REPLYlink written 11 months ago by RamRS20k

Also, don't forget to follow up on your previous threads.

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. Upvote|Bookmark|Accept

ADD REPLYlink written 11 months ago by RamRS20k

I've added highlighting to my post to make it clear what you need to do. If you would like help troubleshooting your script, please post it here and we will be willing to help you.

ADD REPLYlink written 11 months ago by Emily_Ensembl17k
Please log in to add an answer.
The thread is closed. No new answers may be added.

Help
Access

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