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

Warning for seqs > 32 mbp

parent eac34e3d
Branches
No related tags found
No related merge requests found
...@@ -20,8 +20,8 @@ def size_filter(contig_file, outdir, size): ...@@ -20,8 +20,8 @@ def size_filter(contig_file, outdir, size):
if len(seq) < size: if len(seq) < size:
removed_contigs += 1 removed_contigs += 1
if len(seq) > 32000000: if len(seq) > 32000000:
print("Warning! sequence" + title + "is longer than 32000000 BP. If you did not provide a GFF\ print("Warning! sequence" + title + "is longer than 32000000 BP. If you did not provide a GFF"
file, Prodigal will be used to predict genes. Prodigal does not accept sequences > 32 mbp!") "file, Prodigal will be used to predict genes. Prodigal does not accept sequences > 32 mbp!")
print("Total contigs checked:\t" + str(total_contigs), print("Total contigs checked:\t" + str(total_contigs),
'\nNumber of contigs >=' + str(size) + 'bp:\t' + str(kept_contigs), '\nNumber of contigs >=' + str(size) + 'bp:\t' + str(kept_contigs),
...@@ -48,4 +48,4 @@ def split_fasta_taxonomy(fasta_file, outdir, headerfile): ...@@ -48,4 +48,4 @@ def split_fasta_taxonomy(fasta_file, outdir, headerfile):
g.write(">" + title + "\n" + seq + "\n") g.write(">" + title + "\n" + seq + "\n")
else: else:
with open(os.path.join(outdir, "prokaryotes.fasta"), "a") as g: with open(os.path.join(outdir, "prokaryotes.fasta"), "a") as g:
g.write(">" + title + "\n" + seq + "\n") g.write(">" + title + "\n" + seq + "\n")
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment