A snp caller from a multiple alignment using python/biopython
Hi,
I want to call all nucleotide SNP positions from a virus genome multiple sequence alignment of about 15000 nucleotides long. Any one done this? I will prefer to use python I will like the script to return an output of the number of SNP positions detected, actual position in the original alignment and the alignment column the SNP was detected. Gaps to be ignored.
I will appreciate,
Charles
• 2,723 views
•
link
updated
by
Ram
•
written
by
cnyaigoti •
•
written
by
cnyaigoti •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
sequenced noise elimination in searching for snp
written by gagi1993 •
-
Plot nucleotide diversity with multiple sequence alignment
written by asoundi96 •
-
Does Nucleotide Alignment can Translate into Proteins (Amino Acids)?
written by pfee418 •
-
Delete certain nucleotide positions in multiple nucleotide alignment
written by mariasv1 •
-
How do you find the number of SNPs per gene given a genotype and a map file?
written by eternalshadowrose •
-
SNPeff: can I just use hg19?
written by Marvin
-
how to find SNPs within a range of given windows?
written by Ana
-
Obtain heterozyous/IUPAC consensus sequences from multiple sequence alignment
written by Louis Kok •
-
SNP calling,consensus sequence and running Blastn
written by low.weiquan •
-
NJ phylogenetic tree calculation method weighing gaps in the alignment
written by Kame •
Why not just use a premade tool? A quick google for "virus ngs variant call" will yield some useful results.