File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -47,7 +47,7 @@ def TCSwriter(bamfile_name, fasta_d):
4747 position = pileupcolumn .pos
4848 tcov = 0 # Workaround
4949 # Define total coverage (AKA "Tcov")
50- #tcov = pileupcolumn.n #TODO - submit bug for wrong counting
50+ # tcov = pileupcolumn.n #TODO - submit bug for wrong counting
5151
5252 # Define base coverages and qualities
5353 for pileupread in pileupcolumn .pileups :
@@ -139,8 +139,9 @@ def covs_and_quals(bases):
139139
140140
141141def QualityCalc (quals ):
142- """Calculate individual bases qualities, just like mira does, as seen here:
143- http://www.freelists.org/post/mira_talk/Quality-Values,4"""
142+ """Calculate and return individual bases qualities, just like mira does,
143+ as seen here:
144+ http://www.freelists.org/post/mira_talk/Quality-Values,4 ."""
144145 quals .sort ()
145146 min1 = quals [0 ]
146147 if min1 > 0 :
You can’t perform that action at this time.
0 commit comments