Skip to content

Commit

Permalink
Merge pull request #130 from sanger-tol/halfcov
Browse files Browse the repository at this point in the history
fixed findHalfcoverage
  • Loading branch information
DLBPointon authored Aug 10, 2023
2 parents 2a9ff13 + 8675c1d commit a2b1d55
Showing 1 changed file with 3 additions and 5 deletions.
8 changes: 3 additions & 5 deletions bin/findHalfcoverage.py
Original file line number Diff line number Diff line change
Expand Up @@ -72,7 +72,7 @@ def get_cov_peaks(file):
halfPeak = int(peakCov) / 2
qrtPeak = int(peakCov) / 4

print("#Coverage Peak is %s, HalfPeak is %s, QuarterPeak is %s " % (peakCov, halfPeak, qrtPeak))
# print("#Coverage Peak is %s, HalfPeak is %s, QuarterPeak is %s " % (peakCov, halfPeak, qrtPeak))

return (peakCov, halfPeak, qrtPeak)

Expand Down Expand Up @@ -158,15 +158,13 @@ def main():
if coverage > options.covcut:
if totalSize > options.totsize:
print(
"**\t"
+ "\t".join(
"\t".join(
[str(i) for i in [scaffoldName, int(totalSize), int(lowcovSize), "{:.1f}".format(coverage)]]
)
)
else:
print(
"==\t"
+ "\t".join(
"\t".join(
[str(i) for i in [scaffoldName, int(totalSize), int(lowcovSize), "{:.1f}".format(coverage)]]
)
)
Expand Down

0 comments on commit a2b1d55

Please sign in to comment.