Question: Phylogeny based on multiple sequences
0
gravatar for Mehmet
4.3 years ago by
Mehmet540
Japan
Mehmet540 wrote:

Dear All:

I have 113 ortholog groups of 12 species. I want to make a phylogenomics tree based on these orthologs. I was wondering how I could do this. After alignment of each ortholog group separately, should I combine them and then do the phylogenomics tree or do phylogenomics tree of each ortholog groups separately and then combine the each ortholog group`s tree?

Thank you.

alignment next-gen gene genome • 1.3k views
ADD COMMENTlink modified 5 weeks ago by eeb1090 • written 4.3 years ago by Mehmet540
1
gravatar for Nari
4.3 years ago by
Nari880
United States
Nari880 wrote:

You can either align them first and merge alignments to one. or concatenate all representative sequences and align them as one sequence per Ort. group. MUSCLE is my choice for alignment as well as tree generation.

ADD COMMENTlink written 4.3 years ago by Nari880
1

Really fancy way would be to estimate best substitution model for each alignment post gap removal but prior to concatenation and then give that info to RAxML along with the super alignment.

ADD REPLYlink modified 4.3 years ago • written 4.3 years ago by 5heikki8.8k

Hello:

I aligned each ortholog group one by one and merged them by using "cat" command. but I received an error during trimming. The error says sequences are not at the same length. I want to use alignments in nexus format for MrBayes.

ADD REPLYlink written 4.3 years ago by Mehmet540
1

cat will vertically merge them all but you need to merge them horizontally somehow.

ADD REPLYlink written 4.3 years ago by Nari880

I don't know in what format your alignments are but it's almost certain that you can't merge them with cat..

ADD REPLYlink modified 4.3 years ago • written 4.3 years ago by 5heikki8.8k

could you please tell us how to merge alignment files?

ADD REPLYlink written 3.2 years ago by yixinrufeng0

I used a script (FasconCAT), which is available at https://www.zfmk.de/en/research/research-centres-and-groups/fasconcat

ADD REPLYlink written 3.2 years ago by Mehmet540

Instead of cat, paste command is possible to concatenate the sequences side by side. But you have to make sure that the FASTA sequences are in one-liner format before concatenating them using paste. After pasting, remove the delimiter and you will get your concatenated sequences. Remember to pick the delimiter carefully. Pick delimiter that does not exist in your header or any character that might be in your sequences as. Else, it will be a disaster to remove the delimiter later. Eg.

paste -d' ' seq1 seq2 | sed 's/ //g' >> concatenate.fa

The command assigned space as the delimiter and the sed command removed the spaces after pasting. *There are more better way perhaps, for concatenating sequence in command line. This is for sure not the best one, but it somehow work.

ADD REPLYlink modified 5 weeks ago • written 5 weeks ago by eeb1090
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: 1695 users visited in the last hour