Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ete evol: Inconsistent marking between --mark and --internal options #747

Open
lab83bio opened this issue May 8, 2024 · 0 comments
Open

Comments

@lab83bio
Copy link

lab83bio commented May 8, 2024

Hi,

I am using ete3 evol (conda inst. V. 3.1.2) to analyse pseudogenes in primate genomes. According to the documentation if a tree is marked with --mark and two species separated by three commas

the tree is marked from the common ancestor of the surrounding species

Consistently, using the tree ECP_EDN_15.nw of the documentation page with --mark Chimp_EDN,,,Human_EDN, the tree is marked as follows:

$ ete3 evol -t ECP_EDN_15.nw --alg ECP_EDN_15.fasta --model b_free --mark Chimp_EDN,,,Human_EDN -o outdir/ 
...
       marking branches 3, 7, 21

          (((Hylobates_EDN,(Orang_EDN,(Gorilla_EDN,(Chimp_EDN #1,Human_EDN #1) #1))),(Macaq_EDN,(Cercopith_EDN,(Macaq2_EDN,Papio_EDN)))),(Orang_ECP,((Macaq_ECP,Macaq2_ECP),(Goril_ECP,Chimp_ECP,Human_ECP))));

One should expect that this marked tree is also present among the trees generated by the --internal option. However, by using this option, the program generate a tree marked at descendents but not at the common ancestor.

$ ete3 evol -t ECP_EDN_15.nw --alg ECP_EDN_15.fasta --model b_free --internal -o outdir2/
...
       marking branches 3, 7

          (((Hylobates_EDN,(Orang_EDN,(Gorilla_EDN,(Chimp_EDN #1,Human_EDN #1)))),(Macaq_EDN,(Cercopith_EDN,(Macaq2_EDN,Papio_EDN)))),(Orang_ECP,((Macaq_ECP,Macaq2_ECP),(Goril_ECP,Chimp_ECP,Human_ECP))));

I think that the marking by --mark is correct, at variance with the --internal one, because if descendent branches (3, 7) are marked, also the ancestral branch (21) should be marked.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant