except: print 'align failed' Message('align failed',email) sys.exit(1) #======== (4) Convert sam to sorted bam ================== try: sort_bams = sam2bam_sort(map_sam,thread) print 'sort bam files succeed' print 'sort_bams is: ',sort_bams except: print 'sort bam files failed' Message('sort bam files failed',email) sys.exit(1) #======== (5) Markduplicates using picard ================ try: dedup_files = markduplicates(picard,sort_bams) print 'mark duplicates succeed' print 'dedup_files is: ',dedup_files remove(sort_bams) except: print 'mark duplicates failed' Message('mark duplicates failed',email) sys.exit(1) #======== 2. Indel realignment ==================================== #======== (6) Create a target list of intervals=========== try: interval = RealignerTargetCreator(gatk,dedup_files,ref_fa,thread,phaseINDEL,gold_indel) print 'RealignerTarget Creator succeed' print 'interval is: ',interval except: print 'RealignerTarget Creator failed'
sys.stdout.write('sort bam failed\n') Message('sort bam failed',email) sys.exit(1) try: group_bams = addReadGroup(picard,sort_bams,read_group) sys.stdout.write('add group succeed\n') sys.stdout.write('group_bams is: {group}\n'.format(group=group_bams)) except: sys.stdout.write('add group failed\n') Message('add group failed',email) sys.exit(1) #======== (2) mark duplicates ============================ try: dedup_bams = markduplicates(picard,group_bams) sys.stdout.write('mark duplicate succeed\n') sys.stdout.write('dedup_bams is: {dedup}\n'.format(dedup=dedup_bams)) remove(group_bams) except: sys.stdout.write('mark duplicate failed\n') Message('mark duplicate failed',email) sys.exit(1) #======== 3. Split 'N' Trim and reassign mapping qualiteies try: split_bams = splitN(gatk,dedup_bams,ref_fa) sys.stdout.write('split N succeed\n') sys.stdout.write('split N is: {N}\n'.format(N=split_bams)) remove(dedup_bams) except: sys.stdout.write('split N failed\n')
except: print 'align failed' Message('align failed', email) sys.exit(1) #======== (4) Convert sam to sorted bam ================== try: sort_bams = sam2bam_sort(map_sam, thread) print 'sort bam files succeed' print 'sort_bams is: ', sort_bams except: print 'sort bam files failed' Message('sort bam files failed', email) sys.exit(1) #======== (5) Markduplicates using picard ================ try: dedup_files = markduplicates(picard, sort_bams) print 'mark duplicates succeed' print 'dedup_files is: ', dedup_files remove(sort_bams) except: print 'mark duplicates failed' Message('mark duplicates failed', email) sys.exit(1) #======== 2. Indel realignment ==================================== #======== (6) Create a target list of intervals=========== try: interval = RealignerTargetCreator(gatk, dedup_files, ref_fa, thread, phaseINDEL, gold_indel) print 'RealignerTarget Creator succeed' print 'interval is: ', interval except: