Jit
Enthusiast computational biologist with business dream !!!
Tag cloud
Bioinformatics, Computational Biology, Genome, Perl, PostDoc, Education, Journals, Linux, Install, Research, PhD, Modules, CPAN, NGS, Biology, Ancestral, Gene, sequence, reconstruction, Evolution, List, India, Tools, Ubuntu, Assembly, Data, Interview, Job, BioPerl, Yum, Packages, Quick, Perl Modules, GD, Circos, Evolution Highway, MAc Os, MS Window, Enrichments, Pathway, Annotation, Gene Ontology, GO, Ontology, Respect, Work culture, Life, Strategy, Reasons, Deadline
Our Sponsors
Converting a VCF into a FASTA given some reference !
By Jit 2108 days ago
Samtools/BCFtools (Heng Li) provides a Perl script vcfutils.pl which does this, the function vcf2fq (lines 469-528)
This script has been modified by others to convert InDels as well, e.g. this by David Eccles
./vcf2fq.pl -f <input.fasta> <all-site.vcf> > <output.fastq>https://github.com/gringer/bioinfscripts/blob/master/vcf2fq.pl
https://github.com/lh3/samtools/blob/master/bcftools/vcfutils.pl
