Quick Python script to extract contig summary information (lengths and number of reads, as a tabular file) from an ACE assembly file, using the Biopython ACE parser for convenience.

ace_to_contig_stats.py

#!/usr/bin/env python
#Example usage:
#
#    $ python ace_to_contig_stats.py < example.ace > example_stats.tsv
#
import sys
from Bio.Sequencing import Ace

sys.stdout.write("#Contig\tPadded length\tUnpadded length\tReads\n")
for contig in Ace.parse(sys.stdin):
    seq = contig.sequence
    sys.stdout.write("%s\t%i\t%i\t%i\n" % (contig.name,
                                           len(seq),
                                           len(seq) - seq.count("*"),
                                           contig.nreads))