I have Multiple sequence alignment (clustal) file and I want to read this file and arrange nucleotide sequences in such a way that it looks more clear and precise in order.
I am doing this from biopython using AlignIO object. My codes goes like this:
alignment = AlignIO.read("opuntia.aln", "clustal") print "Number of rows: %i" % len(align) for record in alignment: print "%s - %s" % record.id, record.seq)
I have been more fascinated by the output of http://www.ebi.ac.uk/Tools/clustalw2/ project because it looks clear.
My Output -- , it looks messy and long scrolling. What i want to do is print only 50 sequences (nucleotide) in each line and continue till the end of alignment file.
I wish to have output like this from http://www.ebi.ac.uk/Tools/clustalw2/.
Any suggestions, algorithm and sample code is appreciated
Here is the link to the input file-- http://pastebin.com/EaeKsyvg
Thanks in advance