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