Mercurial > repos > estrain > seqsero_v1
view SeqSero/libs/compare_and_change_two_fastq_id.py @ 4:ab0802d77891 draft default tip
Uploaded
author | estrain |
---|---|
date | Thu, 12 Sep 2019 06:46:00 -0400 |
parents | c577b57b7c74 |
children |
line wrap: on
line source
#!/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)