Advertisement
joseleeph

Untitled

Dec 14th, 2020
55
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.96 KB | None | 0 0
  1. from sys import argv, exit
  2. import csv
  3. import re
  4. if len(argv) < 3:
  5. print("Usage: python dna.py data.csv sequence.txt")
  6.  
  7. with open(argv[1],"r") as csvfile, open(argv[2],"r") as file: # arguments switched as suggested
  8. count = 0
  9. contents = file.read()
  10. csvcontents = csv.reader(csvfile)
  11.  
  12.  
  13. header = next(csvcontents)
  14. longestlist = []
  15. for item in header[1:]:
  16. beg = 0
  17. end = len(item)
  18. seqrun = 0
  19. longest = 0
  20. while beg + end <= len(contents):
  21. seqrun = 0
  22. while contents[beg: beg + end] == item:
  23. seqrun += 1
  24. beg += end
  25. if seqrun > longest:
  26. longest = seqrun
  27. beg += 1
  28. longestlist.append(str(longest))
  29. for line in csvcontents:
  30. Match = False
  31. if line[1:] == longestlist:
  32. Match = True
  33. print(line[0])
  34. if Match == False:
  35. print("No match")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement