Hi,
Â
I am new to R and I am
trying to read a file containing sets of multiple sequence alignments as in
the example below. The difference is that for my files, each set of sequences
may contain up to 50 sequences of 10,000bp each.
Â
Example:
 5 10
Seq0Â Â Â Â ATCGTTATTA
Seq1Â Â Â Â
Hi Gwennaël,
Thanks very much!
I did look at what might be reading which I get the following:
50 DNA sequences in binary format stored in a matrix.
All sequences of same length: 1
Labels: Seq0 Seq49 Seq48 Seq47 Seq46 Seq45 ...
Base composition:
   a    c    g    t
Does anyone know of an imlplementation of the RelTime algorithm (
1. Tamura K,
2. Battistuzzi FU,
3. Billing-Ross P,
4. Murillo O,
5. Filipski A,
6. Kumar S
. Estimating divergence times in large molecular phylogenies. Proc Natl
Acad Sci U S A.2012;109:19333-19338.) in R?
Apologies for not replying to the list correctly the first time!
Hi Fabricia,
This may look a little complicated but it's pretty easy to get any string of
characters into a matrix using the s2c function in seqinr and from there it's
simple to get them into ape's DNAbin by running the matrix