Mercurial > repos > charles_s_test > seqsero2
diff libs/compare_and_change_two_fastq_id.py @ 0:6895de35a263 draft
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
author | charles_s_test |
---|---|
date | Thu, 19 Oct 2017 18:16:51 -0400 |
parents | |
children | 3d6680af0bec |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/libs/compare_and_change_two_fastq_id.py Thu Oct 19 18:16:51 2017 -0400 @@ -0,0 +1,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) \ No newline at end of file