def dump_cluster(c):
if os.path.exists("%s_dist_%s_aligned_short.fasta-cluster%d.png" % (prefix, dist, c)):
print """
## Subcluster %d
""" % (c)
if c == 0:
print "Cluster 0 represents isolates that do not cluster with any other isolates within the distance cut-off, i.e. singleton sequences. The sequences presented are unrelated."
print """

""" % (c, prefix, dist, c)
else:
print """
## Subcluster %s
(Tree not shown for clusters with <5 isolates)
Isolates:
""" % (c)
for rec in SeqIO.parse(open("%s_dist_%s_aligned_short.fasta-cluster%d" % (prefix, dist, c)), "fasta"):
print " - %s" % (rec.id)
评论列表
文章目录