Change VCF File Reference Genome to a Local .FA File
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.
• 1,206 views
•
link
updated
by
WouterDeCoster
•
written
by
Arda •
•
written
by
Arda •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Percentage coverage of reference genome by de novo genome
written by Lemonhope •
-
Are we convert U to T in RNA before alignment?
written by ธาดา •
-
Convert multi amino acid sequence file (.faa) to nucleotide sequence file (multifasta; .fa)?
written by goatsrunfaster •
-
hisat2-build multiple fasta and Annotation
written by chimerajit •
-
How to convert .tar file to fasta file?
written by jabbari.parnian •
-
Bowtie alignment using RY format as reference
written by Nalini •
-
BAM/SAM to aligned FASTA file
written by ryan.englander •
-
Speeding up NUCmer for de novo contig to reference alignment
written by lkw222 •
-
VCF file alignment to the reference genome
written by Sathish •
-
How To Create A Bigwig File And Display It
written by rig •
Not completely certain what you are asking but this will likely help: https://gatk.broadinstitute.org/hc/en-us/articles/360037594571-FastaAlternateReferenceMaker
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)