Mercurial > repos > jhliu > fastq_c2t

changeset 0:1dd7e185e9b9 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author jhliu
date Thu, 03 Jun 2021 05:09:43 +0000
parents
children
files fastq_c2t.py
diffstat 1 files changed, 30 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/fastq_c2t.py	Thu Jun 03 05:09:43 2021 +0000
@@ -0,0 +1,30 @@
+#!/usr/bin/env python
+import sys
+import os
+import getopt
+from Bio import SeqIO
+
+opts,args=getopt.getopt(sys.argv[1:],"i:h")
+for op,value in opts:
+	if op == "-i":
+		file = value
+	elif op == "-h":
+		sys.exit()
+
+with open(file,'r') as input:
+	line = input.readline()
+	i = 0
+	while (line):
+		if i == 0:
+			print line.strip("\n")
+			i += 1
+		elif i == 1:
+			print line.strip("\n").replace("C","T")
+			i += 1
+		elif i == 2:
+			print line.strip("\n")
+			i +=1
+		elif i == 3:
+			print line.strip("\n")
+			i = 0
+		line = input.readline()
\ No newline at end of file