how to create a multifasta file from miltisample vcf file?
Hi everyone. I'd like to create a new fasta file containing all sequences for which I called variants in a vcf file. I can only get one sequence, but I would like to rebuild the sequence for all my samples. I know I can do it individually, but is there any tool to do it for all my samples at once? Thank you very much!
1 answer
Thank you very much! This is what I actually did, one by one within an iteration. But it is a long pipeline (split the vcf into individual files, bzip, tabix, then vcf-consensus for each one and finally cat all fasta in one file). Also, I end up with many unuseful (intermediate) files (I'm analyzing 48 samples). The vcf to MAF maybe can help.
Log in to answer this question.
More posts like this
-
How to subtract variants from one VCF file to another?
written by cee28 •
-
Creating a Variant containing FASTA for proteomics search from VCF and genomic FASTA
written by chscho •
-
Building vcf from most common alleles
written by Eugene •
-
how to get a consensus sequence for each sample in a vcf file using a reference sequence?
written by daniellamodukpe •
-
Consensus protein from aavf
written by pennakiza
-
problem in running vcf-consensus command
written by rajesh.msch •
-
How to find variants that are common in one population but rare in others (population differentiati…
written by SOHAIL
-
Splitting vcf files to individual samples
written by Hadeel •
-
Remove positions that are non-variant in a subset of samples from a vcf file
written by Aliaksei Holik
-
Does vcftools or GATK have an option to skip indels when making alternate sequence?
written by severalorks
How about a simple iteration of vcf-consensus -s <sample_name> . The new version also gives an option to use IUPAC ambiguity codes for hetrozygotes or no-calls.
See this post. It may be at least partially helpful.
vcf to MAF to fasta