/usr/share/seqsero/libs/compare_and_change_two_fastq_id.py is in seqsero 1.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 | #!/usr/bin/env python
import os,sys
file1=sys.argv[1]
file2=sys.argv[2]
def compare_and_change_two_fastq_id(file1,file2):
a=os.popen("head "+file1).read().split("\n")
b=os.popen("head "+file2).read().split("\n")
for x in a:
if x.startswith("@"):
a_title=x.split(" ")[0]
for x in b:
if x.startswith("@"):
b_title=x.split(" ")[0]
if a_title==b_title:
pass
else:
print "changing the title of two seperated fastq files..."
print a_title,b_title
os.system("sed "+"-i 's/.1 / /g' "+file1)
print "finished file1"
os.system("sed "+"-i 's/.2 / /g' "+file2)
print "finished file2"
compare_and_change_two_fastq_id(file1,file2)
|