<p>I think you have a bug</p>

<pre><code>raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup -p 10000 -s /Users)/username/Desktop/alignment_file.phy
</code></pre>

<p>Try the above, where username is your username. I'll check my codes later. The most important bit is the space between the -p and word outgroup. You can omit the <code>-p 10000</code> it is only needed  for very geeky maximum likelihood and sets the random number stream. </p>

CC BY-SA 4.0
2019-10-07T01:10:54.413Z
Thank you @Michael G. I tried using this line:  raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup -p 10000 -s /Users/username/Desktop/alignment_file.phy but I get the following error: The file /Users/username/Desktop/alignment_file.phy you want to open for reading does not exist, exiting ...  I also tried using raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup -p 10000 -s /username/Desktop/alignment_file.phy but I got the same error message. As mentioned above, the file name is simply alignment.phy.
Leah
CC BY-SA 4.0
2019-10-07T09:35:20.837Z
I've set up a chat group here, https://chat.stackexchange.com/rooms/99546/phylogeny. Please drop by and let me know when you are around. Tomorrow afternoon your time zone would be good. Please make sure you have access to your Linux/Darwin terminal and your data sets
M__
CC BY-SA 4.0
2019-10-07T18:18:00.463Z
Thanks a lot @Michael G. I will be there on Tuesday October 8 from 17:00 (CDT) onward. Does that work for you? I can connect earlier if it suits you better.
Leah
CC BY-SA 4.0
2019-10-07T18:40:25.153Z
Yes earlier, ie before 17:00 please (6 hour time difference). Must leave by 17:00
M__
CC BY-SA 4.0
2019-10-07T20:39:13.293Z
I see. How about 9:30 am CDT (UTC -5), which would be 15:30 BST (UTC +1)? Is this better?
Leah
CC BY-SA 4.0
2019-10-07T21:20:50.473Z
Ok, see you then
M__
CC BY-SA 4.0
2019-10-05T09:53:57.450Z
2019-10-05T16:30:20.227Z
2019-10-05T16:30:20.227Z
M__
<p>To clear up what may be a possible misunderstanding about RAxML: RAxML (and most/all other maximum likelihood phylogenetic inference programs using <strong>reversible</strong> models of sequence evolution) explicitly infer unrooted trees. This is because, for a reversible model of sequence evolution, changing the position of the root does not affect the likelihood of the tree. Therefore, as you will see in the RAxML manual, specifying an outgroup is explicitly a display option - it doesn't affect the inference in any way. This means that it's usually easier to just infer the tree in RAxML without specifying any outgroup, and then use any tree viewing software (e.g. <a href="http://tree.bio.ed.ac.uk/software/figtree/" rel="nofollow noreferrer">FigTree</a>, <a href="http://phylotree.hyphy.org/" rel="nofollow noreferrer">phylotree.js</a>, or equivalent) to reroot the tree as you need to.  </p>

CC BY-SA 4.0
2019-10-09T09:54:18.990Z
2019-10-09T09:54:18.990Z
NatWH
<p>I building a robust amino acid phylogeny with adequate robustness support (see previous post <a href="https://bioinformatics.stackexchange.com/questions/10365/how-to-make-a-robust-amino-acid-phylogeny-with-adequate-robustness-support">here</a>).</p>
<p>This is a brief summary of what I have done:</p>
<ol>
<li>I performed a BLASTp analysis using a protein sequence from species <em><strong>A</strong></em> as query against the full NCBI protein database,</li>
<li>downloaded all the hits in fasta format and I added the sequences I identified with the query sequence in my <em><strong>target</strong></em> genome, which is closely related to species <em><strong>A</strong></em>.</li>
<li>I then aligned the sequences using the Compute - Alignment - MUSCLE option in Seqotron and I saved my alignment file in Phylip format.</li>
<li>I tried to run this command line in RAxML to build a robust maximum likelihood tree:</li>
</ol>
<p><code>raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup-p 10000 -s ~/Desktop/alignment_file.phy</code></p>
<p>However, the program showed the following error:</p>
<p><em><strong>Error, the outgroup name &quot;outgroup-p&quot; you specified can not be found in the alignment, exiting ....</strong></em></p>
<p><strong>Question</strong> how can I build an outgroup for my phylogenetic tree, and how do I include it in my alignment?</p>
<p>I appreciate your help and suggestions!</p>

CC BY-SA 4.0
2019-10-05T09:47:32.183Z
Can you show us the input file please? Personally I would avoid using a '-' in the filename
M__
CC BY-SA 4.0
2019-10-07T01:08:57.910Z
Sure. The input file is named simply "alignment.phy" to avoid any problems.
Leah
<phylogenetics><phylogeny><duplications>
How to define an outgroup to build a robust amino acid tree
# How to define an outgroup to build a robust amino acid tree

(Asked by: username_0 on 2019-10-05)

I building a robust amino acid phylogeny with adequate robustness support (see previous post [here](https://bioinformatics.stackexchange.com/questions/10365/how-to-make-a-robust-amino-acid-phylogeny-with-adequate-robustness-support)).

This is a brief summary of what I have done:

1. I performed a BLASTp analysis using a protein sequence from species ***A*** as query against the full NCBI protein database,
2. downloaded all the hits in fasta format and I added the sequences I identified with the query sequence in my ***target*** genome, which is closely related to species ***A***.
3. I then aligned the sequences using the Compute - Alignment - MUSCLE option in Seqotron and I saved my alignment file in Phylip format.
4. I tried to run this command line in RAxML to build a robust maximum likelihood tree:

`raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup-p 10000 -s ~/Desktop/alignment_file.phy`

However, the program showed the following error:

***Error, the outgroup name "outgroup-p" you specified can not be found in the alignment, exiting ....***

**Question** how can I build an outgroup for my phylogenetic tree, and how do I include it in my alignment?

I appreciate your help and suggestions!

**Tags:** phylogenetics, phylogeny, duplications
> **Comments:**
> **username_1**: Can you show us the input file please? Personally I would avoid using a '-' in the filename  
> **username_0**: Sure. The input file is named simply "alignment.phy" to avoid any problems.  

---

## 2 answers

### Answer by username_1 on 2019-10-05. Score: 2 | Accepted answer

I think you have a bug
raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup -p 10000 -s /Users)/username/Desktop/alignment_file.phy

Try the above, where username is your username. I'll check my codes later. The most important bit is the space between the -p and word outgroup. You can omit the -p 10000 it is only needed  for very geeky maximum likelihood and sets the random number stream.

> **Comments:**
> **username_0**: Thank you @Michael G. I tried using this line:  raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup -p 10000 -s /Users/username/Desktop/alignment_file.phy but I get the following error: The file /Users/username/Desktop/alignment_file.phy you want to open for reading does not exist, exiting ...  I also tried using raxmlHPC -m PROTGAMMAILG -n output.tre -o outgroup -p 10000 -s /username/Desktop/alignment_file.phy but I got the same error message. As mentioned above, the file name is simply alignment.phy.  
> **username_1**: I've set up a chat group here, https://chat.stackexchange.com/rooms/99546/phylogeny. Please drop by and let me know when you are around. Tomorrow afternoon your time zone would be good. Please make sure you have access to your Linux/Darwin terminal and your data sets  
> **username_0**: Thanks a lot @Michael G. I will be there on Tuesday October 8 from 17:00 (CDT) onward. Does that work for you? I can connect earlier if it suits you better.  
> **username_1**: Yes earlier, ie before 17:00 please (6 hour time difference). Must leave by 17:00  
> **username_0**: I see. How about 9:30 am CDT (UTC -5), which would be 15:30 BST (UTC +1)? Is this better?  
> **username_1**: Ok, see you then  

---

### Answer by username_2 on 2019-10-09. Score: 2

To clear up what may be a possible misunderstanding about RAxML: RAxML (and most/all other maximum likelihood phylogenetic inference programs using **reversible** models of sequence evolution) explicitly infer unrooted trees. This is because, for a reversible model of sequence evolution, changing the position of the root does not affect the likelihood of the tree. Therefore, as you will see in the RAxML manual, specifying an outgroup is explicitly a display option - it doesn't affect the inference in any way. This means that it's usually easier to just infer the tree in RAxML without specifying any outgroup, and then use any tree viewing software (e.g. [FigTree](http://tree.bio.ed.ac.uk/software/figtree/), [phylotree.js](http://phylotree.hyphy.org/), or equivalent) to reroot the tree as you need to.


---

