How to get the number of segregating sites from a multiple alignment?
Hi everyone!
I'm struggling to get the number of segregating sites (S) in a large contig in a fasta format (1.3Mb) with multiple sequences. I've used muscle aligner to get this alignment, and now I need the number, and % of segregating sites across all sequences from my contig. I didn't find info if muscle can provide it.
Any tips? Thanks!
• 899 views
•
link
updated
by
Ram
•
written
by
Flávia Rogério •
•
written
by
Flávia Rogério •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Which alignment for whole bacterial genomes for phylogenetic tree?
written by Dom •
-
Summarise bases from alignment
written by mpm20 •
-
Weird behavior of `--max-depth` argument from `bcftools mpileup` command when its value is too low
written by sbstevenlee •
-
Per position amino acid frequency from mapping short reads vs. long coding amplicon
written by jmartin77777 •
-
DNA input for HMMER
written by bn •
-
[ Mustang ] Problem with Server
written by ReWeeda
-
Obtaining contig stats from Velvet+Oases De Novo Assembly
written by blazer9131 •
-
samtools/bcftools strange (error ?) lack of call
written by mbourgey •
-
python subprocesses and wrappers for Jce tool
written by mossig89 •
-
Retrieving Sequences using NCBI Gene database IDs
written by gavingray1729 •
