Advertisement
Guest User

Untitled

a guest
Nov 24th, 2014
137
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.90 KB | None | 0 0
  1. from Bio.SeqUtils import *
  2. nt_search(chromosome, scaffold)
  3.  
  4. from Bio.Seq import Seq
  5. from Bio.SeqUtils import *
  6.  
  7. def findCoordinates (self):
  8. ''' Performs scaffold string search within chromosomes. Returns scaffold coordinates within said chromosome. '''
  9.  
  10. for chromosome in self.chromosomes.keys():
  11. for scaffold in self.scaffolds.keys():
  12. # search forward strand.
  13. nt_forward = nt_search(self.chromosomes[chromosome], self.scaffolds[scaffold])
  14. if len(nt_forward) > 1:
  15. startCoord = nt_forward[1] + 1
  16. endCoord = (startCoord + len(self.scaffolds[scaffold]))
  17. # save coordinates
  18.  
  19. else: # search reverse strand
  20. scaffold_seq = Seq(self.scaffolds[scaffold])
  21. reverse_seq = scaffold_seq.reverse_complement()
  22. nt_reverse = nt_search(self.chromosomes[chromosome], str(reverse_seq))
  23. if len(nt_reverse) > 1:
  24. startCoord = nt_reverse[1] + 1
  25. endCoord = (startCoord + len(self.scaffolds[scaffold]))
  26. # save coordinates
  27. self.scaffolds[scaffold] = str(scaffold_seq.reverse_complement())
  28.  
  29. Traceback (most recent call last):
  30. File "scaffoldPlacer.py", line 98, in <module>
  31. z.findCoordinates()
  32. File "scaffoldPlacer.py", line 60, in findCoordinates
  33. nt_reverse = nt_search(self.chromosomes[chromosome], str(reverse_seq))
  34. File "/usr/local/lib/python2.7/site-packages/Bio/SeqUtils/__init__.py", line 191, in nt_search
  35. m = re.search(pattern, s)
  36. File "/usr/local/lib/python2.7/re.py", line 142, in search
  37. return _compile(pattern, flags).search(string)
  38. File "/usr/local/lib/python2.7/re.py", line 243, in _compile
  39. p = sre_compile.compile(pattern, flags)
  40. File "/usr/local/lib/python2.7/sre_compile.py", line 523, in compile
  41. groupindex, indexgroup
  42. OverflowError: regular expression code size limit exceeded
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement