Added script for metabinner output conversion

bin/create_metabinner_bins.py

@@ -0,0 +1,36 @@
+#!/usr/bin/env python
+
+import sys
+import os
+from Bio import SeqIO
+
+binning = sys.argv[1]
+fasta = sys.argv[2]
+path = sys.argv[3]
+prefix = sys.argv[4]
+
+try:
+    os.mkdir(path)
+except OSerror as error:
+    pass
+
+bins = {}
+with open(binning, 'r') as b:
+    for line in b:
+        contig, bin = line.strip().split('\t')
+        bins[contig] = bin
+
+for outf in set(bins.values()):
+    name = prefix + outf + ".fa"
+    open(os.path.join(path, name), 'w')
+
+with open(fasta) as handle:
+    for record in SeqIO.parse(handle, "fasta"):
+        if len(record) < 1000:
+            f = "tooShort.fa"
+        elif record.id not in bins:
+            f = "unbinned.fa"
+        else:
+            f = prefix + bins[record.id] + ".fa"
+        with open(os.path.join(path, f), 'a') as out:
+            SeqIO.write(record, out, "fasta")