I would like to align ~15 DNA sequences of 1-2kb, then assess the degree of homology across the alignment using a sliding window. It will be very laborious to do this manually, even using quite large jumps in the windows. Is there a facility to do such analysis automatically with the Clustal Omega tool (or other tool)? Many thanks Catherine Merrick
Too many biologists do indeed say % homology, but that does not mean you should. It is wrong, not just a semantics issue...
Not tested and written on my tiny phone screen, but assuming you can get Biopython running, something like this should generate the input files for your favorite program?
This should print out a bunch of alignments o length window with overlap step.
from Bio import AlignIO window=10 step=5 alignment = AlignIO.read("alignment_in_fasta_format.fas","fasta"): for r in range(0,len(alignment),step): with open("window_"+r+".fas") as out: print >>out, ">%s\n%s" %(alignment[:,r:r+window])