How can we convert a vcf to fasta, so that I can blast some genes against that whole genome sequences?

0

Entering edit mode

6 months ago

Harshita • 0

I want to use my whole genome sequence which is in vcf format to .fa, so that I can find some genes using BLAST.

blast fasta vcf gene • 433 views

ADD COMMENT • link updated 6 months ago by Ram 44k • written 6 months ago by Harshita • 0

2

Entering edit mode

You cannot convert a VCF to a FASTA as they have vastly different information content. Please explain your process - what exactly do you want to "convert"?

ADD REPLY • link 6 months ago by Ram 44k

0

Entering edit mode

well you try to use FastaAlternateReferenceMaker https://gatk.broadinstitute.org/hc/en-us/articles/360037594571-FastaAlternateReferenceMaker but as your VCF was called from a reference genome do you think you'll find something new for a few variants ?

ADD REPLY • link 6 months ago by Pierre Lindenbaum 164k

Login before adding your answer.