annotate fastq_c2t.py @ 0:1dd7e185e9b9 draft default tip

Uploaded
author jhliu
date Thu, 03 Jun 2021 05:09:43 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
1 #!/usr/bin/env python
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
2 import sys
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
3 import os
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
4 import getopt
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
5 from Bio import SeqIO
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
6
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
7 opts,args=getopt.getopt(sys.argv[1:],"i:h")
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
8 for op,value in opts:
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
9 if op == "-i":
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
10 file = value
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
11 elif op == "-h":
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
12 sys.exit()
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
13
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
14 with open(file,'r') as input:
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
15 line = input.readline()
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
16 i = 0
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
17 while (line):
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
18 if i == 0:
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
19 print line.strip("\n")
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
20 i += 1
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
21 elif i == 1:
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
22 print line.strip("\n").replace("C","T")
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
23 i += 1
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
24 elif i == 2:
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
25 print line.strip("\n")
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
26 i +=1
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
27 elif i == 3:
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
28 print line.strip("\n")
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
29 i = 0
1dd7e185e9b9 Uploaded
jhliu
parents:
diff changeset
30 line = input.readline()