Mercurial > repos > lijing > bubio
changeset 8:e93ae4176127 draft
Uploaded a set of tools
author | lijing |
---|---|
date | Thu, 02 Nov 2017 12:48:54 -0400 |
parents | c267c8a8ec12 |
children | 19d3ec566471 |
files | ._MITObim_1.8.pl ._consense.xml ._dnapars.xml ._interleave-fastqgz-MITOBIM.py ._interleave.xml ._mitobim.xml ._seqboot.xml home/lijing/galaxy/tools/ngs_mapping/interleave-fastqgz-MITOBIM.py |
diffstat | 8 files changed, 0 insertions(+), 45 deletions(-) [+] |
line wrap: on
line diff
--- a/home/lijing/galaxy/tools/ngs_mapping/interleave-fastqgz-MITOBIM.py Thu Nov 02 12:45:03 2017 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,45 +0,0 @@ -#!/usr/bin/python -# encoding:utf8 -# authors: Erik Garrison, Sébastien Boisvert -# modified by github@cypridina on 20151104 to work with MITObim -"""This script takes two fastq or fastq.gz files and interleaves them -Usage: - interleave-fasta fasta_file1 fasta_file2 -""" - -import sys,re - -def interleave(f1, f2): - """Interleaves two (open) fastq files. - """ - while True: - line = f1.readline() - if line.strip() == "": - break - print re.sub(r" 1:N.*", "/1",line.strip()) - - for i in xrange(3): - print re.sub(r" 2:N.*","/2",f1.readline().strip()) - - for i in xrange(4): - print re.sub(r" 2:N.*","/2",f2.readline().strip()) - -if __name__ == '__main__': - try: - file1 = sys.argv[1] - file2 = sys.argv[2] - except: - print __doc__ - sys.exit(1) - - if file1[-2:] == "gz": - import gzip - with gzip.open(file1) as f1: - with gzip.open(file2) as f2: - interleave(f1, f2) - else: - with open(file1) as f1: - with open(file2) as f2: - interleave(f1, f2) - f1.close() - f2.close()