From 2f2b3f01a1027d061a7d838263d146bcc5376ce1 Mon Sep 17 00:00:00 2001
From: "Verweij, Stefan" <stefan.verweij@wur.nl>
Date: Fri, 28 Apr 2017 10:21:27 +0200
Subject: [PATCH] Delete 2a_GC.py
---
2a_GC.py | 26 --------------------------
1 file changed, 26 deletions(-)
delete mode 100644 2a_GC.py
diff --git a/2a_GC.py b/2a_GC.py
deleted file mode 100644
index cd6b7e1..0000000
--- a/2a_GC.py
+++ /dev/null
@@ -1,26 +0,0 @@
-#!/usr/bin/python
-
-from Bio.Seq import Seq
-from Bio.Alphabet import IUPAC
-from Bio.SeqUtils import GC
-from Bio import SeqIO
-
-input_file = open('C:\\input_file.Fasta', 'r')
-output_file = open('C:\\output_file.tsv', 'w')
-output_file.write('Gene\tA\tT\tC\tG\tLength\tGC%\n')
-
-for cur_record in SeqIO.parse(input_file, "fasta") :
- gene_name = cur_record.name
- A_count = cur_record.seq.count('A')
- T_count = cur_record.seq.count('T')
- C_count = cur_record.seq.count('C')
- G_count = cur_record.seq.count('G')
- length = len(cur_record.seq)
- GC_content = float(C_count + G_count) / length
- #print GC_content
- output_line = '%s\t%i\t%i\t%i\t%i\t%i\t%f\n' % \
- (gene_name, A_count, T_count, C_count, G_count, length, GC_content)
- output_file.write(output_line)
-
-output_file.close()
-input_file.close()
--
GitLab