This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how can I simplfy a newick file with 10 000 entries into a tree of only 10 entries

Hi, wish to pick 10 species from a tree at random, and rewrite the newick file accordingly. I am programming the code to do that.

How can I do it? can someone clarify the logic?

genome newick
• 2,502 views • link
updated by Ram Ram • written by monstercolorfun • monstercolorfun

1 answer

if you just want to sample 10 random tips from a large tree, this would the way to do it using ETE

from ete3 import Tree
import random
tree = Tree("myNewickFile.nw")
sample_tips = random.sample(tree.get_leaves(), 10)
tree.prune(sample_tips)
print tree.write()
• 0 views • link
written by jhc jhc

You can edit the newick tree file using Sed,Awk, grep commands. As long as you are extracting the ten species out with newick defined values it would seem appropriate..

However..

As you are wanting to redo the tree. I really suggest you download the accession sequences of interest ( your ten) , perform an alignment with a program such as MAAFT then use IQtree. In the manual there are nice examples on how to select the best model, and settings. Between MAAFT and iqtree I would maybe also download an alignment program to manually check your alignments for flanking regions with poor coverage and trim them on the consensus length of all ten sequences. It sounds more complicated than it is, but pretty easy.

• 1 views • link
written by Biogeek Biogeek

Log in to answer this question.

More posts like this