Mercurial > repos > bornea > filter_fasta

changeset 6:5990c4dbfbaa draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author bornea
date Sat, 06 Aug 2016 17:04:57 -0400
parents af454e5a9ef5
children e01b49c112b0
files filter_fasta.py
diffstat 1 files changed, 1 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- 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