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?.
•
written
by
pltbiotech_tkarthi
2 answers
Hey, wonder if you find a solution? I'm also looking for this..
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.
Log in to answer this question.
-
Samtools mpileup for RNA Editing levels generating empty output file
written by haydenshinn •
-
Incorporate SNP/Indel and generate new sequence from original reference gene sequence
written by amar.kant20422 •
-
Extracting amino acid positions from a fasta sequence
written by akang
-
prepare a GFF file for MOCK fasta reference
written by Sam
-
Simple way to extract or mask specific sites only from multi-FASTA alignments
written by Kame •
-
VariantCalling with single sanger read
written by Amirosein
-
A snp caller from a multiple alignment using python/biopython
written by cnyaigoti •
-
Comparing Snps Called Using Two Versions Of A Genome
written by rob234king
-
Determine The Position Of A Snp In A Multiple Sequence Alignment
written by Whetting
-
Mapping Genomic Position To Protein Position
written by John Knowlg •
