Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from sys import argv, exit
- import csv
- import re
- if len(argv) < 3:
- print("Usage: python dna.py data.csv sequence.txt")
- with open(argv[1],"r") as csvfile, open(argv[2],"r") as file: # arguments switched as suggested
- count = 0
- contents = file.read()
- csvcontents = csv.reader(csvfile)
- header = next(csvcontents)
- longestlist = []
- for item in header[1:]:
- beg = 0
- end = len(item)
- seqrun = 0
- longest = 0
- while beg + end <= len(contents):
- seqrun = 0
- while contents[beg: beg + end] == item:
- seqrun += 1
- beg += end
- if seqrun > longest:
- longest = seqrun
- beg += 1
- longestlist.append(str(longest))
- for line in csvcontents:
- Match = False
- if line[1:] == longestlist:
- Match = True
- print(line[0])
- if Match == False:
- print("No match")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement