I have used multiple sequences aligned file of protein data to generate a maximum parsimony tree, then I used the "bootstrap_trees(msa, times, tree_constructor)" function from biopython to generate bootstrap replicates of my MSA data. Since I have my original tree and bootstrap replicates for my tree, I want to ask how do I compare the topology (using python/biopython) for my tree with the bootstrap replicates so that for each edge in original tree the bootstrap support is shown. (Bootstrap support is the number of bootstrap trees that edge appears in.) Any help is greatly appreciated.

  • $\begingroup$ majority_consensus or similar? That said, MP is not a good method for protein trees. You should use a more sophisticated tree builder. $\endgroup$ – user172818 Dec 25 '18 at 3:48
  • $\begingroup$ Actually, I have tried the generating the bootstrap replicates for my msa data. Now I need to count the number of times the same branching pattern appears in the bootstap replicates. How can I do this? $\endgroup$ – Sidra Younas Dec 28 '18 at 14:11

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.