This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SNP calling with fasta files

Hi,

I would like to generate a vcf file listing SNPs present between my bacterial species of interest and the reference genome. However I only have access to FASTA files available on NCBI.

Was wondering if anyone knows whether it's possible to call SNPs with assembled genome sequences?

Thanks, Joyce

snp
• 10,888 views • link
updated by Ram Ram • written by joycewang83 • joycewang83

2 answers

if you want to call SNPs you would need several overlapping sequences (such as NGS reads, or different EST sequences) and the quality of the bases in order to let an algorithm infer whether a base mismatch is a SNP or not.

but if you just want to have the differences between your reference sequence and any FASTA file you may download then I would suggest you to pair align each FASTA sequence with the reference sequence and then look for information on how to get an VCF file from a FASTA alignment, such as in Getting A Vcf File From A Fasta Alignment

• 0 views • link
written by Jorge Amigo Jorge Amigo

If you are dealing with two assembled genomes, you could use a genome aligner such as MAUVE to identify the divergent regions(including SNPs) between them. You can write a simple script to convert this data to VCF format.

• 0 views • link
updated by Ram Ram • written by Felix Francis Felix Francis

Thank you Felix, I have used MAUVE before but am completely inexperienced with script writing!

Jorge, thank you for your explanation and directing me to the other forum. The person who started the forum actually recommended a GUI platform for generating vcf files and I think it might work! Thanks again :)

Joyce

• 0 views • link
updated by Ram Ram • written by joycewang83 • joycewang83

Log in to answer this question.

More posts like this