Mercurial > repos > pjbriggs > amplicon_analysis_pipeline
comparison amplicon_analysis_pipeline.py @ 38:7b2c40450792 draft
planemo upload for repository https://github.com/pjbriggs/Amplicon_analysis-galaxy commit 601f7d1e21b25c583cc6e1a563029bab3f226921-dirty
author | pjbriggs |
---|---|
date | Wed, 17 Oct 2018 05:01:57 -0400 |
parents | e832be1635b4 |
children | d19fca1c009c |
comparison
equal
deleted
inserted
replaced
37:e832be1635b4 | 38:7b2c40450792 |
---|---|
58 for line in message.split('\n'): | 58 for line in message.split('\n'): |
59 sys.stderr.write("* %s%s *\n" % (line,' '*(width-len(line)-4))) | 59 sys.stderr.write("* %s%s *\n" % (line,' '*(width-len(line)-4))) |
60 sys.stderr.write("%s\n\n" % ('*'*width)) | 60 sys.stderr.write("%s\n\n" % ('*'*width)) |
61 | 61 |
62 def clean_up_name(sample): | 62 def clean_up_name(sample): |
63 # Remove trailing "_L[0-9]+_001" from Fastq | 63 # Remove extensions and trailing "_L[0-9]+_001" from |
64 # pair names | 64 # Fastq pair names |
65 split_name = sample.split('_') | 65 sample_name = '.'.join(sample.split('.')[:1]) |
66 split_name = sample_name.split('_') | |
66 if split_name[-1] == "001": | 67 if split_name[-1] == "001": |
67 split_name = split_name[:-1] | 68 split_name = split_name[:-1] |
68 if split_name[-1].startswith('L'): | 69 if split_name[-1].startswith('L'): |
69 try: | 70 try: |
70 int(split_name[-1][1:]) | 71 int(split_name[-1][1:]) |
150 sample_names = [] | 151 sample_names = [] |
151 with open("Final_name.txt",'w') as final_name: | 152 with open("Final_name.txt",'w') as final_name: |
152 fastqs = iter(args.fastq_pairs) | 153 fastqs = iter(args.fastq_pairs) |
153 for sample_name,fqr1,fqr2 in zip(fastqs,fastqs,fastqs): | 154 for sample_name,fqr1,fqr2 in zip(fastqs,fastqs,fastqs): |
154 sample_name = clean_up_name(sample_name) | 155 sample_name = clean_up_name(sample_name) |
156 print "Cleaned up sample name: %s" % sample_name | |
155 r1 = "%s_R1_.fastq" % sample_name | 157 r1 = "%s_R1_.fastq" % sample_name |
156 r2 = "%s_R2_.fastq" % sample_name | 158 r2 = "%s_R2_.fastq" % sample_name |
157 os.symlink(fqr1,r1) | 159 os.symlink(fqr1,r1) |
158 os.symlink(fqr2,r2) | 160 os.symlink(fqr2,r2) |
159 final_name.write("%s\n" % '\t'.join((r1,sample_name))) | 161 final_name.write("%s\n" % '\t'.join((r1,sample_name))) |