# HG changeset patch
# User bornea
# Date 1470517497 14400
# Node ID 5990c4dbfbaa2ab45bc7a111d7a9dc6164015333
# Parent  af454e5a9ef5bc73d67b79aa65d1befe2efb6445
Uploaded

diff -r af454e5a9ef5 -r 5990c4dbfbaa filter_fasta.py
--- a/filter_fasta.py	Sat Aug 06 17:00:41 2016 -0400
+++ b/filter_fasta.py	Sat Aug 06 17:04:57 2016 -0400
@@ -22,7 +22,6 @@
             line = line.strip()
             temp = line.split('\t')
             output.append(temp)
-            print temp
     return output
 def getAccessions(infile): # get list of protein accessions from your data
     data = readtab(infile)
@@ -80,10 +79,9 @@
                 if merged!="":
                     seq.append(merged)
     cnt=0
-    with open("output.fasta","wt") as x:
+    with open("output.fasta","w") as x:
         for i in header:
             x.write(i+'\n'+seq[cnt]+'\n')
-            print i+'\n'+seq[cnt]+'\n'
             cnt+=1
 fasta = sys.argv[1] # fasta file to filter
 data = sys.argv[2] # scaffold report #2 -- filename