https://github.com/fjruizruano/ngs-protocols
Revision 6c3c041f237f2a54076663591a1779db18409d37 authored by fjruizruano on 04 May 2018, 22:30:12 UTC, committed by fjruizruano on 04 May 2018, 22:30:12 UTC
1 parent 14d1b7e
Tip revision: 6c3c041f237f2a54076663591a1779db18409d37 authored by fjruizruano on 04 May 2018, 22:30:12 UTC
adding annot_to_rexp.py
adding annot_to_rexp.py
Tip revision: 6c3c041
extract_reads_blat.py
#!/usr/bin/python
import sys
try:
file = sys.argv[1]
except:
file = raw_input("Blat output file: ")
data = open(file).readlines()
reads = []
for line in data[5:]:
line = line.split()
reads.append(line[9][:-1])
reads_set = set(reads)
w = open(file+".list", "w")
for r in reads_set:
w.write("%s\n%s\n" % (r+"1", r+"2"))
w.close()
Computing file changes ...