How to calculate pairwise ts/tv ratio given a alignment

0

Entering edit mode

9.6 years ago

YMTIO • 0

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 • 2.9k views

ADD COMMENT • link updated 2.4 years ago by Ram 44k • written 9.6 years ago by YMTIO • 0

1

Entering edit mode

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

ADD REPLY • link 9.6 years ago by Dan D 7.4k

0

Entering edit mode

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

ADD REPLY • link updated 2.4 years ago by Ram 44k • written 9.6 years ago by YMTIO • 0

0

Entering edit mode

Show us a sample please.

ADD REPLY • link updated 2.4 years ago by Ram 44k • written 9.6 years ago by Pierre Lindenbaum 164k

0

Entering edit mode

The alignment (fasta) is like:

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

ADD REPLY • link updated 2.4 years ago by Ram 44k • written 9.6 years ago by YMTIO • 0

Login before adding your answer.