Question: how to merge individual vcf to reference genome using python script.
0
gravatar for ahmadpir40
10 months ago by
ahmadpir400
ahmadpir400 wrote:

I need python script to merge the Individual vcf to reference genome to get a new reference. Any suggestion about it how to write it.

ADD COMMENTlink modified 10 months ago by Mark800 • written 10 months ago by ahmadpir400
0
gravatar for Mark
10 months ago by
Mark800
Mark800 wrote:

Super relevant biostars thread that answers this question: New Fasta Sequence From Reference Fasta And Variant Calls File?

Actual answer: Why reinvent the wheel when someone has done the hard work for you? FastaAlternateReferenceMaker is part of GATK, there's vcf-consensus in vcftools, lots of options.

ADD COMMENTlink modified 10 months ago • written 10 months ago by Mark800

Dear Amar thanks for response. i hope it will help.

ADD REPLYlink written 10 months ago by ahmadpir400
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: 1383 users visited in the last hour