How to extract genomic positions of SNPs using snp-sites?

Question

[pltbiotech_tkarthi]

By using snp-sites (https://github.com/sanger-pathogens/snp-sites) I extracted the SNPs with Linux command: snp-sites file name.fasta. This script just extracting the SNPs without position information. I know the genomic positions of my reference sequence. My question is: How to add the genomic position of reference sequence to the aligned fasta sequences as an input to obtain the result of SNPs with genomic positions information as an output. Please let me know an additional script to extract genomic positions of SNPs output from the aligned fasta sequences, besides the SNPs extracting script (snp-sites file name.fasta). Related post: C: How to extract SNPs from multiple alignment fasta file?.

Answer 1

[christine]

Hey, wonder if you find a solution? I'm also looking for this..

Answer 2

[Jorge Amigo]

There's an optional vcf output format:

https://github.com/sanger-pathogens/snp-sites#output

You should be able to get the positions from that resulting vcf file.