Mercurial > repos > artbio > probecoverage
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() |