This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Create Vcf File From A Multiple Sequence Alignments

Dear Biostars,
I have a question concerning the generation of vcf (variant calling format) creation.
Does anyone know of a tool that would allow me to turn a multiple sequence alignment (containing reference and several variants) into a vcf file?
thanks!

EDIT:
I have a multiple sequence alignment of a several cloned papillomaviruses. We know that the sequence of each individual genome are correct. I.e. all variations between the reference and these additional sequences represent naturally occurring SNPs (and not sequencing errors). I would like to extract the SNPs (and indels) from this alignment and create a vcf file. I hope this clarifies the problem! thanks again

alignment vcf variant snp
• 19,823 views • link
updated by Patrick • Patrick • written by Whetting Whetting

SNP calling is a little bit different from what I am looking for. Calling implies a certain threshold before something is considered a SNP and returns a level of confidence for each identified SNP. The sequences i am using are confirmed variants, i.e. I know that each variation is real. I would like to "simply" create a vcf file containing all differences between the files

• 0 views • link
written by Whetting Whetting

Did you figure out a tool that does this? Also do you mean that any multiple sequence alignments using assembled sequences (assuming the assembly is correct) do not have to go through a "variant calling" approach? What about alignment errors?

• 0 views • link
written by Felix Francis Felix Francis

What format is your data in? We need more information to understand what you are trying to do. What and Why = Best answer.

• 0 views • link
written by Zev.Kronenberg Zev.Kronenberg

This was asked 2.6 years ago. :p

• 0 views • link
written by GouthamAtla GouthamAtla

2 answers

I know this is a very old question, but since I got this as first hit in google today, I will add this answer for anyone who also might need this: There is a pretty nice tool now that does this called "snp-sites"

• 0 views • link
written by Patrick • Patrick

There isn't a program I am aware of that does what you want. However, here are the steps I would take:

  1. Inport your MSA (multiple sequence alignment) into a program that can output variant sites only. Paup* can do this.
  2. Output the matrix and map the gene position to the genomic position.
  3. Write a script that will convert these data to VCF format.
• 172 views • link
written by Zev.Kronenberg Zev.Kronenberg

Log in to answer this question.

More posts like this