Change VCF File Reference Genome to a Local .FA File
0
0
Entering edit mode
18 months ago
Arda • 0

I have a single chromosome VCF file that I want to align to a chromosome sequence I have locally in a .fa file, is it possible to convert the VCF file so its reference genome is the sequence I have on my .fa file. If so, how could this be achieved.

VCF fasta • 694 views
ADD COMMENT
0
Entering edit mode

Not completely certain what you are asking but this will likely help: https://gatk.broadinstitute.org/hc/en-us/articles/360037594571-FastaAlternateReferenceMaker

ADD REPLY
0
Entering edit mode

So you mean that the VCF was generated by alignment to another reference genome, and now you want to update the coordinates to a new reference genome? If so, the command you are searching for is liftover, available from Picard (and others)

ADD REPLY

Login before adding your answer.

Traffic: 1595 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6