Dear Amar thanks for response. i hope it will help.
• 0 views
•
link
written
by
ahmadpir40 •
how to merge individual vcf to reference genome using python script.
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.
• 1,464 views
•
link
updated
by
Mark
•
written
by
ahmadpir40 •
•
written
by
ahmadpir40 •
1 answer
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.
Log in to answer this question.
More posts like this
-
How to add reference as new sample to vcf?
written by 沛煒 •
-
mauve for contigs
written by rthapa •
-
Identifying unique SNPs between species
written by Wilber0x •
-
How to write Python script locally?
written by onemore100iq •
-
How to extract gene names for a metabolic pathway from KEGG?
written by bef1 •
-
Merging BAM file before running GATK or merging vcf output
written by pbigbig
-
Python script to find error rate
written by zaidnab •
-
Best way to create a multi-sample VCF
written by Michael •
-
How to get individual chromosome sequence in fasta format from vcf.gz and its vcf.gz.tbi file of 10…
written by David Simmons •
-
check the variations incorporated to consensus fasta by vcf-consensus
written by iraun
