Advertisement
Guest User

Untitled

a guest
Mar 28th, 2015
234
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.35 KB | None | 0 0
  1. def same_seq(seq1, seq2):
  2. # Set placeholder values for larger and smaller sequences
  3. smaller = ""
  4. larger = ""
  5. # Find the larger sequence in the pair
  6. if (len(seq1) > len(seq2)):
  7. larger = seq1
  8. smaller = seq2
  9. else:
  10. larger = seq2
  11. smaller = seq1
  12.  
  13. pair = ""
  14. mypairs = []
  15. for i in range (len(smaller)):
  16. temp = ""
  17. for j in range (i, len(smaller)):
  18. temp = temp + smaller [j]
  19. if (temp in larger):
  20. if (len(temp) > len(pair)):
  21. pair = temp
  22. mypairs = []
  23. mypairs.append(temp)
  24. elif (len(temp) == len(pair) and not temp in mypairs):
  25. mypairs.append(temp)
  26. # mypairs [] contains the largest pairs if any exist
  27. return (mypairs)
  28.  
  29.  
  30.  
  31.  
  32. def get_next_pair (file):
  33. # Read the file line by line
  34. p1 = file.readline()
  35. p2 = file.readline()
  36. return p1, p2
  37.  
  38. def main():
  39. pairs_list = []
  40. # Open the external file
  41. f = open ('dna.txt' , 'r')
  42. num_pairs = int (f.readline())
  43. # Set a placeholding value for each pair in a sequence
  44. p1 = ""
  45. p2 = ""
  46. print ("Longest Common Sequences")
  47. print (" ")
  48. temp = ""
  49. for i in range (num_pairs):
  50. temp = print ("Pair",str(i+1)+":")
  51. p1, p2 = get_next_pair (f)
  52. pairs_list = same_seq (p1, p2)
  53. if (len(pairs_list) == 0):
  54. print ("No Common Sequence Found")
  55. else:
  56. for j in range (len(pairs_list)):
  57. print (pairs_list[j])
  58. print ()
  59.  
  60.  
  61. main()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement