Skip to content

Commit 3653b5b

Browse files
authored
Merge pull request #24 from dpellow/master
Bug fix for SPAdes latest version reverse complement node ordering
2 parents b792015 + 535ec99 commit 3653b5b

File tree

1 file changed

+5
-1
lines changed

1 file changed

+5
-1
lines changed

bin/make_fasta_from_fastg.py

Lines changed: 5 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -20,9 +20,13 @@ def parse_user_input():
2020

2121
def parse_lines(fastg, ofile):
2222
fp = open(fastg, 'r')
23+
seen = set() ##
2324
for name,seq,qual in readfq(fp):
2425
name = re.sub('[:,]'," ", name[:-1]).split(" ")[0]
25-
if name[-1] == "'": continue # only need one of forward and reverse complement
26+
check_seen = name
27+
if check_seen[-1] == "'": check_seen = check_seen[:-1]
28+
if check_seen in seen: continue
29+
else: seen.add(check_seen)
2630
line = ">"+name+"\n"+seq+"\n"
2731
ofile.write(line)
2832

0 commit comments

Comments
 (0)