Mercurial > repos > pjbriggs > amplicon_analysis_pipeline
comparison amplicon_analysis_pipeline.py @ 39:d19fca1c009c draft
planemo upload for repository https://github.com/pjbriggs/Amplicon_analysis-galaxy commit c3970ae16b02774148dbf950880f07255734a9a8-dirty
author | pjbriggs |
---|---|
date | Wed, 17 Oct 2018 08:14:27 -0400 |
parents | 7b2c40450792 |
children |
comparison
equal
deleted
inserted
replaced
38:7b2c40450792 | 39:d19fca1c009c |
---|---|
94 fp.write("---- %d: %s\n" % (ix, | 94 fp.write("---- %d: %s\n" % (ix, |
95 os.path.relpath(f,results_dir))) | 95 os.path.relpath(f,results_dir))) |
96 # Close output file | 96 # Close output file |
97 if filen is not None: | 97 if filen is not None: |
98 fp.close() | 98 fp.close() |
99 | |
100 def find_executable(name): | |
101 # Locate executable on PATH | |
102 for p in os.environ['PATH'].split(os.pathsep): | |
103 exe = os.path.normpath( | |
104 os.path.abspath(os.path.join(p,name))) | |
105 if os.path.isfile(exe) and os.access(exe,os.X_OK): | |
106 return exe | |
107 return None | |
108 | 99 |
109 if __name__ == "__main__": | 100 if __name__ == "__main__": |
110 # Command line | 101 # Command line |
111 print "Amplicon analysis: starting" | 102 print "Amplicon analysis: starting" |
112 p = argparse.ArgumentParser() | 103 p = argparse.ArgumentParser() |
147 os.symlink(categories_file,"Categories.txt") | 138 os.symlink(categories_file,"Categories.txt") |
148 print "-- made symlink to Categories.txt" | 139 print "-- made symlink to Categories.txt" |
149 | 140 |
150 # Link to FASTQs and construct Final_name.txt file | 141 # Link to FASTQs and construct Final_name.txt file |
151 sample_names = [] | 142 sample_names = [] |
143 print "-- making Final_name.txt" | |
152 with open("Final_name.txt",'w') as final_name: | 144 with open("Final_name.txt",'w') as final_name: |
153 fastqs = iter(args.fastq_pairs) | 145 fastqs = iter(args.fastq_pairs) |
154 for sample_name,fqr1,fqr2 in zip(fastqs,fastqs,fastqs): | 146 for sample_name,fqr1,fqr2 in zip(fastqs,fastqs,fastqs): |
155 sample_name = clean_up_name(sample_name) | 147 sample_name = clean_up_name(sample_name) |
156 print "Cleaned up sample name: %s" % sample_name | 148 print " %s" % sample_name |
157 r1 = "%s_R1_.fastq" % sample_name | 149 r1 = "%s_R1_.fastq" % sample_name |
158 r2 = "%s_R2_.fastq" % sample_name | 150 r2 = "%s_R2_.fastq" % sample_name |
159 os.symlink(fqr1,r1) | 151 os.symlink(fqr1,r1) |
160 os.symlink(fqr2,r2) | 152 os.symlink(fqr2,r2) |
161 final_name.write("%s\n" % '\t'.join((r1,sample_name))) | 153 final_name.write("%s\n" % '\t'.join((r1,sample_name))) |