# HG changeset patch # User bornea # Date 1470517241 14400 # Node ID af454e5a9ef5bc73d67b79aa65d1befe2efb6445 # Parent 489b560129c8192b7dfef403ceb7b5fe519cb6cf Uploaded diff -r 489b560129c8 -r af454e5a9ef5 filter_fasta.py --- a/filter_fasta.py Sat Aug 06 16:57:14 2016 -0400 +++ b/filter_fasta.py Sat Aug 06 17:00:41 2016 -0400 @@ -1,4 +1,4 @@ -# -*- coding: utf-8 -*- + """ Python-code: Merge Scaffold Samples Report files @author = Brent Kuenzi @@ -22,6 +22,7 @@ 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) @@ -82,6 +83,7 @@ with open("output.fasta","wt") 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