Question: Phylogenetic Tree With Dna, Binary And Multi Data - Raxml Update
1
gravatar for Pavid
8.3 years ago by
Pavid160
Pavid160 wrote:

Hey!

I've made a question about supertrees.
My goal is to integrate data from SNPs (DNA), MIRUs (multi -range around 1 to 25) and Spoligotype (binary) and eventually RFLPs, into a single tree.

@Dan Gaston suggested to use RAxML to obtain the same output.

I've been reading and for I what I understood I need the give a file with all the data (example file that I've created to test) and a partition file.
(I'm testing only with binary and dna data, for now)

In the command line:

./raxmlHPC-PTHREADS-SSE3 -m BINGAMMA -q partition_file.txt -s phylip_file.phy -n Test -T 2

but it gives me a error:

ERROR: Bad base (2) at site 47 of sequence 1
Problem reading alignment file

UPDATE

phylip file with the three types
partition file

command line and error:

./raxmlHPC-PTHREADS-SSE3 -m MULTIGAMMA -q test_part_all.txt -s test_all_types.phy -K ORDERED -n Multi -T 2

Multi State Error, characters must be used in the order they are available, i.e.
0, 1, 2, 3, 4, 5, 6, 7, 8, 9, A, B, C, D, E, F, G, H, I, J, K, L, M, N, O, P, Q, R, S, T, U, V   
You are using the following caharcters: 
0 1 2 3 4 A C G T

This is probably something easy that I'm not seeing.

Thank you guys!

phylogenetics • 5.7k views
ADD COMMENTlink modified 8.3 years ago by Dan Gaston7.1k • written 8.3 years ago by Pavid160
1

well I think so. IN my previous question I spoke about doing a supertree cause I though it was the best solution. I'm working my genotypes from M. tuberculosis. And we have SNP, MIRU, RFLP and spoligotype. So I wanted to integrate these methods into a single tree, as @Gaston suggested I've tried RAxML but I'm doing something wrong

ADD REPLYlink written 8.3 years ago by Pavid160

just as a clarification, is consensifying different trees obtained with different methods into one final tree an option for you?

ADD REPLYlink written 8.3 years ago by 2184687-1231-83-5.0k

Dear Pavid, Could you suggest me any tool for generating partition file. I have a alignment file which constitute 2000 genes and need to find the position of each gene from the alignment file. Is it possible to do by any tool?

ADD REPLYlink written 5 months ago by Dineshkumar K30
1
gravatar for Dan Gaston
8.3 years ago by
Dan Gaston7.1k
Canada
Dan Gaston7.1k wrote:

Your first partition looks like it is multi-state data, not binary. Binary data can only be 0, 1, and ? (maybe - I'm not sure). Anyway you have other state characters in there which is what it is complaining about, the first incorrect character for binary data is the character 2 at position 47 of the first sequence.

ADD COMMENTlink written 8.3 years ago by Dan Gaston7.1k
1

yes, you're right. Its multi data. but for multi data this command should work no? ./raxmlHPC-PTHREADS-SSE3 -m MULTIGAMMA -q test_part.txt -s test_bin_dna.phy -K MK -n Test -T 2 it gives me the same error

ADD REPLYlink written 8.3 years ago by Pavid160
1

ok, I've corrected the positions and integrate them all in one file. now it gives me this error: Multi State Error, characters must be used in the order they are available, i.e. 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, A, B, C, D, E, F, G, H, I, J, K, L, M, N, O, P, Q, R, S, T, U, V You are using the following caharcters: 0 1 2 3 4 A C G T

ADD REPLYlink written 8.3 years ago by Pavid160

do I have to make a script to transform the data?

ADD REPLYlink written 8.3 years ago by Pavid160

Your partitions are marked wrong in the partitions file. Your Multistate data doesn't go all the way to 99. Double check your counts and how you have it set up.

ADD REPLYlink written 8.3 years ago by Dan Gaston7.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1771 users visited in the last hour