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?.
Question: How to extract genomic positions of SNPs using snp-sites?
19 months ago by
pltbiotech_tkarthi • 180
pltbiotech_tkarthi • 180 wrote:
ADD COMMENT • link •
4 weeks ago by
Jorge Amigo ♦ 12k
Santiago de Compostela, Spain
Jorge Amigo ♦ 12k wrote:
There's an optional vcf output format:
You should be able to get the positions from that resulting vcf file.
Please log in to add an answer.
Powered by Biostar version 2.3.0
Traffic: 1959 users visited in the last hour