Fixed naming scheme for 3rd column

loosolab · Jan 9, 2019 · f6ddce3 · f6ddce3
1 parent b97841f
commit f6ddce3
Show file tree

Hide file tree

Showing 2 changed files with 16 additions and 14 deletions.
diff --git a/.gitignore b/.gitignore
@@ -203,6 +203,4 @@ venv.bak/
 
 # mypy
 .mypy_cache/
-
-Würde bin/3.1_create_gtf/data/ löschen
-Würde data/ löschen
+/bin/3.1_create_gtf/data/
diff --git a/bin/3.1_create_gtf/Modules/ucsc/ucsc.py b/bin/3.1_create_gtf/Modules/ucsc/ucsc.py
@@ -58,17 +58,21 @@ def read_gff_to_gtf(self):
         with open(self.output, 'r') as csvfile:
             tsvreader = csv.reader(csvfile, delimiter='\t')
             for row in tsvreader:
-                sequence = []
-                sequence.append(row[0])
-                sequence.append("UCSC")
-                sequence.append(row[3].lower().replace(' ', '_'))
-                sequence.append(row[1])
-                sequence.append(row[2])
-                sequence.append(".")
-                sequence.append(row[5])
-                sequence.append(".")
-                sequence.append('; '.join([self.find_ID(''.join(row[11:])), 'activity \"'+", ".join(self.get_activity(''.join(row[11:]))) + '"'])+";")
-                gtf_lines.append(sequence)
+                if row[9] not in ["region", "sequence_feature",
+                                  "CAAT_signal", "stem_loop",
+                                  "sequence_secondary_structure"]:
+
+                    sequence = []
+                    sequence.append(row[0])
+                    sequence.append("UCSC")
+                    sequence.append(row[9].lower().replace(' ', '_'))
+                    sequence.append(row[1])
+                    sequence.append(row[2])
+                    sequence.append(".")
+                    sequence.append(row[5])
+                    sequence.append(".")
+                    sequence.append('; '.join([self.find_ID(''.join(row[11:])), 'activity \"'+", ".join(self.get_activity(''.join(row[11:]))) + '"'])+";")
+                    gtf_lines.append(sequence)
 
             return gtf_lines