Added simple script that converts a single-line fasta to a line-wrapped fasta

Hector Galvez · Hector Galvez · commit 4c2b60b609bd · 2019-12-16T12:22:54.000-05:00
diff --git a/python/wrapFasta.py b/python/wrapFasta.py
@@ -0,0 +1,21 @@
+#!/usr/bin/env python
+
+from Bio.Seq import Seq
+from Bio import SeqIO 
+import sys
+
+'''
+Simple script that takes in a fasta file and returns a line-warped fasta. 
+'''
+
+infile = sys.argv[1]
+
+if infile[-3:] == '.fa': 
+    outfile = infile.replace('.fa', '.wrap.fa') 
+elif infile[-6:] == '.fasta': 
+    outfile = infile.replace('.fasta', '.wrap.fasta')
+else: 
+    print("Unsopported filetype")
+
+SeqIO.convert(infile, 'fasta', outfile, 'fasta') 
+
