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
17 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 • 666 views
ADD COMMENTlink modified 17 months ago by Emily_Ensembl18k • written 17 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 17 months ago by RamRS23k

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 17 months ago by Pierre Lindenbaum122k

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 16 months ago by Michael Dondrup46k
1
gravatar for Emily_Ensembl
17 months ago by
Emily_Ensembl18k
EMBL-EBI
Emily_Ensembl18k 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 17 months ago • written 17 months ago by Emily_Ensembl18k

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

ADD REPLYlink written 17 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 17 months ago by RamRS23k

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 17 months ago by RamRS23k

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 17 months ago by Emily_Ensembl18k
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: 1052 users visited in the last hour