This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to calculate pairwise ts/tv ratio given a alignment

Hi Everyone

Could anyone tell me how to calculate pairwise ts/tv ratio given a alignment (the 4-fold degenerate sites extracted from a CDS alignment)? which software or script can do this?

One example for the alignment (fasta) is like:

>taxa001
TATATATTTTTTTTAAAAATAATAATATTA
>taxa002
TATATATTCTTACTAAAAATAATAATATTA
>taxa003
TATATATTCTTACTAAAAATAATAATATTA
>taxa004
TATATATTCTTACTAAAAATAATAATATTA
>taxa005
TATATATTTTTTTTAAAAATAATAATATTA
>taxa006
TATATATTCTTACTAAAAATAATAATATTA
substitution software alignment

Can you define what you mean by "alignments?" Typically for NGS data you're going to pull ts/tv from the variant (VCF) data.

The alignment is the 4-fold degenerate sites extracted from a CDS alignment

The alignment (fasta) is like:

>taxa001
TATATATTTTTTTTAAAAATAATAATATTA
>taxa002
TATATATTCTTACTAAAAATAATAATATTA
>taxa003
TATATATTCTTACTAAAAATAATAATATTA
>taxa004
TATATATTCTTACTAAAAATAATAATATTA
>taxa005
TATATATTTTTTTTAAAAATAATAATATTA
>taxa006
TATATATTCTTACTAAAAATAATAATATTA

0 answers

No answers yet.

Log in to answer this question.