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))) |
