Skip to content
Snippets Groups Projects
Commit 0beb3cbb authored by lotte.pronk's avatar lotte.pronk
Browse files

Removed a few lines of codes I used for testing, which caused fasta output to...

Removed a few lines of codes I used for testing, which caused fasta output to contain max 11 sequences.
parent 93337acd
No related branches found
No related tags found
No related merge requests found
......@@ -48,13 +48,8 @@ def split_fasta_taxonomy(fasta_file, outdir):
prok_headers_list = prok_headers.split("\n") # list of all prokaryote headers
prok_file.close()
counter = 0
with open(fasta_file, "r") as f:
for title, seq in SimpleFastaParser(f):
if counter > 10:
print("breaking")
break
counter += 1
title_id = title.split(" ")[0]
if title_id in euk_headers_list:
with open(os.path.join(outdir, "eukaryotes.fasta"), "a") as euk:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment