This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to Generate VCF from two FASTA genomes?

I am wondering if there is a way/tool to align two multi-chromosome genomes in FASTA format and call variants into a VCF. I am trying to determine variants between two strains of Plasmodium falciparum (malaria parasite), which should be vastly identical with relatively few variants, and we have assembled genomes for both stains.

fasta variant alignment vcf

You could use minimap2 to align the data and then call variants. If chromosomes are well defined in each genome then you could do one to one alignments in order to prevent all sorts of secondary alignments elsewhere,

0 answers

No answers yet.

Log in to answer this question.