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