RAxML outgroup seemingly not working
0
0
Entering edit mode
5.1 years ago
dthorbur ★ 1.9k

I'm using RAxML to create gene trees. A call I use is simplified below (I've taken out the full paths, bootstraps, and seeds):

raxmlHPC -f d -m GTRGAMMA  -o LG14_Gar_ENSLOCT00000011644_ENSGACT00000016797_CDS -s All_ENSGACG00000015146_chrI_CDS_aln.fa -n All_ENSGACG00000015146_chrI_CDS

The outgroup I have specified as the sequence in the fasta file, but without the > character. The result I get is something like this:

enter image description here

I'm at a loss why there is neither an error message, or why the specified outgroup doesn't appear to be an outgroup. The rest of the sequences are intra-specific but from different populations, so naturally we would expect the inter-specific sequences to be the most different. Yet RAxML doesn't appear to agree with that. The sequences have been aligned with MAFFT using a codon alignment algorithm.

Any help would be greatly appreciated.

RAxML outgroup • 1.7k views
ADD COMMENT

Login before adding your answer.

Traffic: 2405 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6