[Biopython] suggestion for a little change in the ACE cookbook
    Fungazid 
    fungazid at yahoo.com
       
    Sun Jul  5 22:57:08 UTC 2009
    
    
  
Hi,
About the cookbook here
http://biopython.org/wiki/ACE_contig_to_alignment
instead of:
def cut_ends(read, start, end):
  return (start-1) * '-' + read[start-1:end] + (end +1) * '-'
I think it is better to write:
def cut_ends(self,read, start, end):
    return (start-1) * 'x' + read[start-1:end-1] + (len(read)-end) * 'x'
The 2 changes are:
1) correcting the coordinates of the clipped 5' region
2) adding 'x' instead of '-' to separate the clipped region from the gaps
      
    
    
More information about the Biopython
mailing list