The BioPython Bio.AlignIO module will read FASTA output generated using the -m 10 option.

For the fasta parser I wrote a code in python took the idea from Pierre Lindenbaum  in his post  
http://www.biostars.org/post/show/19426/counting-ns-within-fasta/#19439

def fastaio(fh):
        """
        it can take any file handler as input
       eg. fh=open(filename) , sys.stdin  , import gzip; fh=gzip.open(filename)
        """
        buff=[]
        header=[]
        while 1:
                c=fh.read(1)
                if not c or c==">":
                        if len(header)!=0:
                                yield (''.join(header) ,  ''.join(buff))
                        if not c: break
                        header[:]=[] ; buff[:]=[]
                        while 1:
                                hc=fh.read(1)
                                if hc=='\n':break
                                header.append(hc)
                else:
                        buff.append(c)

fh=open("myfile.fasta") for name, seq in fastaio(fh): print ">"+name print seq.replace("\n\s\r","")