# HG changeset patch # User jhliu # Date 1622696983 0 # Node ID 1dd7e185e9b94ed93630eaa9ed8aafdbed51d9b9 Uploaded diff -r 000000000000 -r 1dd7e185e9b9 fastq_c2t.py --- /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