comparison multicov.py @ 6:7c4feda2d9c7 draft

"planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/probecoverage commit 62a8b073b9ac98b0231641e5266768e7f8b80b89"
author artbio
date Tue, 07 Jan 2020 06:08:56 -0500
parents 35d2db3753d9
children 0adb846ca056
comparison
equal deleted inserted replaced
5:2ec6529e555d 6:7c4feda2d9c7
24 continue 24 continue
25 fields = line[:-1].split('\t') 25 fields = line[:-1].split('\t')
26 chr = fields[0] 26 chr = fields[0]
27 start = fields[1] 27 start = fields[1]
28 end = fields[2] 28 end = fields[2]
29 coverage = bam_object.count_coverage(reference=chr, 29 coverage = bam_object.count_coverage(contig=chr,
30 start=int(start)-1, 30 start=int(start)-1,
31 end=int(end), 31 stop=int(end),
32 quality_threshold=quality) 32 quality_threshold=quality)
33 """ Add the 4 coverage values """ 33 """ Add the 4 coverage values """
34 coverage = [sum(x) for x in zip(*coverage)] 34 coverage = [sum(x) for x in zip(*coverage)]
35 coverage_column.append(numpy.mean(coverage)) 35 coverage_column.append(numpy.mean(coverage))
36 bed_object.close() 36 bed_object.close()