Don't like ads? PRO users don't see any ads ;-)
Public Pastes
Guest

Untitled

By: a guest on Aug 22nd, 2012  |  syntax: None  |  size: 0.40 KB  |  hits: 7  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. #! /usr/bin/env python
  2.  
  3. from dendropy.interop import genbank
  4. from dendropy.interop import muscle
  5. from dendropy.interop import raxml
  6. gb = genbank.GenBankDna(id_range=(98332, 98359), prefix="AF0")
  7. data = gb.generate_char_matrix(
  8.         label_components=["accession", "organism"])
  9. data = muscle.muscle_align(data)
  10. rr = raxml.RaxmlRunner()
  11. tree = rr.estimate_tree(data, ['-N', '250'])
  12. print tree.as_ascii_plot()
create a new version of this paste RAW Paste Data